Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,447,869	Δ1 bp	coding (3758/3780 nt)	tamB →	translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,447,867	0	T	.	100.0%	329.0 / NA	79	coding (3756/3780 nt)	tamB	translocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (42/37);  total (42/37)

CAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTA  >  NC_000913/4447729‑4448000                                                                                                                                           |                                                                                                                                      cAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATATGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                                                                                                       <  1:1546228/139‑1 (MQ=255) cAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                                                                                                       <  1:1172186/139‑1 (MQ=255)       cTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATCTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGctct                                                                                                                                >  1:1138969/1‑139 (MQ=255)              ttGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGtt                                                                                                                         <  2:750366/139‑1 (MQ=255)                    ggTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGtt                                                                                                                   >  1:1264242/1‑139 (MQ=255)                     gTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGttt                                                                                                                  >  1:796704/1‑139 (MQ=255)                      tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGtttt                                                                                                                 >  1:1042825/1‑139 (MQ=255)                        tGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaa                                                                                                               >  1:593111/1‑138 (MQ=255)                        tGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaa                                                                                                               >  2:1366791/1‑138 (MQ=255)                                       gCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttg                                                                                                                                     <  2:1345358/103‑2 (MQ=255)                                       gCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttg                                                                                                                                     >  1:1345358/1‑102 (MQ=255)                                          tGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCt                                                                                             >  1:696082/1‑139 (MQ=255)                                                        tctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGc                                                                               >  2:287756/1‑139 (MQ=255)                                                          tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                             >  2:19739/1‑139 (MQ=255)                                                              gCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                         <  2:696082/139‑1 (MQ=255)                                                                  cacTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                                     >  1:812622/1‑139 (MQ=255)                                                                  cacTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                                     >  2:1568688/1‑139 (MQ=255)                                                                    cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTGACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   >  2:962690/1‑139 (MQ=255)                                                                    cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   >  1:474003/1‑139 (MQ=255)                                                                    cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   <  2:1042825/139‑1 (MQ=255)                                                                    cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   >  1:1332193/1‑139 (MQ=255)                                                                      cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTGACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  1:962690/139‑1 (MQ=255)                                                                      cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 >  2:668711/1‑139 (MQ=255)                                                                      cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  2:1138969/139‑1 (MQ=255)                                                                           tacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                            <  2:1264242/139‑1 (MQ=255)                                                                            acgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaa                                                                                                               <  1:413583/87‑2 (MQ=255)                                                                            acgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaa                                                                                                               >  2:413583/1‑86 (MQ=255)                                                                            acgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAgg                                                                    >  2:1485037/1‑130 (MQ=255)                                                                            acgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAgg                                                                    <  1:1485037/130‑1 (MQ=255)                                                                             cgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                          <  1:287756/139‑1 (MQ=255)                                                                                  tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaac                                                                                                              <  1:535110/82‑3 (MQ=255)                                                                                  tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaac                                                                                                              >  2:535110/1‑80 (MQ=255)                                                                                         aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAAt                                              >  2:173339/1‑139 (MQ=255)                                                                                          tGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATg                                             <  1:1568688/139‑1 (MQ=255)                                                                                            ccTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGATTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                           <  1:1366791/139‑1 (MQ=255)                                                                                            ccTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                           >  1:982151/1‑139 (MQ=255)                                                                                               aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGtt                                                                                                                         <  1:873960/58‑1 (MQ=255)                                                                                               aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGtt                                                                                                                         >  2:873960/1‑58 (MQ=255)                                                                                                aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                             <  1:898084/117‑1 (MQ=255)                                                                                                aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                             >  2:898084/1‑117 (MQ=255)                                                                                                aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                           >  1:1163770/1‑135 (MQ=255)                                                                                                aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                           <  2:1163770/135‑1 (MQ=255)                                                                                                      atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                           <  2:254896/129‑1 (MQ=255)                                                                                                      atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                           >  1:254896/1‑129 (MQ=255)                                                                                                        cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                               <  1:668711/139‑1 (MQ=255)                                                                                                          ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGa                             >  2:269512/1‑139 (MQ=255)                                                                                                           gAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                            >  2:590935/1‑139 (MQ=255)                                                                                                           gAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTAAGGCAGTTTACGCCACAATCAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                            >  1:1037816/1‑139 (MQ=255)                                                                                                            aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                       >  1:330528/1‑95 (MQ=255)                                                                                                            aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                       <  2:330528/95‑1 (MQ=255)                                                                                                             aGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                          <  1:1426941/107‑1 (MQ=255)                                                                                                             aGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                          >  2:1426941/1‑107 (MQ=255)                                                                                                               ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                     >  1:279139/1‑110 (MQ=255)                                                                                                               ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                     <  2:279139/110‑1 (MQ=255)                                                                                                               ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                               >  1:416616/1‑116 (MQ=255)                                                                                                               ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                               <  2:416616/116‑1 (MQ=255)                                                                                                               ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                        <  1:269512/139‑1 (MQ=255)                                                                                                                cGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGtt                                                                                                                         <  2:649948/41‑1 (MQ=255)                                                                                                                cGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGtt                                                                                                                         >  1:649948/1‑41 (MQ=255)                                                                                                                  tgtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaac                                                                                                              <  2:595586/50‑3 (MQ=255)                                                                                                                  tgtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaac                                                                                                              >  1:595586/1‑48 (MQ=255)                                                                                                                  tgtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                     <  2:796704/139‑1 (MQ=255)                                                                                                                   gtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                    <  1:19739/139‑1 (MQ=255)                                                                                                                     cTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCg                  <  2:1037816/139‑1 (MQ=255)                                                                                                                       ggTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt                <  2:593111/139‑1 (MQ=255)                                                                                                                           tAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATCACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt            <  2:982151/139‑1 (MQ=255)                                                                                                                            aGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                           <  2:1331914/43‑1 (MQ=255)                                                                                                                            aGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                           >  1:1331914/1‑43 (MQ=255)                                                                                                                              aCCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAcc         <  1:428984/139‑1 (MQ=255)                                                                                                                                  ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                             <  2:1289508/83‑1 (MQ=255)                                                                                                                                  ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                             >  1:1289508/1‑83 (MQ=255)                                                                                                                                  ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                               >  1:765203/1‑113 (MQ=255)                                                                                                                                  ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                               <  2:765203/113‑1 (MQ=255)                                                                                                                                  ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                    <  2:61485/124‑1 (MQ=255)                                                                                                                                  ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                    >  1:61485/1‑124 (MQ=255)                                                                                                                                  ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg       <  1:888188/137‑1 (MQ=255)                                                                                                                                  ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg       >  2:888188/1‑137 (MQ=255)                                                                                                                                   gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 >  1:968839/1‑78 (MQ=255)                                                                                                                                   gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  2:968839/78‑1 (MQ=255)                                                                                                                                   gCACTGGA‑TTGCTCTATCAGATCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTg    >  2:1093748/1‑139 (MQ=255)                                                                                                                                     aCTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                         <  1:370346/68‑1 (MQ=255)                                                                                                                                     aCTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                         >  2:370346/1‑68 (MQ=255)                                                                                                                                     aCTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta  >  2:620761/1‑139 (MQ=255)                                                                                                                                           |                                                                                                                                      CAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTA  >  NC_000913/4447729‑4448000

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑