Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,771,257 A→G S41S (AGT→AGC yibH ← inner membrane protein YibH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,771,2570AG80.6% 137.4 / 22.1 73S41S (AGT→AGCyibHinner membrane protein YibH
Reads supporting (aligned to +/- strand):  ref base A (2/12);  new base G (26/32);  total (28/45)
Fisher's exact test for biased strand distribution p-value = 6.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

GCACGCAGCGTCTTTATATTATGCGTCGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGAAATTCGGTCAGAAAGTGCGGCGATTTCCGGGTCAGTCGCTTTCTCCTGCGACGC  >  NC_000913/3771015‑3771501
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
gcacgcgGCGTCTTTATATTATGCGTCGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAg                                                                                                                                                                                                                                              <  2:222523/251‑1 (MQ=255)
 cacgcaGCGTCTTTATATTATGCGTCGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCAcc                                                                                                                                                                                                                                                >  2:292942/1‑248 (MQ=255)
       gCGTCTTTATATTATGCGTCGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACg                                                                                                                                                                                                                                       <  1:412114/251‑1 (MQ=255)
         gTCTTTATATTATGCGTCGCCGTCATCAGGTCAGCCTTAAGGCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCAGCAACAAAATCAAACCGCTCTCCAGAAACACGcc                                                                                                                                                                                                                                     <  2:430218/251‑1 (MQ=255)
          tCTTTATATTATGCGTCGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGccc                                                                                                                                                                                                                                    >  1:561852/1‑250 (MQ=255)
             ttATATTATGCGTCGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCTCTCACCAGAAACACGCCTcc                                                                                                                                                                                                                                  <  2:293827/250‑1 (MQ=255)
              tataTTATGCGTCGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTccc                                                                                                                                                                                                                                 <  1:120284/250‑1 (MQ=255)
                 attatGCGTCGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAgc                                                                                                                                                                                                                              <  1:374511/250‑1 (MQ=255)
                    atGCGTCGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGcgcc                                                                                                                                                                                                                           <  1:613939/250‑1 (MQ=255)
                    atGCGTCGCCGTCAGCAGGTCAGACTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCCCCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGGGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGATCACCAGAAACACGCCTCCCATcgcc                                                                                                                                                                                                                           <  2:142280/250‑1 (MQ=255)
                      gCGTCGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCgccgcc                                                                                                                                                                                                                        <  1:199921/251‑1 (MQ=255)
                             cGTCATCAGGTCAGCCTGAAGTCTGTCACCTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGACTAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGccc                                                                                                                                                                                                                 >  2:55645/1‑250 (MQ=255)
                             cGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTTGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCa                                                                                                                                                                                                                 <  1:20191/251‑1 (MQ=255)
                               tcatcaGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCCGCGCCGCCGTCGcccgc                                                                                                                                                                                                               >  1:95011/1‑248 (MQ=255)
                                    aGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTGTTGAATAAGCTGAGTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCGCCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGTTGGTGGTAGGTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCGAGCGTCCa                                                                                                                                                                                                          <  2:288845/251‑1 (MQ=255)
                                              gAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTa                                                                                                                                                                                                  <  2:482058/249‑1 (MQ=255)
                                              gAAGTCTGGCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTGTGGAATAAGCTGATTATTCTTGGCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTAc                                                                                                                                                                                                 <  2:570034/250‑1 (MQ=255)
                                               aaGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCg                                                                                                                                                                                               >  1:16388/1‑251 (MQ=255)
                                               aaGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCg                                                                                                                                                                                               >  1:172010/1‑251 (MQ=255)
                                                     gTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCGTTTGAATAAGCTGATTATTCTTTTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGTCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCAGCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAAt                                                                                                                                                                                           <  1:575433/249‑1 (MQ=255)
                                                       cAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCg                                                                                                                                                                                        <  1:71538/250‑1 (MQ=255)
                                                                ccTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGTGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTa                                                                                                                                                                              <  2:446930/251‑1 (MQ=255)
                                                                   ggTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTGCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAg                                                                                                                                                                           <  1:412491/251‑1 (MQ=255)
                                                                       aCGAACCGGTGCGAGCTTAAAAAGCACCTCGCCCGGTTGAATAAGCTGATTATTCTTGTAAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACGGCCTGTTGCGCGGTAAAAGTGTAAGGGTTGTTTTAGATCAGCAACAAAATCAAACCGCTCACCAGAAACACGCCACCCAGAGCCGCCGTCGCCAGCGTCCACTGATGGACCGGAATGCGGAAGAATTGAAAGTCCg                                                                                                                                                                       <  2:621819/251‑1 (MQ=255)
                                                                        cGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGc                                                                                                                                                                      >  1:588449/1‑251 (MQ=255)
                                                                         gAACCGGGTCGAGCTTAAAAAGCCCCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGcc                                                                                                                                                                     <  1:12492/251‑1 (MQ=255)
                                                                         gAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGcc                                                                                                                                                                     <  1:20116/251‑1 (MQ=255)
                                                                          aaCCGTGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCAGCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGccc                                                                                                                                                                    <  2:294716/251‑1 (MQ=255)
                                                                           aCCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTGCAGTAACAATTCCCGTCACCTGTGGCGGGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCATCGTCCACTGATTTACCGTAATGCGGAAGATTTTAAAGACCGCCCa                                                                                                                                                                   <  2:620126/251‑1 (MQ=255)
                                                                                        tAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAAAAAAATCACACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGACGCCCGCGCCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCatgccctcg                                                                                                                                                      >  2:334514/1‑243 (MQ=255)
                                                                                         aaaaaGCACCTCGTCCTTTTGAATAATCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCCCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACg                                                                                                                                                      <  2:461272/250‑1 (MQ=255)
                                                                                                  ctagCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGTCGTGCTAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAAACGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGAAAGTTaaa                                                                                                                                             >  2:193364/4‑250 (MQ=255)
                                                                                                               aTAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCCGCGCCACGTAAGTTAAAAcaatcaat                                                                                                                                    >  2:71205/1‑246 (MQ=255)
                                                                                                                  aGCTGAGTATTCTTGGCAGTGACTTCAGTAACAATTCCCGTCACCTTTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGTTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACCCGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCAt                                                                                                                            <  1:5449/251‑1 (MQ=255)
                                                                                                                   gCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATa                                                                                                                           <  1:513891/251‑1 (MQ=255)
                                                                                                                   gCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGGGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAAGCAAACCACTCCCCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGGAAGTTAAAACAATCAATAGATCCATa                                                                                                                           <  2:343757/251‑1 (MQ=255)
                                                                                                                    cTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATa                                                                                                                           >  2:191629/1‑250 (MQ=255)
                                                                                                                      gATTATTCTTGTCATTGAATTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCTCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCTGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATaa                                                                                                                          <  1:86283/249‑1 (MQ=255)
                                                                                                                      gATTATTCTTGTCAGGGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAAt                                                                                                                         <  1:246744/250‑1 (MQ=255)
                                                                                                                       attattCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAAtt                                                                                                                        >  1:279990/1‑250 (MQ=255)
                                                                                                                         tattCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACGGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACGGATTTACCGGAATGCGGAAGATTTTAAAGACCTCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATATATCCATAATTAg                                                                                                                      <  2:248506/250‑1 (MQ=255)
                                                                                                                            tCTTGTCAGTGACTGCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGGGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCGCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTGACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAAGAGATCCATAATTAGAtc                                                                                                                   <  2:611969/250‑1 (MQ=255)
                                                                                                                              ttGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTcc                                                                                                                >  2:144246/1‑251 (MQ=255)
                                                                                                                               tGTCAGGGACTTCAGTAACAATTCCCGTCACCTGGGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGGTAAAACAATCAATAGAGCCATAATTAGAtctc                                                                                                                 <  2:78781/249‑1 (MQ=255)
                                                                                                                                 tCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAATTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAAAAATCAATAGATCCATAATTAGATCTCCggg                                                                                                             >  1:231643/1‑251 (MQ=255)
                                                                                                                                  cAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACAAGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGCTTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCggg                                                                                                             >  2:146158/1‑250 (MQ=255)
                                                                                                                                   attGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTTTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACTCGCCTCCCATCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGaa                                                                                                           <  1:354643/249‑1 (MQ=255)
                                                                                                                                    gTGACTTCAGTAACAATGCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAAc                                                                                                          <  2:304369/251‑1 (MQ=255)
                                                                                                                                            aGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAACCCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAgc                                                                                                                                                            <  1:389343/193‑1 (MQ=255)
                                                                                                                                               aaCAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCGTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCAGCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGAGTTTAAATACCGCCCACGCCAGCTCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGAGCTCCTTGAACGTGGAATagtc                                                                                               <  2:609801/251‑1 (MQ=255)
                                                                                                                                                  aaTTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATagtcagt                                                                                            >  1:60010/1‑251 (MQ=255)
                                                                                                                                                      cccGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTCCCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCttt                                                                                        <  1:287061/251‑1 (MQ=255)
                                                                                                                                                      cccGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGATAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCttt                                                                                        >  1:543143/1‑251 (MQ=255)
                                                                                                                                                       ccGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCtttt                                                                                       <  1:68164/251‑1 (MQ=255)
                                                                                                                                                        cGTCACCTGTGGCGTGATAGGGATCGCTATCCCTGCCTTTTGCGCGGTAAAAGCGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCttttt                                                                                      >  1:7836/1‑251 (MQ=255)
                                                                                                                                                                  ggCGTTATAGGGATCGCTATCACGGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCCCCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCggcggc                                                                             <  1:114789/250‑1 (MQ=255)
                                                                                                                                                                   gCGTGATAGGGATCGCTACCACTGCCTTTGGCGCGGTAAAAGTTTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGTAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAACTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTACTCGGCGGCGa                                                                           <  2:64171/251‑1 (MQ=255)
                                                                                                                                                                          aGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGAt                                                                     <  1:278859/250‑1 (MQ=255)
                                                                                                                                                                                 gCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCCCCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAAGAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGc                                                             <  2:94341/251‑1 (MQ=255)
                                                                                                                                                                                       aCTGCCTTTTGCGCGGTAATAGGGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACTTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGGGGAATAGTCATTCTTTTTCTCGGCGGCGAGGTGAGGCTCCAGCCGGGaa                                                       <  2:112335/251‑1 (MQ=255)
                                                                                                                                                                                                      gTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGAAATTCGGTCAGAAAGt                                        >  1:562166/1‑251 (MQ=255)
                                                                                                                                                                                                        aaaaGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACCCGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGAAATTCGGTCAGAAAGTg                                       <  1:507605/250‑1 (MQ=255)
                                                                                                                                                                                                        aaaaGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGAAATTCGGTCAGACAGTgg                                      >  1:409081/1‑250 (MQ=255)
                                                                                                                                                                                                        aaaaGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGAAATTCGGTCAGAAAGTgc                                      >  1:9496/1‑251 (MQ=255)
                                                                                                                                                                                                             tgtcaGGGTGGTTGAAGTTCATCAACCAACTCAAACCGCTCACCAGAAACACGCTTCCCAGAGCCGCCGTCGCCAGCGTCCACTGCTTTCCCGGCCTGCGGCCGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGCTCTCCGGGAACGTGGACTAGTCAGTCTTTTTCTCGGCGGCGCGTTGATGCTCCAGCCGGGACATTCGGTACGAAAGTgcggcg                                  >  2:564566/5‑250 (MQ=255)
                                                                                                                                                                                                               tAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCCGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGAAATTCGGTCAGAAAGTGAGGCGAtt                               >  2:520804/1‑251 (MQ=255)
                                                                                                                                                                                                                    gtggtTGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACTCGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGAAATTCGGTCAGAAAGTGCGGCGATTTCCgg                          <  2:60010/251‑1 (MQ=255)
                                                                                                                                                                                                                       gtagTAGTTCATCAACAAAATCAAACCACTCCCCAGACACACGCCTCCCAGAGACGCCGTCGCCAGCGTCCACTGCTTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCCCGTAAGTTAAAAAAATCAATAGATCCATAATTAGATCTCCGGGAAAGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGCTGCTCCAGCCGGGAACTTCGGTCAGAACGAGCGGCGATTTCCGGGTc                       >  2:418880/4‑251 (MQ=255)
                                                                                                                                                                                                                        ttGTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGACTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGAAATTCGGTCAGAAAGTGCGGCGATTTCCGGGTc                       >  1:226446/1‑250 (MQ=255)
                                                                                                                                                                                                                          gTAGTTCATCAACAAAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGAAATTCGGTCAGAAAGTGCGGCGATTTCCGGGTCAg                     <  1:71205/250‑1 (MQ=255)
                                                                                                                                                                                                                                      caaAATCAAACCGCTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGACATTCGGTCCGAAAGTGCGGCGATTTCCGGGTCAGTCGCTTTCTCCtt        >  2:310948/1‑250 (MQ=255)
                                                                                                                                                                                                                                         aaTCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGAAATTCGGTCAGAAAGTGCGGCGATTTCCGGGTCAGTCGCTTTCTCCTGCg      >  1:313815/1‑250 (MQ=255)
                                                                                                                                                                                                                                             aaaCCGCTCACCAGAAACACGACTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGGGTTGCTGCTCCAGCCGGGAAATTCGGTCAGACAGTGCGGCGATTTCCGGGTCAGTCGCTTTCTCCTGCGACGc  >  2:277876/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
GCACGCAGCGTCTTTATATTATGCGTCGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCACGCCAGCGCCACGTAAGTTAAAACAATCAATAGATCCATAATTAGATCTCCGGGAACGTGGAATAGTCAGTCTTTTTCTCGGCGGCGAGTTGATGCTCCAGCCGGGAAATTCGGTCAGAAAGTGCGGCGATTTCCGGGTCAGTCGCTTTCTCCTGCGACGC  >  NC_000913/3771015‑3771501

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: