Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,194,922 G→A T10T (ACC→ACT)  rluE ← 23S rRNA pseudouridine(2457) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,194,9220GA97.9% 142.6 / ‑6.1 48T10T (ACC→ACT) rluE23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (25/22);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAGTTGTTCCGGCACAAGGCAAGAAGATCAC  >  NC_000913/1194676‑1195166
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
gcCGTTATTTGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAAACGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGCAGAACCTCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTAt                                                                                                                                                                                                                                                  >  1:628129/1‑251 (MQ=255)
 ccGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTCCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTAtt                                                                                                                                                                                                                                                 >  1:325094/1‑251 (MQ=255)
        tGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGACCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTGTCaga                                                                                                                                                                                                                                           <  1:543372/250‑1 (MQ=255)
        tGGTCAGCACCAGCAACCCTACGCTATCGCGGTCAAGGCGACCTGCTGCATAAACGCCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTGCTAAAAGAAGTTTTTTGCATAGTATTTTCaga                                                                                                                                                                                                                                           <  1:411666/250‑1 (MQ=255)
           tCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGAGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGAtt                                                                                                                                                                                                                                       >  2:669276/1‑251 (MQ=255)
            cAGCACCATCAACCCTACGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAGGGACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGCTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACATGGTGTTTTATAAACGAAGTTTTGTGCATAGTATTTTCAGCGAttt                                                                                                                                                                                                                                      >  2:529969/1‑250 (MQ=255)
                                          aaGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATCGCCTAATAAcgcg                                                                                                                                                                                                        >  1:330074/1‑251 (MQ=255)
                                                  ctgctgCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTACCCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCAt                                                                                                                                                                                                                         >  2:375726/1‑226 (MQ=255)
                                                  ctgctgCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCAt                                                                                                                                                                                                                         <  1:375726/226‑1 (MQ=255)
                                                         aTAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACgg                                                                                                                                                                                         >  1:227003/1‑251 (MQ=255)
                                                                   tGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAAAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCCGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAAt                                                                                                                                                                                 >  1:51951/1‑249 (MQ=255)
                                                                      tccGGGATGAAGTACTTTAATGTTGTGTGGCAGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGGAGGGTTGAGTGAAAAGGATCACACGCGTGGGCTGGTTTGCAGGTGTACGCCTGGTAGAACGTTGCGCGCTGAATCTCTTAACCTGGTGATTTCTTAAAGAAGTTTTTTGCGTAGTATTTGCAGAGATTATGAATTTCCGCAGTATAGCCTAATAACGCGCAGTTTTAATGACGGccaatactaggg                                                                                                                                                                             <  2:616796/249‑12 (MQ=255)
                                                                                   cTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGc                                                                                                                                                               >  1:70087/1‑251 (MQ=255)
                                                                                     tttAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAg                                                                                                                                                                <  1:119531/248‑1 (MQ=255)
                                                                                     tttAATGTGTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACTCGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGcc                                                                                                                                                              <  2:102062/250‑1 (MQ=255)
                                                                                         aTGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTa                                                                                                                                                         <  1:305877/251‑1 (MQ=255)
                                                                                              ttgtgTCCGGCTTCATCGGTGTGC‑CGCGGAAGAACATCGTAGGGGTTATTGAACAGGATCACACGCGTGTGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACaaa                                                                                                                                                     <  2:557647/246‑1 (MQ=255)
                                                                                               tgtgTCCTGCTTCATCG‑TGTACCTGCTGAAGAACATCGGAGGGTTTATTGAACATGATCCCACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTCGTACGTTGCTAGCTGAATCTCTTATTCTGGTGATTTCTAAAAGAAGGTTTTTGCATAGTATTTTCAGAGATTATGAATGGCCGCATTAGAGCCTAATAACGCGCATCTAACATGCCGGCAAACAATAGGGTAGTATTTACAAGCCAATTACAAAt                                                                                                                                                    <  2:320449/247‑1 (MQ=255)
                                                                                                         ttCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACaaa                                                                                                                                          >  1:359362/1‑250 (MQ=255)
                                                                                                          tCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACaaaa                                                                                                                                         >  1:80951/1‑250 (MQ=255)
                                                                                                           cATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACaaaaa                                                                                                                                        >  1:239775/1‑250 (MQ=255)
                                                                                                            aTCGGTGAAC‑TGCGGCAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGAGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGTTTATGAATTACCGCATTATAACCTAATAACGAGCATCTTTCATGACGGCAAACACTCGGGAAGTATTGACAAGCAAATTACAAATCATTAACAAAAAAt                                                                                                                                      >  2:372236/1‑251 (MQ=255)
                                                                                                               ggTGAAA‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGCTCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGc                                                                                                                                   >  2:473532/1‑251 (MQ=255)
                                                                                                                  gAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGctct                                                                                                                                <  1:671599/251‑1 (MQ=255)
                                                                                                                   aaC‑TGCGGAAGATGATCGTAGTGTTCATGGAAAAGGATCACACGCGTGTGCTGTTTTTCTGGTTTTCGCCTGGTAGTAGGTTCGGAGATGAATCGCTTAACCTGGTTATTGCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTACTAACACGCATCTTTCATGAAGGCAAAAAATAGGGTAGTATTGACAAACAAATTACAAATCATTAACAAAAAAGTGCTCTa                                                                                                                               <  2:478569/251‑1 (MQ=255)
                                                                                                                           gaagaaCATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATc                                                                                                                        <  1:669276/251‑1 (MQ=255)
                                                                                                                              gaaaaTCGGAGGGTTTATTGAACAGGAGCACTCGCGTGGGCTGGTTGTCAGTTTTACGCCTGGTAGAACTGTGCGATCTGAATCGCTTAACCTGTTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTAGTGACAAGCTAATTACAAATCATTAACAAAAAATTGCTCTAAAGTATCCGt                                                                                                                     <  2:76699/247‑1 (MQ=255)
                                                                                                                                                tGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAACAAATTGCTCTAAAGCATACGTATCGCAGGACGCAAACGc                                                                                                   >  2:658539/1‑251 (MQ=255)
                                                                                                                                                          cccacgctggggCTGGTTTTCAGGTGTACGCCGGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAAc                                                                                          <  2:679739/241‑1 (MQ=255)
                                                                                                                                                           acacGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTCGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTACCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTCCAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACgg                                                                                        >  2:615101/1‑250 (MQ=255)
                                                                                                                                                                     gatGGTTTTCAGGTTTACGCCTGGTAGAACTTTGCGAGTTTAATCGCTTAACCTTGTGATTTCTAAAAGAAGTTTTTTGCATAGTATGTTCAGAGATTATTAATTGCCGCATTATAGCCTAAGAACGCGCATCTTTTATGACGGCAAACAATAGGGTAGTATTGACAAGCAAATTACAAATCATTAACAAAAAAGTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGGGGTGGAAGaa                                                                              <  1:392070/249‑2 (MQ=255)
                                                                                                                                                                        ggTTGTCAGGATTACGCCTGGTAGAACGGTGCGAGCGGAATCACCTAACCTGGTGTTTTCTAAAAGAAGTTTGTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGACAATTACAAATCATTAACAAAAAATTGCTCCAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGaa                                                                           <  2:546498/251‑2 (MQ=255)
                                                                                                                                                                                                   gTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGt                                                <  1:538386/251‑1 (MQ=255)
                                                                                                                                                                                                   gTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAggg                                                >  1:444095/1‑250 (MQ=255)
                                                                                                                                                                                                      gCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGAATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGa                                              <  2:359362/250‑1 (MQ=255)
                                                                                                                                                                                                              aaTCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagt                                     <  1:615101/251‑1 (MQ=255)
                                                                                                                                                                                                                  gCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGAGGCAGACGAGCAAACCAGTAGCGCTAGAAGGAGAGGTGAATGGaaagtaa                                   >  2:153976/1‑249 (MQ=255)
                                                                                                                                                                                                                    ttAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagtaaagt                                >  1:52734/1‑250 (MQ=255)
                                                                                                                                                                                                                      aaCCTGGTGATTTCTAAAAGCAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGCCAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGAATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAAGTGAATGGgaaataaagaagg                             >  2:635774/1‑247 (MQ=255)
                                                                                                                                                                                                                      aaCCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAagtagt                             <  1:372236/251‑1 (MQ=255)
                                                                                                                                                                                                                           ggTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAGTTGTTc                        >  1:209193/1‑251 (MQ=255)
                                                                                                                                                                                                                            gTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTACCGCATTATAGCTTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATAAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATTCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTGGAAGGGGATGAGAGTGGAAAGTAAAGTAGTTGTTcc                       >  2:351810/1‑251 (MQ=255)
                                                                                                                                                                                                                             tGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGATAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGCCGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAGTTGTTCCg                      >  1:290906/1‑251 (MQ=255)
                                                                                                                                                                                                                                tttCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACACTAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAg                                                                                                                            >  2:473219/1‑146 (MQ=255)
                                                                                                                                                                                                                                tttCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAg                                                                                                                            <  1:473219/146‑1 (MQ=255)
                                                                                                                                                                                                                                      aaaGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAGTTGTTCCGGCACAAGGc             <  1:473532/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ttttttACATAGTATTTTCAGAGATTATGAATTACCGCATTCTAGCCTAATAAAGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGAAAAGCCCATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCCGAAGCGCAGGAATGAGAGGGGAAGGGAAAGTAAAGAAGTTGTTCCGGcaccagcaacga         >  2:432515/1‑239 (MQ=255)
                                                                                                                                                                                                                                                  tGCATAGTATTTTCAGATATTATGAATTGCCGTATTATAGCCTAAGAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGTCGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAGTTGTTCCGGTACAATGCAAGACGATCAc  <  2:317488/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAAC‑TGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAGTTGTTCCGGCACAAGGCAAGAAGATCAC  >  NC_000913/1194676‑1195166

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: ‑