Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,712,130 T→A 6.7% Q212L (CAG→CTG)  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,1300TA 6.7% 502.2 / 17.7 162Q212L (CAG→CTG) yfiCtRNA m(6)A37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (86/65);  new base A (6/5);  total (92/70)
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

Only 100 of 162 total aligned reads displayed.

TGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCC  >  NC_000913/2711984‑2712278
                                                                                                                                                  |                                                                                                                                                    
tGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCg                                                                                                                                                   >  1:813114/1‑150 (MQ=255)
tGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCg                                                                                                                                                   >  2:1427929/1‑150 (MQ=255)
tGCTCACCCATTCCGGTCCGTCGATCCTTTCCCACGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCg                                                                                                                                                   <  2:2336798/150‑1 (MQ=255)
 gCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCgg                                                                                                                                                  <  1:249810/150‑1 (MQ=255)
  cTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCggg                                                                                                                                                 >  1:1937591/1‑150 (MQ=255)
   tCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCgggg                                                                                                                                                <  1:1893819/150‑1 (MQ=255)
    cACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGga                                                                                                                                               >  1:60234/1‑150 (MQ=255)
      cccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgaga                                                                                                                                             >  1:930479/1‑150 (MQ=255)
      cccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgaga                                                                                                                                             >  1:1580353/1‑150 (MQ=255)
      cccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGACTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGACCACGAATGACTAAGCGATCGCTAAAGCATTCTCACGCCTGCGGGgaga                                                                                                                                             >  2:285213/1‑150 (MQ=255)
        cATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAAt                                                                                                                                           >  2:47373/1‑150 (MQ=255)
           tCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGcc                                                                                                                                        <  1:1441456/150‑1 (MQ=255)
             cGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGcttt                                                                                                                                      <  2:1962270/150‑4 (MQ=255)
             cGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCGGGGTCAGCGCCGTATACGCGTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcag                                                                                                                                      <  2:1831511/150‑1 (MQ=255)
              ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                     >  2:2054392/1‑150 (MQ=255)
              ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                     >  2:2045938/1‑150 (MQ=255)
              ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                     <  2:1765823/150‑1 (MQ=255)
              ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                     >  2:1542696/1‑150 (MQ=255)
               gTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagca                                                                                                                                    <  1:2125131/150‑1 (MQ=255)
               gTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagca                                                                                                                                    >  1:472571/1‑150 (MQ=255)
               gTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagca                                                                                                                                    <  1:1107631/150‑1 (MQ=255)
               gTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATAAGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGcacgca                                                                                                                                    >  2:1726660/1‑145 (MQ=255)
                    tCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCAcc                                                                                                                               >  1:2237413/1‑150 (MQ=255)
                    tCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCAcc                                                                                                                               >  1:435414/1‑150 (MQ=255)
                     cGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCAccc                                                                                                                              >  2:1640910/1‑150 (MQ=255)
                     cGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCAccc                                                                                                                              >  2:2474276/1‑150 (MQ=255)
                                cccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                                   >  1:409699/1‑150 (MQ=255)
                                 ccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTc                                                                                                                  >  1:2428920/1‑150 (MQ=255)
                                 ccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTc                                                                                                                  <  2:1345164/150‑1 (MQ=255)
                                     cAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCgcgc                                                                                                              >  1:332094/1‑150 (MQ=255)
                                          aTTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGt                                                                                                         <  1:779979/150‑1 (MQ=255)
                                             aCATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTTGTCTGGCCCACGAATGAGTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGGCGCGCTTCGtttt                                                                                                      <  2:1148743/150‑1 (MQ=255)
                                             aCATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGtttt                                                                                                      >  2:771035/1‑150 (MQ=255)
                                             aCATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGtttt                                                                                                      <  2:182890/150‑1 (MQ=255)
                                             aCATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGtttt                                                                                                      >  1:841679/1‑150 (MQ=255)
                                               aTAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTcc                                                                                                    >  2:9160/1‑150 (MQ=255)
                                                tAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCg                                                                                                   >  1:462079/1‑150 (MQ=255)
                                                 aaaCAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGc                                                                                                  >  1:17561/1‑150 (MQ=255)
                                                 aaaCAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGc                                                                                                  >  2:1198469/1‑150 (MQ=255)
                                                  aaCAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGcc                                                                                                 <  2:472571/150‑1 (MQ=255)
                                                    cAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCcac                                                                                               <  1:2260854/150‑1 (MQ=255)
                                                    cAGATAAAAACCCTGGGTCAGCTCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACTAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTTCGGGGAGAATGCCAGCAGCACCCCATGCGGCAGTCGCGCTTCGTTTTCCGCcaa                                                                                               >  2:1047255/1‑149 (MQ=255)
                                                      gATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAt                                                                                             >  1:806462/1‑150 (MQ=255)
                                                           aaaGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtccgtc                                                                                        <  2:1117916/150‑1 (MQ=255)
                                                             aGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCGCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtccgtccg                                                                                      <  2:2142989/150‑1 (MQ=255)
                                                              gCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGc                                                                                     <  1:1639405/150‑1 (MQ=255)
                                                                cTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCaa                                                                                   >  1:2476168/1‑150 (MQ=255)
                                                                cTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCaa                                                                                   >  1:1237567/1‑150 (MQ=255)
                                                                   ggTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACg                                                                                >  2:1561891/1‑150 (MQ=255)
                                                                   ggTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACACGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACg                                                                                >  2:687336/1‑150 (MQ=255)
                                                                        gcgcCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAAt                                                                           >  1:644991/1‑150 (MQ=255)
                                                                         cgcCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATg                                                                          <  1:2054392/150‑1 (MQ=255)
                                                                         cgcCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATg                                                                          <  2:2132878/150‑1 (MQ=255)
                                                                         cgcCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATg                                                                          <  1:2045938/150‑1 (MQ=255)
                                                                           ccGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGTTCGCTAAAGCATTCTCCCGGCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTGGCGCTTCGTTTTCGGCCACATCCGTCCGCAAACGTAAATGcc                                                                        <  1:1535265/150‑1 (MQ=255)
                                                                           ccGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGcc                                                                        <  2:813114/150‑1 (MQ=255)
                                                                             gTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACACCCGTCCGCAAACGTAAATGCCAg                                                                      <  2:2483138/150‑1 (MQ=255)
                                                                               ataCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGcc                                                                    <  2:2237413/150‑1 (MQ=255)
                                                                                taCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAACTGCCAGccc                                                                   >  1:2297835/1‑150 (MQ=255)
                                                                                taCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATAGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGccc                                                                   >  2:1308776/1‑150 (MQ=255)
                                                                                 aCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCa                                                                  >  1:1359222/1‑150 (MQ=255)
                                                                                   gCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATg                                                                <  1:1814/150‑1 (MQ=255)
                                                                                   gCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATg                                                                >  1:1160864/1‑150 (MQ=255)
                                                                                   gCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATg                                                                >  2:64405/1‑150 (MQ=255)
                                                                                       cGGAATAGTTCTGGTCTGGCCCACTAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTa                                                            <  1:724933/150‑1 (MQ=255)
                                                                                       cGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTa                                                            <  2:201223/150‑1 (MQ=255)
                                                                                          aaTAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATg                                                         <  1:2082093/150‑1 (MQ=255)
                                                                                            tAGTTCTGGTCTGGCCCACGAATGCCTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCCCCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGcc                                                       >  2:2342386/1‑150 (MQ=255)
                                                                                              gTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAg                                                     >  1:14399/1‑150 (MQ=255)
                                                                                                  tggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTcc                                                 <  1:1561891/150‑1 (MQ=255)
                                                                                                    gtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGt                                               >  1:1713277/1‑150 (MQ=255)
                                                                                                        ggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTaaaa                                           >  1:1482368/1‑150 (MQ=255)
                                                                                                        ggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTaaaa                                           >  1:1473210/1‑150 (MQ=255)
                                                                                                        ggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTaaaa                                           >  1:1656277/1‑150 (MQ=255)
                                                                                                         gCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAAc                                          <  1:1829325/150‑1 (MQ=255)
                                                                                                         gCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTACGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAAc                                          >  2:1914473/1‑150 (MQ=255)
                                                                                                          cccACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAAGGCCAGCTCCTTaaccca                                         >  2:842160/1‑146 (MQ=255)
                                                                                                           ccACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCAAATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCa                                        >  1:1969511/1‑150 (MQ=255)
                                                                                                              cGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTa                                     <  1:1769370/150‑1 (MQ=255)
                                                                                                               gAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTAc                                    >  1:1805705/1‑150 (MQ=255)
                                                                                                                aaTGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTAcc                                   <  1:145729/150‑1 (MQ=255)
                                                                                                                  tGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCaa                                 >  1:308376/1‑150 (MQ=255)
                                                                                                                    aCTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATc                               >  2:1631765/1‑150 (MQ=255)
                                                                                                                     cTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCt                              >  2:1061886/1‑150 (MQ=255)
                                                                                                                          cGATCGCTAAAGCATTCTCCCGGCTTCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTAGCAATCTGctct                         <  2:1014171/150‑1 (MQ=255)
                                                                                                                              cGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCCGCTCCGTAAAACCATTACCAATCTGCTCTGgca                     >  2:366322/1‑150 (MQ=255)
                                                                                                                              cGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgca                     <  1:1075137/150‑1 (MQ=255)
                                                                                                                              cGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgca                     <  2:806462/150‑1 (MQ=255)
                                                                                                                              cGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgca                     >  2:1811264/1‑150 (MQ=255)
                                                                                                                              cGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgca                     >  2:1763066/1‑150 (MQ=255)
                                                                                                                                cTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGTCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcagc                   >  1:1369726/1‑150 (MQ=255)
                                                                                                                                cTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCGCTGgcagc                   >  1:221162/1‑150 (MQ=255)
                                                                                                                                 gctggCATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcagca                  >  1:754479/5‑150 (MQ=255)
                                                                                                                                  aaaGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCAc                 <  2:415205/150‑1 (MQ=255)
                                                                                                                                  aaaGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCAc                 >  2:924649/1‑150 (MQ=255)
                                                                                                                                    ggCATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTa               >  2:2432458/2‑150 (MQ=255)
                                                                                                                                      cATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACg             <  1:1058489/150‑1 (MQ=255)
                                                                                                                                      cATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACg             >  2:932313/1‑150 (MQ=255)
                                                                                                                                          ctcCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGa         >  1:622717/1‑150 (MQ=255)
                                                                                                                                                 cTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATcc  >  1:2436241/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCC  >  NC_000913/2711984‑2712278

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: