Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,071 A→T 10.0% noncoding (85/1542 nt) rrsG ← 16S ribosomal RNA of rrnG operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,731,0710AT10.0% 144.1 / 5.8 60noncoding (85/1542 nt)rrsG16S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base T (3/3);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGA  >  NC_000913/2730924‑2731211
                                                                                                                                                   |                                                                                                                                            
cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCaa                                                                                                                                             <  2:448306/149‑1 (MQ=40)
   ttgcGATGGCAGGGGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTCATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgct                                                                                                                                          <  1:312963/146‑1 (MQ=255)
   aTCCGATGGCAAGGGGCCCGAGGGTCCCGCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTCTTCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGTGAAGCAGCAAgct                                                                                                                                          <  2:385202/149‑1 (MQ=11)
   aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgct                                                                                                                                          <  1:641505/149‑1 (MQ=39)
   aTCCAATGGCAAGAGGCCCGAAGTTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgct                                                                                                                                          <  2:570210/149‑1 (MQ=255)
    tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTGCAAGTAGTTATCCCCCCCCATCAGGCAGCGACCGAGACATCACTCACCCGTCCGCCACTCGTCAGCGAAGGAGCAAgctc                                                                                                                                         >  2:750518/1‑148 (MQ=11)
        aTGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAGCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTc                                                                                                                                     <  2:221296/149‑1 (MQ=255)
         tGGCAAGAGGCCCGAAGGTCCCGCTCGTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcg                                                                                                                                    <  2:162009/149‑2 (MQ=255)
          cgCAAGAGGCCCGAAGGCCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCt                                                                                                                                   <  1:564452/148‑1 (MQ=255)
           gCAAGAGGCCCGAAGGTCCCCCTCTTGGGTCTGGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTg                                                                                                                                  <  2:483496/149‑1 (MQ=255)
            cAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGt                                                                                                                                 >  1:591310/1‑149 (MQ=40)
             aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGtt                                                                                                                                >  1:236533/1‑149 (MQ=40)
             aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGtt                                                                                                                                >  1:484780/1‑149 (MQ=40)
                    ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTc                                                                                                                         <  1:430780/149‑1 (MQ=255)
                             ccccTCTTTGGTCTTGCGACATTGTGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAt                                                                                                                >  2:717487/1‑149 (MQ=255)
                                 tctTTGGTTGTGCTACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATgtgt                                                                                                            <  1:488322/149‑1 (MQ=255)
                                 tctTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATgtgt                                                                                                            <  2:26898/149‑1 (MQ=40)
                                 tctTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATgtgt                                                                                                            <  2:471941/149‑1 (MQ=255)
                                 tctTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGGAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATgtgt                                                                                                            <  1:697381/149‑1 (MQ=255)
                                             cGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCCTGTGTTAGGCCTgccgc                                                                                                >  2:181718/1‑149 (MQ=255)
                                               aCATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCATCTTCCCAGACATTACTCACCCGTCCGCCCCTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCACGGTCAAGTCCTGCCGCCa                                                                                              <  2:677428/149‑1 (MQ=14)
                                               aCATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCa                                                                                              <  1:323014/149‑1 (MQ=255)
                                               aCATTATCCGGTATTAGCTCCCATTTCCAGCAGTTAGCCCCCTCCACCAGGCAGCCTCCCAGCCACTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCa                                                                                              <  1:212939/149‑1 (MQ=14)
                                                    aTGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGtt                                                                                         <  2:71002/149‑1 (MQ=39)
                                                                    cGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGAt                                                                         <  2:345059/149‑1 (MQ=38)
                                                                      tttCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGCGCCATGATCa                                                                       >  1:387476/1‑149 (MQ=255)
                                                                      tttCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCACCGTCCAATCTGCGCCATGATCa                                                                       >  1:242301/1‑149 (MQ=255)
                                                                              atttATCCCCCCCCACCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTc                                                               <  2:523102/147‑1 (MQ=255)
                                                                                  aTCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAAtt                                                           >  2:745572/1‑149 (MQ=38)
                                                                                   ccccccTCCATCCGGCCGCTTCCCAGACATTACTCCCCCGGCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGCTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAAttt                                                          <  2:740589/148‑1 (MQ=11)
                                                                                   ccccccTCCACCCGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAAttt                                                          <  2:340807/148‑1 (MQ=255)
                                                                                        cTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAg                                                     >  2:715708/1‑149 (MQ=38)
                                                                                           cctCAGGCAGCTTCCCAGACATTACCCACCCGTCCGCCACTCGTCAGCGAGACGGCGAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGttt                                                  <  1:607327/147‑1 (MQ=11)
                                                                                            aTCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATAAAACACTTCAATTAAAAAGttag                                                 >  1:258896/1‑147 (MQ=255)
                                                                                            aTCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTCCCGCCAGCGTTCCATCTGAGCCATGATCAAACTCTTACATTTAAAAGATTg                                                 >  2:145290/1‑149 (MQ=11)
                                                                                             ccAGGCCGCTTCCCAGCCATCACCCGCCCGTCCGCCACTCGTGAGCGATGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCGGCGTTCATTCTGAGCCATGATCGAACTCTTCAATTTAAAAGTTTGa                                                <  2:585415/148‑1 (MQ=11)
                                                                                                        cccAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAgca                                     <  1:745572/149‑3 (MQ=21)
                                                                                                         ccAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAAt                                    >  2:530528/1‑149 (MQ=35)
                                                                                                            gACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGTTCAAACTCTTCATTTTAAAGGTTTGATGTTCAAAGAATTaa                                 >  2:60045/1‑149 (MQ=255)
                                                                                                             aCATTACTCCCCCGTCCGCCACCAGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTaaa                                <  2:350946/149‑1 (MQ=255)
                                                                                                                    tCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTa                         >  1:704473/1‑149 (MQ=35)
                                                                                                                      aCCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGCCTTGCATGTGTTAGGCCTGACGCCAGCGTTCAATCTGAGCCATTTTCATACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTCAt                       >  1:614731/1‑149 (MQ=11)
                                                                                                                      aCCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACGCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAAt                       >  2:391912/1‑149 (MQ=21)
                                                                                                                        ccGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGGTCAAAGAATTAAACTTCGTAATGa                     >  1:763607/1‑149 (MQ=25)
                                                                                                                        ccGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGa                     <  1:518404/149‑1 (MQ=34)
                                                                                                                         cGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCCTGGTCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGaa                    >  2:90237/1‑149 (MQ=25)
                                                                                                                         cGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGaa                    >  2:697281/1‑149 (MQ=35)
                                                                                                                                cACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACgtg             >  2:465620/1‑149 (MQ=35)
                                                                                                                                    cGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCa         >  1:52624/1‑149 (MQ=35)
                                                                                                                                     taacaggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGCGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCAc        >  2:281827/7‑149 (MQ=11)
                                                                                                                                     taacaggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCAc        >  2:701146/7‑149 (MQ=255)
                                                                                                                                      tCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAGCTCTTCAATTTAAAAGTTTGATGCTCAGAGAATTGAGCTTCGTAATGAATTACGTGTTCAct       >  1:530529/1‑149 (MQ=255)
                                                                                                                                      tCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCAct       >  2:334946/1‑149 (MQ=35)
                                                                                                                                      tCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCAct       >  1:279020/1‑149 (MQ=35)
                                                                                                                                       cggCGAAGCAGCTAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTGAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAGCTTCGTAATGAATTACGTGTTCActc      <  2:178212/147‑1 (MQ=255)
                                                                                                                                         gCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCtt    <  2:259027/149‑1 (MQ=35)
                                                                                                                                         aggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCtt    >  1:607772/3‑149 (MQ=255)
                                                                                                                                          tgAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCTATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTg   <  1:749631/148‑1 (MQ=25)
                                                                                                                                          cGAAGCAGCAAGCTGCTTTCGGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTg   <  1:191983/149‑1 (MQ=25)
                                                                                                                                          cGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTg   <  1:744545/149‑1 (MQ=35)
                                                                                                                                           gAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTga  <  2:159039/149‑1 (MQ=34)
                                                                                                                                                   |                                                                                                                                            
CACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGA  >  NC_000913/2730924‑2731211

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