Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,445,428	A→G	Q439Q (CAA→CAG)	tamB →	translocation and assembly module subunit TamB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,445,428	0	A	G	100.0%	245.8 / NA	67	Q439Q (CAA→CAG)	tamB	translocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/36);  total (31/36)

GCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATGGA  >  NC_000913/4445291‑4445562                                                                                                                                          |                                                                                                                                       gCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAGCTCAAGGCGTTGCTCGACTGgcagca                                                                                                                                      <  1:975980/140‑1 (MQ=255) gCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGgcagca                                                                                                                                      >  1:521223/1‑140 (MQ=255)     tGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGcc                                                                                                                                  >  1:800834/1‑140 (MQ=255)     tGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGcc                                                                                                                                  >  2:857554/1‑140 (MQ=255)            aCAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGtt                                                                                                                           >  1:814651/1‑140 (MQ=255)                gagaTCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGgcg                                                                                                                       <  1:156047/140‑1 (MQ=255)                gagaTCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGgcg                                                                                                                       <  2:814651/140‑1 (MQ=255)                 agaTCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGgcgc                                                                                                                      >  2:507650/1‑140 (MQ=255)                  gaTCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGcgcg                                                                                                                     <  1:507650/140‑1 (MQ=255)                   aTCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCgg                                                                                                                    >  1:962302/1‑140 (MQ=255)                    tCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGGGGt                                                                                                                   <  2:587026/140‑1 (MQ=255)                     cccgcccgCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTg                                                                                                                  >  2:169682/1‑140 (MQ=255)                     cccgcccgCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTg                                                                                                                  >  2:198883/1‑140 (MQ=255)                        gcccgcGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGc                                                                                                               <  2:931448/140‑1 (MQ=255)                        gcccgcGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGc                                                                                                               >  1:70523/1‑140 (MQ=255)                          ccgcGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGgcg                                                                                                                       >  1:896194/1‑130 (MQ=255)                          ccgcGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGgcg                                                                                                                       <  2:896194/130‑1 (MQ=255)                           cgcgACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTaa                                                                                                            <  2:395656/140‑1 (MQ=255)                               aCTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCt                                                                                                        >  1:301434/1‑140 (MQ=255)                                   ttACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAAc                                                                                                    <  2:585254/140‑1 (MQ=255)                                   ttACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAAc                                                                                                    <  1:198883/140‑1 (MQ=255)                                    tACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACg                                                                                                   >  2:725882/1‑140 (MQ=255)                                         tCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCAtt                                                                                              <  2:521223/140‑1 (MQ=255)                                          cGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTa                                                                                                      >  1:821048/1‑131 (MQ=255)                                          cGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTa                                                                                                      <  2:821048/131‑1 (MQ=255)                                          cGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTa                                                                                             >  1:67718/1‑140 (MQ=255)                                           gACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTaa                                                                                            <  1:574360/140‑1 (MQ=255)                                             cGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAAc                                                                                                    >  2:392430/1‑130 (MQ=255)                                             cGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAAc                                                                                                    <  1:392430/130‑1 (MQ=255)                                             cGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGGATTAaca                                                                                          <  2:70523/140‑1 (MQ=255)                                             cGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAaca                                                                                          >  2:42502/1‑140 (MQ=255)                                             cGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAaca                                                                                          <  1:857554/140‑1 (MQ=255)                                             cGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAaca                                                                                          <  2:706124/140‑1 (MQ=255)                                             cGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAaca                                                                                          <  2:106845/140‑1 (MQ=255)                                                cAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCg                                                                                       <  2:962302/140‑1 (MQ=255)                                                           aaCAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTc                                                                            <  2:362753/140‑1 (MQ=255)                                                            aCAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTcc                                                                           <  2:432646/140‑1 (MQ=255)                                                               gcagGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAag                                                                                 >  1:1038519/1‑131 (MQ=255)                                                               gcagGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAag                                                                                 <  2:1038519/131‑1 (MQ=255)                                                               gcagGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCgg                                                                        >  2:819412/1‑140 (MQ=255)                                                                  ggTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACt                                                                     >  1:784457/1‑140 (MQ=255)                                                                   gTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTg                                                                    <  2:67718/140‑1 (MQ=255)                                                                        tctcGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGt                                                               >  1:516765/1‑140 (MQ=255)                                                                           cGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGTCAAAGAGTTCCCGGACTGGCCGTCGa                                                            <  2:799399/140‑1 (MQ=255)                                                                            gACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGaa                                                           <  2:99177/140‑1 (MQ=255)                                                                            gACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGaa                                                           <  1:601740/140‑1 (MQ=255)                                                                             aCAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGt                                                               >  2:194334/1‑135 (MQ=255)                                                                             aCAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGt                                                               <  1:194334/135‑1 (MQ=255)                                                                                      ccGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGc                                                 >  1:845249/1‑140 (MQ=255)                                                                                       cGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCt                                                <  1:826622/140‑1 (MQ=255)                                                                                               cgcTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTaaa                                        >  2:154461/1‑140 (MQ=255)                                                                                                    gAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCcg                                   >  2:153311/1‑140 (MQ=255)                                                                                                    gAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCcg                                   <  1:54245/140‑1 (MQ=255)                                                                                                      aGGGAAAACTGAACTCAAGGCGTTGCTCGGCTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCcgcg                                 <  1:935355/140‑1 (MQ=255)                                                                                                      aGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCcgcg                                 <  1:111886/140‑1 (MQ=255)                                                                                                      aGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCcgcg                                 <  2:49707/140‑1 (MQ=255)                                                                                                         gAAAACTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTa                              >  2:642319/1‑140 (MQ=255)                                                                                                           aaaCTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGc                            <  2:128427/140‑1 (MQ=255)                                                                                                           aaaCTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGc                            <  1:169682/140‑1 (MQ=255)                                                                                                           aaaCTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGc                            >  1:166871/1‑140 (MQ=255)                                                                                                           aaaCTGAACTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGc                            <  1:154461/140‑1 (MQ=255)                                                                                                                   cTCAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTAcgg                    >  1:636234/1‑140 (MQ=255)                                                                                                                     cAAGGCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTAcggcg                  <  2:166871/140‑1 (MQ=255)                                                                                                                        ggCGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACggcggcg               >  1:854772/1‑139 (MQ=255)                                                                                                                          cGTTGCTCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACggcggcgcc             <  2:800834/140‑4 (MQ=255)                                                                                                                                tCGACTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAg       >  2:573505/1‑140 (MQ=255)                                                                                                                                   aCTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATg    >  2:803484/1‑140 (MQ=255)                                                                                                                                   aCTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATg    >  1:795187/1‑140 (MQ=255)                                                                                                                                    cTGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATgg   >  2:614406/1‑140 (MQ=255)                                                                                                                                     tGGCAGCAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGa  <  2:845249/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       GCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATGGA  >  NC_000913/4445291‑4445562

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