Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,712,128 C→G 7.8% A213P (GCG→CCG)  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,1280CG 7.8% 284.9 / 10.6 86A213P (GCG→CCG) yfiCtRNA m(6)A37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (42/37);  new base G (4/3);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

ACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATC  >  NC_000913/2711982‑2712277
                                                                                                                                                  |                                                                                                                                                     
aCTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTg                                                                                                                                                    >  1:728586/1‑150 (MQ=255)
  tGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCg                                                                                                                                                  >  1:2993/1‑150 (MQ=255)
   gCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCgg                                                                                                                                                 >  2:1909551/1‑150 (MQ=255)
    cTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCggg                                                                                                                                                >  2:185725/1‑150 (MQ=255)
     tCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCgggg                                                                                                                                               <  1:1331212/150‑1 (MQ=255)
      cACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGga                                                                                                                                              >  2:892227/1‑150 (MQ=255)
      cACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGga                                                                                                                                              <  2:326443/150‑1 (MQ=255)
      cACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGga                                                                                                                                              >  1:537103/1‑150 (MQ=255)
      cACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATACTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGga                                                                                                                                              >  2:109708/1‑150 (MQ=255)
      cACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATAATTACATAAACAGATAAAAAGCATGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGATTGAATAAGCGATCGCTAAAGCATTCTCCCGCATGCGGGga                                                                                                                                              >  2:898956/1‑150 (MQ=255)
       aCCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgag                                                                                                                                             <  1:454764/150‑1 (MQ=255)
       aCCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGGCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgag                                                                                                                                             <  2:1970750/150‑1 (MQ=255)
        cccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgaga                                                                                                                                            <  1:1697427/150‑1 (MQ=255)
          cATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACTCATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGACTGGCCAACGAATGACTAAGCGATCGCTAAAGCATTATACCGCCTGCGGGGAGAAt                                                                                                                                          >  2:1714898/1‑150 (MQ=255)
          cATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAAt                                                                                                                                          >  1:264881/1‑150 (MQ=255)
              ccGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCca                                                                                                                                      <  1:1617064/150‑1 (MQ=255)
               cGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGcttt                                                                                                                                     <  2:1644853/150‑4 (MQ=255)
                ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                    <  2:1041520/150‑1 (MQ=255)
                ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                    >  1:1188061/1‑150 (MQ=255)
                 gTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagca                                                                                                                                   >  2:931850/1‑150 (MQ=255)
                  tCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagcag                                                                                                                                  <  2:264881/150‑1 (MQ=255)
                  tCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagcag                                                                                                                                  <  1:1247935/150‑1 (MQ=255)
                      tCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGTCAGCAGCAcc                                                                                                                              <  1:1909551/150‑1 (MQ=255)
                       cGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCAccc                                                                                                                             <  1:193810/150‑1 (MQ=255)
                       cGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGGCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCAccc                                                                                                                             <  1:1186252/150‑1 (MQ=255)
                                 ccccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCCGCCCCCCCCGATGCGGcgg                                                                                                                   >  1:1192386/1‑148 (MQ=255)
                                   ccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGc                                                                                                                                        >  2:797260/1‑127 (MQ=255)
                                   ccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGc                                                                                                                                        <  1:797260/127‑1 (MQ=255)
                                    cGGCAATCATTACATAAACAGATAAAAAGCCTGGGGCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAATGCATTCTCCCGCCTGCGGGGAGAATGCCTGCAGCACCCGATGCGGCAGTcg                                                                                                                <  1:1376911/150‑1 (MQ=255)
                                            aTTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGt                                                                                                        >  1:1921747/1‑150 (MQ=255)
                                              tACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGttt                                                                                                      <  1:441709/150‑1 (MQ=255)
                                               aCATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGtttt                                                                                                     <  1:373915/150‑1 (MQ=255)
                                                cATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTc                                                                                                    >  2:923572/1‑150 (MQ=255)
                                                  tAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCg                                                                                                  >  1:39438/1‑150 (MQ=255)
                                                   aaaCAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGc                                                                                                 <  2:1492395/150‑1 (MQ=255)
                                                          tAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtcc                                                                                          >  2:1768014/1‑150 (MQ=255)
                                                           aaaaaGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtccg                                                                                         >  2:1967974/1‑150 (MQ=255)
                                                            aaaaGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtccgt                                                                                        >  1:1234657/1‑150 (MQ=255)
                                                                  cTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCaa                                                                                  >  1:1222270/1‑150 (MQ=255)
                                                                        cAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTaa                                                                            >  1:1647787/1‑150 (MQ=255)
                                                                                tataCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGc                                                                    >  2:1570900/1‑150 (MQ=255)
                                                                                tataCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGc                                                                    <  2:728586/150‑1 (MQ=255)
                                                                                 ataCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGcc                                                                   <  1:931850/150‑1 (MQ=255)
                                                                                   aCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCa                                                                 <  1:1967974/150‑1 (MQ=255)
                                                                                    cGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCAt                                                                >  1:14927/1‑150 (MQ=255)
                                                                                     gCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATg                                                               <  1:604251/150‑1 (MQ=255)
                                                                                        tCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCACGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCtt                                                            >  1:1968908/1‑150 (MQ=255)
                                                                                         cGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCACCTGCGGGGAGAATGCCAGCCGCACCCGATACGGCAGTCGCGCTTCGTTTTCAGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTa                                                           >  1:675302/1‑150 (MQ=255)
                                                                                         cGGAATAGTTCTGGGCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTa                                                           <  2:1012385/150‑1 (MQ=255)
                                                                                            aaTAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATg                                                        <  2:1222270/150‑1 (MQ=255)
                                                                                            aaTAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCGTAATg                                                        <  1:923572/150‑1 (MQ=255)
                                                                                             aTAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGc                                                       >  2:414387/1‑150 (MQ=255)
                                                                                              tAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGcc                                                      <  1:1779169/150‑1 (MQ=255)
                                                                                               aGTTCTGGTCTGGCCCACGAATGACTAAGCGAACGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCa                                                     >  2:1491855/1‑150 (MQ=255)
                                                                                                  tctggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCt                                                  >  1:1721872/1‑150 (MQ=255)
                                                                                                   ctggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTc                                                 >  1:773075/1‑150 (MQ=255)
                                                                                                     ggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCg                                               >  1:1209909/1‑150 (MQ=255)
                                                                                                     ggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCg                                               >  2:153969/1‑150 (MQ=255)
                                                                                                     ggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCg                                               <  2:1773981/150‑1 (MQ=255)
                                                                                                       tctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGTGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTa                                             <  1:1913150/150‑1 (MQ=255)
                                                                                                           gCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAAc                                         <  1:898956/150‑1 (MQ=255)
                                                                                                           gCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAAc                                         <  1:892227/150‑1 (MQ=255)
                                                                                                             ccACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCa                                       >  2:1831558/1‑150 (MQ=255)
                                                                                                             ccACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCa                                       <  1:153969/150‑1 (MQ=255)
                                                                                                              cACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCAt                                      <  2:53883/150‑1 (MQ=255)
                                                                                                               aCGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCAtt                                     >  2:1219061/1‑150 (MQ=255)
                                                                                                                  aaTGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTAcc                                  <  1:1736967/150‑1 (MQ=255)
                                                                                                                   aTGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATACGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCa                                 >  2:1290350/1‑150 (MQ=255)
                                                                                                                    tGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCaa                                >  1:1716285/1‑150 (MQ=255)
                                                                                                                    tGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACCTAAATGCCAGCCCATGCTTAATGCAAGCTCCGAAAAACCATTACCaa                                >  1:846912/1‑150 (MQ=255)
                                                                                                                            cGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCTGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGctct                        <  1:185725/150‑1 (MQ=255)
                                                                                                                             gATCGATAAAGCATTATCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCACTCGCGCTTCGTTTTCCGCCACATCCGACCGCAAACGTAAATGCAAGACCAAGCTTAATGCCAGCTCCCTAAAACCATTACCAATCTGCTATg                       >  2:1863917/1‑150 (MQ=255)
                                                                                                                               tCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgc                     <  2:1188061/150‑1 (MQ=255)
                                                                                                                               tCGCTAAAGCAATCTCCCGGCTGCGGGGTGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgc                     <  2:1192386/150‑1 (MQ=255)
                                                                                                                                 gCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcag                   >  2:946882/1‑150 (MQ=255)
                                                                                                                                  tgctggCATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcagc                  >  2:735114/6‑150 (MQ=255)
                                                                                                                                   gctggCATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcagcc                 >  1:1940786/5‑149 (MQ=255)
                                                                                                                                    ctggCATTCTCCCCGCAGGCGGGAGAATGCTAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCAc                >  1:1885244/4‑150 (MQ=255)
                                                                                                                                    aaaGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAACCCATTACCAATCTGCTCTGGCAGCAc                >  1:1355876/1‑150 (MQ=255)
                                                                                                                                      aGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTa              <  2:1643694/150‑1 (MQ=255)
                                                                                                                                        cATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACg            >  2:263845/1‑150 (MQ=255)
                                                                                                                                         aTTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGc           <  2:39438/150‑1 (MQ=255)
                                                                                                                                          ttCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgca                    >  2:577265/1‑140 (MQ=255)
                                                                                                                                          ttCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgca                    <  1:577265/140‑1 (MQ=255)
                                                                                                                                            ctcCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCACGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGa        <  2:1968908/150‑1 (MQ=255)
                                                                                                                                              cccGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCAGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGaaa      >  2:1024767/1‑150 (MQ=255)
                                                                                                                                                  ccTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATc  <  2:537103/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATC  >  NC_000913/2711982‑2712277

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