Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8610GT38.5% 36.3 / 17.7 39P259Q (CCA→CAA) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base G (8/16);  new base T (15/0);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 2.52e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CACAGAACAAAATTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTT  >  NC_000913/2434624‑2435108
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                           
cacaGAACAAAATTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAcc                                                                                                                                                                                                                                                <  1:197245/250‑1 (MQ=255)
  caGAACAAAATTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGAGATCATActcgc                                                                                                                                                                                                                                             >  2:107923/1‑247 (MQ=255)
            ttAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCCGGGTAAtt                                                                                                                                                                                                                                  >  1:83297/1‑250 (MQ=255)
            ttAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATc                                                                                                                                                                                                                                  >  1:229187/1‑251 (MQ=255)
               aGAACAAAATGGCATAAACCCAGAAGCCGCACTCCGCGACCATCTCTGAAAGATCCACGTCAATGAGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGCAGATctgaacgggtcagggaattcac                                                                                                                                                                                                                                 >  2:348083/1‑230 (MQ=255)
                 aaCAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGTAGATCATACCGGCCAGGGGAATACAccc                                                                                                                                                                                                                             >  2:279113/1‑250 (MQ=255)
                   cAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATACCGGTCAGGGGAAACCACcgcg                                                                                                                                                                                                                           >  1:47354/1‑251 (MQ=255)
                       aTGGCATAAACCCAGACGCCGTACTCCGCGCCCAGCTAAGCCCGATGCACGTCACTGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAAATCACACCCGTCAGGGGAAaccaccccgcacc                                                                                                                                                                                                                       >  1:332083/1‑244 (MQ=255)
                                      aCGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCCCCACCAGATAAAGCCCTTccc                                                                                                                                                                                                        >  1:100376/1‑250 (MQ=255)
                                       cGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCAGTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGc                                                                                                                                                                                                       <  2:306427/251‑1 (MQ=255)
                                           gTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAAACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCtttt                                                                                                                                                                                                  >  1:224103/1‑251 (MQ=255)
                                                        aGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGGCAGGGGAATCCACCGCGACCAGATAAAGCCCCTCCCCTTTTTCCGttttttccc                                                                                                                                                                                     >  1:174267/1‑244 (MQ=255)
                                                                            aaTGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCCCCTTTTCCCGTTCTGCCCCCCGCGGTCTGTCCttttcc                                                                                                                                                                 >  2:293089/1‑248 (MQ=255)
                                                                                   aGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCa                                                                                                                                                          <  2:339566/251‑1 (MQ=255)
                                                                                        aaaaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAActc                                                                                                                                                      <  1:152356/250‑1 (MQ=255)
                                                                                          aaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGc                                                                                                                                                   <  2:68009/251‑1 (MQ=255)
                                                                                             cGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTTAAGATCCTACCGGt                                                                                                                                                                                                                                             >  1:235138/1‑160 (MQ=255)
                                                                                                          aatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGCGATCCTACCCGTCCGGGTAATCCACCCCGACCCGATACCGCCCTTCCCCTTTTTCCGTTGCTTCCCCCCGCGTTCTTTCCTTTGCCGCCCGCCACTCTGCTTTCCCCCTCTcccgc                                                                                                                                   >  1:109250/1‑248 (MQ=255)
                                                                                                                     tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCCTACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCGCCCACGGTCTGTCCTTTGCCGCCAGCAACTCTTCTATCCCGCTCTcccggttggtgtgggg                                                                                                                        >  1:269308/1‑237 (MQ=255)
                                                                                                                              aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGAACGGGGAAATCCACCCGACCAGATACAGCCCCTCCCCTCTTTCCGTTGCTGCCCCCAGCGGTCTGTCCTTTGCCGCCAGCACATCTGCTATCCTGCCCTccctctgggttggggccccgacatc                                                                                                               >  2:63777/1‑228 (MQ=255)
                                                                                                                                  cTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAAACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCCCTTTTTCCGTTGCTGCCGCCCGCGCTCTGTCCTTTGCCGCCCGCCACTCTGCCATTCCGCCCCCCCGCTGGGTGTGGGCACCGACtctcacc                                                                                                           >  2:204162/1‑245 (MQ=255)
                                                                                                                                           gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCCGGGGGTTCCACCGCGCCCAAATACAGCCCTTTCGCTTTTTCCGTTGCTGCCGCCCGCGTTCTTTCCTTTTCCGCCCGCAACTCTTCTATCCAGCTCTCCCGCCGGTTGTGGGGACCCACTTCCCTCCGGCTTCCg                                                                                                  >  1:284183/1‑251 (MQ=255)
                                                                                                                                           gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCCCCGCGACCCGAAACAGCCCCTCCCCTTTTTCCGTTTCTGCCGCCCGCGTTCTGTCCTTTGCCGCCAGCCACTCCCCTTTCCCGCCCTCCCGCTTGGTGTGGGGACCCACTTCCCTCATGccgccg                                                                                                  >  2:281983/1‑251 (MQ=255)
                                                                                                                                                    cAGGCTGTAACTACGTTGGCTACATCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc                                                                                         <  1:99084/251‑1 (MQ=255)
                                                                                                                                                    cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACGGGCAGGGGAATCCACCGCGACCAGATAAAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCCGCTGGTTGTGGGCACCCACTTCCCTCAGGCTGCCGACAATATTc                                                                                         >  1:334651/1‑251 (MQ=255)
                                                                                                                                                           tAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCCTAACCGCCAGGGGAATCCCCCGCGACCAGCTACAGCCCTTCCCCTTTTGCCGTTGCTGCCCCCAGCGTTCTTTCCTTTTTCGCCCGCAACTCTGCTATCCCGCTCTCCGGCTGGTTGTGGGGCCCCACTTCCATCAGGCCGCCGACCATATTCCTGAccc                                                                                  >  1:170416/1‑250 (MQ=255)
                                                                                                                                                                aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTACCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATatg                                                                              <  2:256285/250‑1 (MQ=255)
                                                                                                                                                                                 aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAATCCACCGCGGCCCGATACAGCCCCTCCGCCTTTTCCGTTGCTGCCCCCAGCGGTTTCTTCTTTTCCGCCCGCCACTCTGCTATCCCGCTCTCCGGGCGGTTGTGGGCGCACACTTTCCTCCGGCTGCCGCCCATAATCCTGACCACATGATGTTACaaagaatttgc                                                            >  1:346914/1‑243 (MQ=255)
                                                                                                                                                                                     gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCttt                                                        <  2:343395/250‑1 (MQ=255)
                                                                                                                                                                                        ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAata                                                    <  1:313743/251‑1 (MQ=255)
                                                                                                                                                                                          gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATc                                                  <  2:164544/251‑1 (MQ=255)
                                                                                                                                                                                                         tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCCTACCGGGCAGGGGAAACCACCGCGACCAGATACAACCCTTCCGCTTTTTCCGTTGCTGCCCCCCGCGTTTTGTCCTTTTCCCCCCGCAACTCTTCTATCCAGCTCTCCGGCTGGGTGTGGGGACCGACTTCCCTCAGGCTGCCGACCATATTTCTGACCATATGATGTACAAAAACATTTGCTTTAATATCTACCACCACATAAgg                                   >  1:192509/1‑251 (MQ=255)
                                                                                                                                                                                                                  aGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCCGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTGCCGCCCGCAACTCTTCTATCCAGCTCTCCGGCTGGTTGTGGGCCCCGACTTCCATCAGGCTGCCGACAAAATTCCTGACCATATGATGTACAAAAGCATTCGCTTTTATATTTACCACCACCTAAAGACCGTGAcg                          >  1:82446/1‑251 (MQ=255)
                                                                                                                                                                                                                  aGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGAcg                          <  1:40978/251‑1 (MQ=255)
                                                                                                                                                                                                                            ggCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCCCCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTACTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTa                <  2:47354/251‑1 (MQ=255)
                                                                                                                                                                                                                              cccATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAat              <  2:276557/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGt   <  1:281983/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCTGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAATGACCGTGACGCGTGACATTAATATGCATAACGtt  <  1:331862/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGtt  <  1:279113/250‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                           
CACAGAACAAAATTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTT  >  NC_000913/2434624‑2435108

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