Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,745,0890AC31.9% 123.4 / 12.0 69V84G (GTG→GGG) trmDtRNA m(1)G37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (43/4);  new base C (0/22);  total (43/26)
Fisher's exact test for biased strand distribution p-value = 2.54e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCC  >  NC_000913/2744842‑2745335
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
tGCCGAGGCTGCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGAGCAAGCTTGCGTCCCTGCGGTGACAGATAAATCAcc                                                                                                                                                                                                                                                       <  2:401715/250‑1 (MQ=255)
    gAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAACTGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCACCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGGCGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGCTAAATCACCTTTgg                                                                                                                                                                                                                                                  >  1:182139/1‑250 (MQ=255)
      ggCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTgcgc                                                                                                                                                                                                                                               >  1:240011/1‑251 (MQ=255)
       gCTTCATGTCCCAGAACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATCGACCCTTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGCATCACTTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTAGATCAAGCTTTCGTCCCTGTGGTGCACGATAACTCACCTTTGCGcc                                                                                                                                                                                                                                              >  2:291540/1‑251 (MQ=255)
         ttCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGATTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCtt                                                                                                                                                                                                                                            >  2:265311/1‑251 (MQ=255)
           cATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAAGCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCTCACACCAGAATCAATTTTTGATTCGTTGCCAGGTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCa                                                                                                                                                                                                                                          <  1:448141/251‑1 (MQ=255)
                   aGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACACACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATAACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGCTTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAACGCACCTTTGCGCCTTCACCcgccgc                                                                                                                                                                                                                                  >  2:293165/1‑251 (MQ=255)
                        tCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCtt                                                                                                                                                                                                                             >  2:132774/1‑251 (MQ=255)
                             gAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTtgct                                                                                                                                                                                                                        >  1:477439/1‑251 (MQ=255)
                                 aaaCCGGAAAACGGAGTCAATCAGCGTCATTGCTGGAAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCTTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCAt                                                                                                                                                                                                                    <  1:300043/251‑1 (MQ=255)
                                   ccctggaaacggattaacccAGCGTCATTGCTGGTCCCTCACCCCCACTGAGACCGTAACCGCCGATTCGCCCCTCTTCGTCATTTTCGTTTTGGATCACGCCCTCACCAATCCCTTCGTCGCTACCCCCCCCCAGAGTCAATTTTTGATTCGTTGCCAGTCCGCCTTCGCCCTCTTGCCCAAGTTGGCGTCCTTGTGGTAACAGATAAATCACCTTTTCGCCTTCACCCGCCGCGGCTTTTGCTGCATGa                                                                                                                                                                                                                  <  1:470519/232‑1 (MQ=255)
                                           aCGGAGTCAATCAGCGTCATTGCTGGTACCTCCCCACCCCTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGACCCCGCGCTCATCTATACCTTCGTAGCGACCGCACCCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACCCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTc                                                                                                                                                                                                          <  1:296569/251‑1 (MQ=255)
                                           aCGGAGACAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAAAGGAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCACGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGCATGGCGTc                                                                                                                                                                                                          >  1:190124/1‑251 (MQ=255)
                                                     ccagcgtcttggctgtacaccccccaccacTNACAACGTAATCGCCGTTTGACCATTCTGCGTCAATTTGGGTTTGGATCACGCGCTCTTTTATACCTTCGTAGCGCCCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTCCGCTGTCGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAACCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGt                                                                                                                                                                                                <  2:247460/228‑1 (MQ=255)
                                                      cggcgacATTGCTGGTACCACCCCCCCCCTGNGAACGTAACCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCCGGCGCTCATCTATCCCTTCGTAGCGACCGCCCACCAGAATCAATTTTGGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGtt                                                                                                                                                                                               <  2:232702/245‑1 (MQ=255)
                                                      cAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGtt                                                                                                                                                                                               >  1:290510/1‑251 (MQ=255)
                                                       aGCGTCATTGCTGGTAACGCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAATACACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGTGTCCCGCAAgcgtgg                                                                                                                                                                                              >  2:372651/1‑246 (MQ=255)
                                                          gTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGCTGCAc                                                                                                                                                                                           >  2:261563/1‑251 (MQ=255)
                                                            cATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAcc                                                                                                                                                                                          >  1:366200/1‑250 (MQ=255)
                                                             aTTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTNCATGAATGGCGTCCCGCACGGCTTGCACcat                                                                                                                                                                                        >  2:231699/1‑251 (MQ=255)
                                                                     tAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAAc                                                                                                                                                                                 >  1:120551/1‑250 (MQ=255)
                                                                                     agaactaaaTCGCCGATGGACCTTTCGTCGGCAATTGCGGTTGGATTGACGCCCTCTTCTATCCCTTCGTAGCGACCGCACACCAGAATCATTTTTTGATTCGTTGCCATTGCGCTGACGCCCGCTTGATCAAGCTGGCGTCCCTGTGGTGACAGAAAAATCCCCGTTGCGCCTCCACCCGCCGCGGCTTTTGCTTCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgccgcc                                                                                                                                                                <  2:203773/244‑1 (MQ=255)
                                                                                            aacgcccGATTGCCCATTCTCCGTCATTTTCGGTGTGGATCCCGCGCTCATCTATCCCTTCTTAGCGACCGCCCCCCAGAATCATTTTTTGATGCGTTGCCAGTCCGCTGCCGCCCGCTTGACCAAGCTTGCGCCTCTGGGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGa                                                                                                                                                         <  2:53885/246‑1 (MQ=255)
                                                                                                    aTTGACCATTCTTCGTCCATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTc                                                                                                                                                 >  1:396539/1‑251 (MQ=255)
                                                                                                       gccctttttccGTCATTTTCGTTTTGGACCACGCGCTCATCTTTCCCTTCGTAGCTCCCCCCCCCAAAATCAATTTTTTGTTTCGTTGCCGGTTCGCTGCCGCCCGCTTAACCAAGCTTTCTTCCCTGTGTTGATAGAAAAACCCCCTTTGCGCTTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGTTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCac                                                                                                                                              <  1:410366/241‑1 (MQ=255)
                                                                                                          cATTCTTCGTCAATTTCGGTTTGGATCCCGCGCTCATCTATCCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCTTTGCCAGTTCGCTGACGCCCGCTTGATCAAGTTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggt                                                                                                                                           <  1:372402/251‑1 (MQ=255)
                                                                                                           aTTCTTTGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggt                                                                                                                                           >  1:260981/1‑250 (MQ=255)
                                                                                                           aTTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggt                                                                                                                                           >  1:199798/1‑250 (MQ=255)
                                                                                                                tcGTCAATTTCGGTTTGGATCCCGCGCTCATCTATACCTTCGTAGCGCCCGCACACCAGAATCAATTTTTGATTCGTTGCCATTTCGCTGCCGCCCGCTTGATCAAGCTTGCGCCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggtacg                                                                                                                                        <  1:256204/248‑1 (MQ=255)
                                                                                                                  gTCAATTTCGGTTTGCATAACGCGCTAATCAATACCTTCGTAGCGACCGCAAACAAAAATCAATTTTTGATTCGTTGCCAGATAGCTGACGCCCGCTTGATCAAGCTTGAGTCCCTGTGGGGACAGATAAATCACCTCTGAGCCTTCACCCGCAGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGCTCTGCACCATCATTCACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGAACGGTGcc                                                                                                                                   >  2:318362/1‑251 (MQ=255)
                                                                                                                    cAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCgg                                                                                                                                 >  1:136632/1‑251 (MQ=255)
                                                                                                                     aaTTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCttt                                                                                                                                                                                                                                                    <  2:292990/136‑1 (MQ=255)
                                                                                                                     aaTTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCttt                                                                                                                                                                                                                                                    >  1:292990/1‑136 (MQ=255)
                                                                                                                        gccgttttGGATCACGCGTCCATCTATCCCTTCTTACCGCCCGCCCCCCAGAACCATTTTTTGTTCCTTTCCCAGTCCGATACCGCCCGCTTGATCAGGCTTGCGCTCTTGTGTTGACAAAAAAACCCCTTTGGCCCCTCCCCCCCCCGCTGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGTTCGTCCACGGTACGGTGCCGGTCAt                                                                                                                             <  1:397920/246‑1 (MQ=255)
                                                                                                                                       cgcgCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTcgcg                                                                                                                >  1:398147/1‑250 (MQ=255)
                                                                                                                                             cATCTATACCTTCTTAGCGCCCGCACCCCAGAATAAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCTTTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTGGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTc                                                                                                         <  1:420128/251‑1 (MQ=255)
                                                                                                                                             cATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGACAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGCTCATGCGTGAAGTCGCGAGGACTc                                                                                                         >  1:66359/1‑251 (MQ=255)
                                                                                                                                              aTCTATACCTTCTTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGATTCACTGTGGTGTCAGATAAATCACCTTTGCGCCTTCACCCGCCACGGCTTTTGCGGCATGAATCGCGTCCCGCAAGGGATGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATGGCCCACGGTAAGGTTCCGGTCATGCGTGAAGGCGCAAGGGCTcc                                                                                                        >  2:43834/1‑251 (MQ=255)
                                                                                                                                                   taCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCTGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGAAAGGACGATCGTCCACGGCACGGTTCCGGTCATGCGGGAAGtggagaggcctcacgccc                                                                                                   >  2:255792/1‑234 (MQ=255)
                                                                                                                                                    aCCTTCGTAGAGAGAGCCCACCAGAATCCATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCACTGTGGTGACAGATAAATCAACTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGATGCACCATCATTAACATCCCCGGTCCGCCGCCGTACGGACGATCGTCCCCGGTACGGTGCCGNTCATGCGTGAAGTCGCGAGGACTCCAGctct                                                                                                  >  2:224529/1‑251 (MQ=255)
                                                                                                                                                        tcttagcgccccccccccAGAACAATTTTTGGATTCTTTCCCGGTTCTCTGACCCCCGCTTGACCAACCTGGCGTCCCTGTGGTGACAGAAAAACCCCCTTTGCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGa                                                                                               <  1:369375/234‑1 (MQ=255)
                                                                                                                                                         cGTAGCGCCCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGCCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGAt                                                                                              <  1:275307/250‑1 (MQ=255)
                                                                                                                                                          gtagcgcccgccacccAGAATCATTTTTGGGTTCTTTGCCAGTTCGCTGCCCCCCGCTTGATAAAGTTGGCGTCCTGGGGGTGACAGATAAACCCCCTTGGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGc                                                                                            <  2:290510/237‑1 (MQ=255)
                                                                                                                                                          gTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGc                                                                                            >  1:131241/1‑251 (MQ=255)
                                                                                                                                                              cGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGAGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTc                                                                                          >  1:475947/1‑249 (MQ=255)
                                                                                                                                                               gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGAACGGTGc                                                                                                                                    >  1:396704/1‑206 (MQ=255)
                                                                                                                                                                 cccccccccAGAATCAATTTTTGATTCGTTGCCTGTCCGCTGACGCCCGCTTGATCAATCTTGCGCCCCTGTGGTGACAAAAAACCCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTcagcag                                                                                     <  1:100754/244‑1 (MQ=255)
                                                                                                                                                                    cacaCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGcc                                                                                  >  1:351050/1‑251 (MQ=255)
                                                                                                                                                                        ccAGAAACAATTTTTGATTCGTTGCCATTTCGCTGACGCCCGCTTGATCACGCTTGCGTCCCTGTGGTGACAGAAAAACCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCAttt                                                                              <  1:334505/251‑1 (MQ=255)
                                                                                                                                                                        cagaaacaaaTTTTTGATTCGTTCCGTGTTCGCTGACGCCCGCTGTATCAACCTTGCTCCCCTGTGGTGACAAATAAACCCCCTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACACCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCAttt                                                                              <  1:132774/243‑1 (MQ=255)
                                                                                                                                                                                  cgccgctctccctatccaggTTCGCGGCCGCCCGCTTGCTCAAGCTTGCGTCCGTGTGGTGACAGATAAATCCCCTTTGCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGc                                                                      <  2:166303/232‑1 (MQ=255)
                                                                                                                                                                                      tGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGcc                                                                >  1:87914/1‑251 (MQ=255)
                                                                                                                                                                                           ttCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGt                                                             >  2:68839/1‑251 (MQ=255)
                                                                                                                                                                                                   cAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGc                                                                      >  1:166303/1‑234 (MQ=255)
                                                                                                                                                                                                     gTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGCGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCAGGCCAGTTCCCCCGTa                                                    >  2:303440/1‑250 (MQ=255)
                                                                                                                                                                                                           tGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGt                                              >  1:3757/1‑250 (MQ=255)
                                                                                                                                                                                                             aCGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCCCGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTact                                           >  2:342933/1‑249 (MQ=255)
                                                                                                                                                                                                                 ccGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTgcg                                       >  1:184429/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTgcg                                       >  1:143058/1‑251 (MQ=255)
                                                                                                                                                                                                                       gaccaatcTTGCGCCCCTGTGGTGACAAAAAAACCCCCTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGaa                                  <  1:336529/243‑1 (MQ=255)
                                                                                                                                                                                                                            aGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATAGGTAATTGCGCGGAACATTTc                            >  1:240776/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cccTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGCGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAAc                     >  2:285586/1‑249 (MQ=255)
                                                                                                                                                                                                                                     cccTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGCTGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAAc                     >  1:203695/1‑249 (MQ=255)
                                                                                                                                                                                                                                     cccTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAAc                     >  1:78896/1‑249 (MQ=255)
                                                                                                                                                                                                                                              aaCAGATAAATCCCCTGTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTAt          <  1:16070/250‑1 (MQ=255)
                                                                                                                                                                                                                                              aaCAGATAAACCCCCTTGGCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTAt          <  1:354217/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaaaTCCCCTTTGCCCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAAt    <  1:427576/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaTCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTTAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTCATTATGCCAATcc  >  2:114711/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          cACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGCAATCGGTAATTGCGCGGAACATTTCAGGAACCAGGCTAATTAt          >  2:61591/1‑240 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
TGCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCC  >  NC_000913/2744842‑2745335

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