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.

TGCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCC  >  NC_000913/2744842‑2745335
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
tGCCGAGGCTGCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGAGCAAGCTTGCGTCCCTGCGGTGACA‑GATAAATCAcc                                                                                                                                                                                                                                                       <  2:401715/250‑1 (MQ=255)
    gAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAACTGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCACCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGGCGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GCTAAATCACCTTTgg                                                                                                                                                                                                                                                  >  1:182139/1‑250 (MQ=255)
      ggCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTg‑cgc                                                                                                                                                                                                                                               >  1:240011/1‑251 (MQ=255)
       gCTTCATGTCCCAGAACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATCGACCCTTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGCATCACTTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTAGATCAAGCTTTCGTCCCTGTGGTG‑CACGATAACTCACCTTTG‑CGcc                                                                                                                                                                                                                                              >  2:291540/1‑251 (MQ=255)
         ttCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGATTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCtt                                                                                                                                                                                                                                            >  2:265311/1‑251 (MQ=255)
           cATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAAGCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCTCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGGTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCa                                                                                                                                                                                                                                          <  1:448141/251‑1 (MQ=255)
                   aGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACACACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑AACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑CTTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAACGCACCTTTG‑CGCCTTCACCcgccgc                                                                                                                                                                                                                                  >  2:293165/1‑251 (MQ=255)
                        tCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCtt                                                                                                                                                                                                                             >  2:132774/1‑251 (MQ=255)
                             gAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTtgct                                                                                                                                                                                                                        >  1:477439/1‑251 (MQ=255)
                                 aaaCCGGAAAACGGAGTCAATCAGCGTCATTGCTGGAAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCTTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCCCCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCAt                                                                                                                                                                                                                    <  1:300043/251‑1 (MQ=255)
                                   ccctggaaacggattaacccAGCGTCATTGCTGGTCCCTCACCCCCACTGAGACCGTAACCGCCGATTCGCCCCTCTTCGTCATTTTCGTTTTGGAT‑CACGCCCTCACCAATCCCTTCGTCGCTACCCCCCCCCAGAGTCAATTTTT‑G‑ATTCGTTGCCAGTCCGCCTTCGCCCTCTTGCCCAAGTTGGCGTCCTTGTGGTAACA‑GATAAATCACCTTTT‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGa                                                                                                                                                                                                                  <  1:470519/232‑1 (MQ=255)
                                           aCGGAGTCAATCAGCGTCATTGCTGGTACCTCCCCACCCCTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAC‑CCCGCGCTCATCTATACCTTCGTAGCGACCGCACCCCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACCCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCCCCTTTG‑CGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTc                                                                                                                                                                                                          <  1:296569/251‑1 (MQ=255)
                                           aCGGAGACAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAAAGGAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCACGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGCATGGCGTc                                                                                                                                                                                                          >  1:190124/1‑251 (MQ=255)
                                                     ccagcgtcttggctgtacaccccccaccacTNACAACGTAATCGCCGTTTGACCATTCTGCGTCAATTTGGGTTTGGAT‑CACGCGCTCTTTTATACCTTCGTAGCGCCCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTCCGCTGTCGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAACCCCCTTTG‑CGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGt                                                                                                                                                                                                <  2:247460/228‑1 (MQ=255)
                                                      cggcgacATTGCTGGTACCACCCCCCCCCTGNGAACGTAACCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CCGGCGCTCATCTATCCCTTCGTAGCGACCGCCCACCAGAATCAATTTTG‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCCCCTTTG‑CGCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGtt                                                                                                                                                                                               <  2:232702/245‑1 (MQ=255)
                                                      cAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGtt                                                                                                                                                                                               >  1:290510/1‑251 (MQ=255)
                                                       aGCGTCATTGCTGGTAACGCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAATACACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGTGTCCCGCAAgcgtgg                                                                                                                                                                                              >  2:372651/1‑246 (MQ=255)
                                                          gTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGCTGCAc                                                                                                                                                                                           >  2:261563/1‑251 (MQ=255)
                                                            cATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAcc                                                                                                                                                                                          >  1:366200/1‑250 (MQ=255)
                                                             aTTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTNCATGAATGGCGTCCCGCACGGCTTGCACcat                                                                                                                                                                                        >  2:231699/1‑251 (MQ=255)
                                                                     tAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAAc                                                                                                                                                                                 >  1:120551/1‑250 (MQ=255)
                                                                                     agaactaaaTCGCCGATGGACCTTTCGTCGGCAATTGCGGTT‑GGATTGACGCCCTCTTCTATCCCTTCGTAGCGACCGCACACCAGAATCATTTTTT‑G‑ATTCGTTGCCATTGCGCTGACGCCCGCTTGATCAAGCTGGCGTCCCTGTGGTGACA‑GAAAAATCCCCGTTG‑CGCCTCCACCCGCCGCGGCTTTTGCTTCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgccgcc                                                                                                                                                                <  2:203773/244‑1 (MQ=255)
                                                                                            aacgcccGATTGCCCATTCTCCGTCATTTTCGGTGTGGAT‑CCCGCGCTCATCTATCCCTTCTTAGCGACCGCCCCCCAGAATCATTTTTT‑G‑ATGCGTTGCCAGTCCGCTGCCGCCCGCTTGACCAAGCTTGCGCCTCTGGGGTGACA‑GATAAATCCCCTTTG‑CGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGa                                                                                                                                                         <  2:53885/246‑1 (MQ=255)
                                                                                                    aTTGACCATTCTTCGTCCATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTc                                                                                                                                                 >  1:396539/1‑251 (MQ=255)
                                                                                                       gccctttttccGTCATTTTCGTTTTGGAC‑CACGCGCTCATCTTTCCCTTCGTAGC‑TCCCCCCCCCAAAATCAATTTTTTG‑TTTCGTTGCCGGTTCGCTGCCGCCCGCTTAACCAAGCTTTCTTCCCTGTGTTGATA‑GAAAAACCCCCTTTG‑CGCTTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGTTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCac                                                                                                                                              <  1:410366/241‑1 (MQ=255)
                                                                                                          cATTCTTCGTCAATTTCGGTTTGGAT‑CCCGCGCTCATCTATCCCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCTTTGCCAGTTCGCTGACGCCCGCTTGATCAAGTTTGCGTCCCTGTGGTGACA‑GATAAATCCCCTTTG‑CCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggt                                                                                                                                           <  1:372402/251‑1 (MQ=255)
                                                                                                           aTTCTTTGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggt                                                                                                                                           >  1:260981/1‑250 (MQ=255)
                                                                                                           aTTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggt                                                                                                                                           >  1:199798/1‑250 (MQ=255)
                                                                                                                tcGTCAATTTCGGTTTGGAT‑CCCGCGCTCATCTATACCTTCGTAGCGCCCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCATTTCGCTGCCGCCCGCTTGATCAAGCTTGCGCCCCTGTGGTGACA‑GATAAATCCCCTTTG‑CGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggtacg                                                                                                                                        <  1:256204/248‑1 (MQ=255)
                                                                                                                  gTCAATTTCGGTTTGCAT‑AACGCGCTAATCAATACCTTCGTAGCGACCGCAAACAAAAATCAATTTTT‑G‑ATTCGTTGCCAGATAGCTGACGCCCGCTTGATCAAGCTTGAGTCCCTGTGGGGACA‑GATAAATCACCTCTG‑AGCCTTCACCCGCAGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGCTCTGCACCATCATTCACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGAACGGTGcc                                                                                                                                   >  2:318362/1‑251 (MQ=255)
                                                                                                                    cAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCgg                                                                                                                                 >  1:136632/1‑251 (MQ=255)
                                                                                                                     aaTTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCttt                                                                                                                                                                                                                                                    <  2:292990/136‑1 (MQ=255)
                                                                                                                     aaTTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCttt                                                                                                                                                                                                                                                    >  1:292990/1‑136 (MQ=255)
                                                                                                                        gccgttttGGAT‑CACGCGTCCATCTATCCCTTCTTACCGCCCGCCCCCCAGAACCATTTTTT‑G‑TTCCTTTCCCAGTCCGATACCGCCCGCTTGATCAGGCTTGCGCTCTTGTGTTGACA‑‑AAAAAACCCCTTTGGCCCCTCCCCCCCCCGCTGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGTTCGTCCACGGTACGGTGCCGGTCAt                                                                                                                             <  1:397920/246‑1 (MQ=255)
                                                                                                                                       cgcgCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTcgcg                                                                                                                >  1:398147/1‑250 (MQ=255)
                                                                                                                                             cATCTATACCTTCTTAGCGCCCGCACCCCAGAATAAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCTTTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCCCCTTGG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTc                                                                                                         <  1:420128/251‑1 (MQ=255)
                                                                                                                                             cATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGACAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGCTCATGCGTGAAGTCGCGAGGACTc                                                                                                         >  1:66359/1‑251 (MQ=255)
                                                                                                                                              aTCTATACCTTCTTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGATTCACTGTGGTGTCA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCACGGCTTTTGCGGCATGAATCGCGTCCCGCAAGGGATGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATGGCCCACGGTAAGGTTCCGGTCATGCGTGAAGGCGCAAGGGCTcc                                                                                                        >  2:43834/1‑251 (MQ=255)
                                                                                                                                                   taCCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCTGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGAAAGGACGATCGTCCACGGCACGGTTCCGGTCATGCGGGAAGtggagaggcctcacgccc                                                                                                   >  2:255792/1‑234 (MQ=255)
                                                                                                                                                    aCCTTCGTAGAGAGAGCCCACCAGAATCCATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCACTGTGGTGACA‑GATAAATCAACTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGATGCACCATCATTAACATCCCCGGTCCGCCGCCGTACGGACGATCGTCCCCGGTACGGTGCCGNTCATGCGTGAAGTCGCGAGGACTCCAGctct                                                                                                  >  2:224529/1‑251 (MQ=255)
                                                                                                                                                        tcttagcgccccccccccAGAA‑CAATTTTT‑GGATTCTTTCCCGGTTCTCTGACCCCCGCTTGACCAACCTGGCGTCCCTGTGGTGACA‑GAAAAACCCCCTTTG‑CCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGa                                                                                               <  1:369375/234‑1 (MQ=255)
                                                                                                                                                         cGTAGCGCCCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGCCCCTGTGGTGACA‑GATAAATCCCCTTTG‑CGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGAt                                                                                              <  1:275307/250‑1 (MQ=255)
                                                                                                                                                          gtagcgcccgccacccAGAATCATTTTTG‑G‑GTTCTTTGCCAGTTCGCTGCCCCCCGCTTGATAAAGTTGGCGTCCTGGGGGTGACA‑GATAAACCCCCTTGG‑CGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGc                                                                                            <  2:290510/237‑1 (MQ=255)
                                                                                                                                                          gTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGc                                                                                            >  1:131241/1‑251 (MQ=255)
                                                                                                                                                              cGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGAGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTc                                                                                          >  1:475947/1‑249 (MQ=255)
                                                                                                                                                               gACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGAACGGTGc                                                                                                                                    >  1:396704/1‑206 (MQ=255)
                                                                                                                                                                 cccccccccAGAATCAATTTTT‑G‑ATTCGTTGCCTGTCCGCTGACGCCCGCTTGATCAATCTTGCGCCCCTGTGGTGACA‑AAAAACCCCCCTTTG‑CCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTcagcag                                                                                     <  1:100754/244‑1 (MQ=255)
                                                                                                                                                                    cacaCCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGcc                                                                                  >  1:351050/1‑251 (MQ=255)
                                                                                                                                                                        ccAGAAACAATTTTT‑G‑ATTCGTTGCCATTTCGCTGACGCCCGCTTGATCACGCTTGCGTCCCTGTGGTGACA‑GAAAAACCCCCTTTG‑CGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCAttt                                                                              <  1:334505/251‑1 (MQ=255)
                                                                                                                                                                        cagaaacaaaTTTTT‑G‑ATTCGTTCCGTGTTCGCTGACGCCCGCTGTATCAACCTTGCTCCCCTGTGGTGACA‑AATAAACCCCCTTTG‑CGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACACCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCAttt                                                                              <  1:132774/243‑1 (MQ=255)
                                                                                                                                                                                  cgccgctctccctatccaggTTCGCGGCCGCCCGCTTGCTCAAGCTTGCGTCCGTGTGGTGACA‑GATAAATCCCCTTTG‑CCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGc                                                                      <  2:166303/232‑1 (MQ=255)
                                                                                                                                                                                      t‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGcc                                                                >  1:87914/1‑251 (MQ=255)
                                                                                                                                                                                           ttCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGt                                                             >  2:68839/1‑251 (MQ=255)
                                                                                                                                                                                                   cAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGc                                                                      >  1:166303/1‑234 (MQ=255)
                                                                                                                                                                                                     gTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGCGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCAGGCCAGTTCCCCCGTa                                                    >  2:303440/1‑250 (MQ=255)
                                                                                                                                                                                                           tGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGt                                              >  1:3757/1‑250 (MQ=255)
                                                                                                                                                                                                             aCGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCCCGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTact                                           >  2:342933/1‑249 (MQ=255)
                                                                                                                                                                                                                 ccGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTgcg                                       >  1:184429/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTgcg                                       >  1:143058/1‑251 (MQ=255)
                                                                                                                                                                                                                       gaccaatcTTGCGCCCCTGTGGTGACA‑AAAAAACCCCCTTTG‑CGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGaa                                  <  1:336529/243‑1 (MQ=255)
                                                                                                                                                                                                                            aGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATAGGTAATTGCGCGGAACATTTc                            >  1:240776/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cccTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGCGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAAc                     >  2:285586/1‑249 (MQ=255)
                                                                                                                                                                                                                                     cccTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGCTGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAAc                     >  1:203695/1‑249 (MQ=255)
                                                                                                                                                                                                                                     cccTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAAc                     >  1:78896/1‑249 (MQ=255)
                                                                                                                                                                                                                                              aaCA‑GATAAATCCCCTGTG‑CCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTAt          <  1:16070/250‑1 (MQ=255)
                                                                                                                                                                                                                                              aaCA‑GATAAACCCCCTTGG‑CCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTAt          <  1:354217/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaaaTCCCCTTTG‑CCCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAAt    <  1:427576/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaTCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTTAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTCATTATGCCAATcc  >  2:114711/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          cACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGCAATCGGTAATTGCGCGGAACATTTCAGGAACCAGGCTAATTAt          >  2:61591/1‑240 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
TGCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGAT‑CACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTT‑G‑ATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACA‑GATAAATCACCTTTG‑CGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCC  >  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: ‑