Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,712,133 G→C 5.6% P211R (CCG→CGG)  yfiC ← tRNA m(6)A37 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,712,1330GC 5.6% 278.3 / 2.9 83P211R (CCG→CGG) yfiCtRNA m(6)A37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (38/40);  new base C (2/3);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTC  >  NC_000913/2711984‑2712281
                                                                                                                                                     |                                                                                                                                                    
tGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGTCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGGCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCg                                                                                                                                                      <  1:1613258/150‑1 (MQ=255)
 gCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCgg                                                                                                                                                     <  1:839235/150‑1 (MQ=255)
    cACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGga                                                                                                                                                  >  1:582421/1‑150 (MQ=255)
    cACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGga                                                                                                                                                  <  2:384150/150‑1 (MQ=255)
      cccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACCAATGACTAACCGATCGCTAAAGCATTCTCCCGCCTGCGGGgaga                                                                                                                                                >  1:1422845/1‑150 (MQ=255)
      cccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCATGGGTCAGCGCCTTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTATCCCGCCTGCGGGgaga                                                                                                                                                >  2:1118248/1‑150 (MQ=255)
        cATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAAt                                                                                                                                              >  1:1125103/1‑150 (MQ=255)
         aTTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATg                                                                                                                                             <  2:568153/150‑1 (MQ=255)
          ttCCGGTTCGTCGATCCTTTCCCCCTGCAATCATTAGATAAACAGATAAAAAGCCTTGGTCAGGGCCGTTTACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGc                                                                                                                                            <  2:1154377/150‑1 (MQ=255)
           tCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGcc                                                                                                                                           >  2:692712/1‑150 (MQ=255)
            tctgTCCGTATATCCTTTCCCCCTTCAATCATTACATAAACATATAAAAAGCCTGGTTCAGCGCCGTATACGATTTGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCTCTAAAGCATTCTCCCAGCAGTCGGGAGAATGctt                                                                                                                                          <  2:30859/147‑3 (MQ=255)
             cGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcag                                                                                                                                         <  1:585982/150‑1 (MQ=255)
              ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                        >  2:144734/1‑150 (MQ=255)
                 ccGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCagcagc                                                                                                                                     <  1:736636/150‑1 (MQ=255)
                      gATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                        >  2:167663/1‑142 (MQ=255)
                      gATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                        <  1:167663/142‑1 (MQ=255)
                          cTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGc                                                                                                                            >  1:419198/1‑150 (MQ=255)
                          cTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGc                                                                                                                            >  1:232097/1‑150 (MQ=255)
                          cTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGc                                                                                                                            >  1:724481/1‑150 (MQ=255)
                          cTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGc                                                                                                                            >  2:396132/1‑150 (MQ=255)
                           cttCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCg                                                                                                                           <  2:504637/149‑1 (MQ=255)
                             tCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGc                                                                                                                         <  1:1639501/150‑1 (MQ=255)
                             tCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGc                                                                                                                         <  1:1400703/150‑1 (MQ=255)
                                cccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                                      <  2:1499938/150‑1 (MQ=255)
                                    gCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTcgcg                                                                                                                  >  2:1390917/1‑150 (MQ=255)
                                      aaTCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCt                                                                                                                <  1:999870/150‑1 (MQ=255)
                                         cATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCg                                                                                                             >  1:71124/1‑150 (MQ=255)
                                               aTAATCAGTTAATATGTCTGGGTCAGCGGCGTATACGCTTCGGAATAGTTCTTGTCTGGCCCACGTATGTCTAAGCGATCGCTAAAGGAGTCTCCCGCCTGCGGGGTGAATGCGAGGAGCTGCCGATGCGGCAGTCGCGCTGCGGTTTcc                                                                                                       <  2:923195/150‑1 (MQ=255)
                                               aTAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTcc                                                                                                       <  1:390901/150‑1 (MQ=255)
                                                tAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGaa                                                                                                                                               <  2:658399/109‑1 (MQ=255)
                                                tAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGaa                                                                                                                                               >  1:658399/1‑109 (MQ=255)
                                                 aaaCAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGc                                                                                                     >  2:1611878/1‑150 (MQ=255)
                                                  aaCAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGcc                                                                                                    >  1:238753/1‑150 (MQ=255)
                                                     aGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGACCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCACCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCcaca                                                                                                 <  2:741004/150‑1 (MQ=255)
                                                      gATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAt                                                                                                <  2:34075/150‑1 (MQ=255)
                                                       aTAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCTTTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtc                                                                                               <  2:218223/150‑1 (MQ=255)
                                                        tAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtcc                                                                                              >  1:462725/1‑150 (MQ=255)
                                                          aaaaGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtccgt                                                                                            >  1:1081912/1‑150 (MQ=255)
                                                          aaaaGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtccgt                                                                                            >  1:248514/1‑150 (MQ=255)
                                                              gCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGc                                                                                        <  2:1288082/150‑1 (MQ=255)
                                                                   ggTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACg                                                                                   <  1:273979/150‑1 (MQ=255)
                                                                    gTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGt                                                                                  >  1:1501182/1‑150 (MQ=255)
                                                                     tCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTa                                                                                 <  1:1572692/150‑1 (MQ=255)
                                                                     tCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTa                                                                                 <  2:462725/150‑1 (MQ=255)
                                                                        gcgcCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAAt                                                                              <  1:1153980/150‑1 (MQ=255)
                                                                             gTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAg                                                                         <  2:41715/150‑1 (MQ=255)
                                                                                 aCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCa                                                                     <  2:238753/150‑1 (MQ=255)
                                                                                    cTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGc                                                                  >  2:433665/1‑150 (MQ=255)
                                                                                    cTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGc                                                                  <  2:1125103/150‑1 (MQ=255)
                                                                                      tCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCtt                                                                <  1:1390917/150‑1 (MQ=255)
                                                                                       cGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTa                                                               <  2:995140/150‑1 (MQ=255)
                                                                                       cGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTa                                                               <  2:71124/150‑1 (MQ=255)
                                                                                          aaTAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATg                                                            <  1:144734/150‑1 (MQ=255)
                                                                                          aaTAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATg                                                            >  2:229351/1‑150 (MQ=255)
                                                                                            tAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCAGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGcc                                                          <  1:1095974/150‑1 (MQ=255)
                                                                                              gTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAg                                                        >  2:292300/1‑150 (MQ=255)
                                                                                                   ggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCg                                                   >  1:1009051/1‑150 (MQ=255)
                                                                                                   ggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCg                                                   <  1:5061/150‑1 (MQ=255)
                                                                                                    gtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGt                                                  <  2:248514/150‑1 (MQ=255)
                                                                                                     tctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTa                                                 >  1:1472060/1‑150 (MQ=255)
                                                                                                       tggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTaaa                                               >  2:352877/1‑150 (MQ=255)
                                                                                                           ccACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCa                                           <  1:378601/150‑1 (MQ=255)
                                                                                                              cGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTa                                        <  1:433665/150‑1 (MQ=255)
                                                                                                                  tGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCaa                                    <  2:1081912/150‑1 (MQ=255)
                                                                                                                    aCTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGTAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATc                                  >  2:310448/1‑150 (MQ=255)
                                                                                                                    aCTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATc                                  <  1:292300/150‑1 (MQ=255)
                                                                                                                     cTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCt                                 >  1:266421/1‑150 (MQ=255)
                                                                                                                          cGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGctct                            >  1:91453/1‑150 (MQ=255)
                                                                                                                             tCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCGACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgc                         <  2:665568/150‑1 (MQ=255)
                                                                                                                              cGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgca                        <  2:525623/150‑1 (MQ=255)
                                                                                                                                  aaaGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCAc                    >  2:581000/1‑150 (MQ=255)
                                                                                                                                     gCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTAc                 >  2:1025804/1‑150 (MQ=255)
                                                                                                                                      cATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACg                >  1:136202/1‑150 (MQ=255)
                                                                                                                                      cATTCTCCCCGCAGGCGGGAGAATGCTAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACg                >  1:48688/1‑150 (MQ=255)
                                                                                                                                      cATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACg                >  1:1054914/1‑150 (MQ=255)
                                                                                                                                          ctcCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGa            >  2:920311/1‑150 (MQ=255)
                                                                                                                                            cccGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGaaa          <  1:1463130/150‑1 (MQ=255)
                                                                                                                                                ccTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATc      <  1:1611878/150‑1 (MQ=255)
                                                                                                                                                ccTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATc      >  1:1459826/1‑150 (MQ=255)
                                                                                                                                                 cTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATcc     <  2:419198/150‑1 (MQ=255)
                                                                                                                                                  tGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCCCATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCCGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATccc    >  2:159291/1‑150 (MQ=255)
                                                                                                                                                  tGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATccc    <  2:724481/150‑1 (MQ=255)
                                                                                                                                                   gCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCct   >  1:964655/1‑150 (MQ=255)
                                                                                                                                                    cGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCctc  >  1:426017/1‑150 (MQ=255)
                                                                                                                                                    cGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCctc  >  1:1011964/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTC  >  NC_000913/2711984‑2712281

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