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→ACTrluE23S 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

GCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAGTTGTTCCGGCACAAGGCAAGAAGATCAC  >  NC_000913/1194676‑1195166
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
gcCGTTATTTGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAAACGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGCAGAACCTCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTAt                                                                                                                                                                                                                                                  >  1:628129/1‑251 (MQ=255)
 ccGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTCCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTAtt                                                                                                                                                                                                                                                 >  1:325094/1‑251 (MQ=255)
        tGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGACCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTGTCaga                                                                                                                                                                                                                                           <  1:543372/250‑1 (MQ=255)
        tGGTCAGCACCAGCAACCCTACGCTATCGCGGTCAAGGCGACCTGCTGCATAAACGCCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTGCTAAAAGAAGTTTTTTGCATAGTATTTTCaga                                                                                                                                                                                                                                           <  1:411666/250‑1 (MQ=255)
           tCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGAGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGAtt                                                                                                                                                                                                                                       >  2:669276/1‑251 (MQ=255)
            cAGCACCATCAACCCTACGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAGGGACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGCTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACATGGTGTTTTATAAACGAAGTTTTGTGCATAGTATTTTCAGCGAttt                                                                                                                                                                                                                                      >  2:529969/1‑250 (MQ=255)
                                          aaGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATCGCCTAATAAcgcg                                                                                                                                                                                                        >  1:330074/1‑251 (MQ=255)
                                                  ctgctgCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTACCCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCAt                                                                                                                                                                                                                         >  2:375726/1‑226 (MQ=255)
                                                  ctgctgCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCAt                                                                                                                                                                                                                         <  1:375726/226‑1 (MQ=255)
                                                         aTAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACgg                                                                                                                                                                                         >  1:227003/1‑251 (MQ=255)
                                                                   tGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAAAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCCGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAAt                                                                                                                                                                                 >  1:51951/1‑249 (MQ=255)
                                                                      tccGGGATGAAGTACTTTAATGTTGTGTGGCAGGCTTCATCGGTGAACTGCGGAAGAACATCGGAGGGTTGAGTGAAAAGGATCACACGCGTGGGCTGGTTTGCAGGTGTACGCCTGGTAGAACGTTGCGCGCTGAATCTCTTAACCTGGTGATTTCTTAAAGAAGTTTTTTGCGTAGTATTTGCAGAGATTATGAATTTCCGCAGTATAGCCTAATAACGCGCAGTTTTAATGACGGccaatactaggg                                                                                                                                                                             <  2:616796/249‑12 (MQ=255)
                                                                                   cTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGc                                                                                                                                                               >  1:70087/1‑251 (MQ=255)
                                                                                     tttAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAg                                                                                                                                                                <  1:119531/248‑1 (MQ=255)
                                                                                     tttAATGTGTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACTCGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGcc                                                                                                                                                              <  2:102062/250‑1 (MQ=255)
                                                                                         aTGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTa                                                                                                                                                         <  1:305877/251‑1 (MQ=255)
                                                                                              ttgtgTCCGGCTTCATCGGTGTGCCGCGGAAGAACATCGTAGGGGTTATTGAACAGGATCACACGCGTGTGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACaaa                                                                                                                                                     <  2:557647/246‑1 (MQ=255)
                                                                                               tgtgTCCTGCTTCATCGTGTACCTGCTGAAGAACATCGGAGGGTTTATTGAACATGATCCCACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTCGTACGTTGCTAGCTGAATCTCTTATTCTGGTGATTTCTAAAAGAAGGTTTTTGCATAGTATTTTCAGAGATTATGAATGGCCGCATTAGAGCCTAATAACGCGCATCTAACATGCCGGCAAACAATAGGGTAGTATTTACAAGCCAATTACAAAt                                                                                                                                                    <  2:320449/247‑1 (MQ=255)
                                                                                                         ttCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACaaa                                                                                                                                          >  1:359362/1‑250 (MQ=255)
                                                                                                          tCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACaaaa                                                                                                                                         >  1:80951/1‑250 (MQ=255)
                                                                                                           cATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACaaaaa                                                                                                                                        >  1:239775/1‑250 (MQ=255)
                                                                                                            aTCGGTGAACTGCGGCAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGAGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGTTTATGAATTACCGCATTATAACCTAATAACGAGCATCTTTCATGACGGCAAACACTCGGGAAGTATTGACAAGCAAATTACAAATCATTAACAAAAAAt                                                                                                                                      >  2:372236/1‑251 (MQ=255)
                                                                                                               ggTGAAATGCGGAAGAACATCGTAGGGTTTATTGAACAGGCTCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGc                                                                                                                                   >  2:473532/1‑251 (MQ=255)
                                                                                                                  gAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGctct                                                                                                                                <  1:671599/251‑1 (MQ=255)
                                                                                                                   aaCTGCGGAAGATGATCGTAGTGTTCATGGAAAAGGATCACACGCGTGTGCTGTTTTTCTGGTTTTCGCCTGGTAGTAGGTTCGGAGATGAATCGCTTAACCTGGTTATTGCTAAAAGAAGTTTTTTGCATAGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTACTAACACGCATCTTTCATGAAGGCAAAAAATAGGGTAGTATTGACAAACAAATTACAAATCATTAACAAAAAAGTGCTCTa                                                                                                                               <  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAGTTGTTCCGGCACAAGGCAAGAAGATCAC  >  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: