Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,466,211 G→A 16.9% intergenic (‑31/+88) treB ← / ← treR trehalose‑specific PTS enzyme: IIB and IIC component/trehalose 6‑phosphate‑inducible trehalose regulon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,466,2110GA16.9% 257.2 / 37.9 118intergenic (‑31/+88)treB/treRtrehalose‑specific PTS enzyme: IIB and IIC component/trehalose 6‑phosphate‑inducible trehalose regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (54/41);  new base A (12/8);  total (68/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 129 total aligned reads displayed.

AATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAGGACAGGGT  >  NC_000913/4466112‑4466310
                                                                                                   |                                                                                                   
aaTATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTttggttg                                                                                                     <  1:139334/100‑1 (MQ=255)
aaTATTGCCGCGCCCGCCGACCAGTTCAAGCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTttggttg                                                                                                     <  1:771248/100‑1 (MQ=255)
  ccgatctgcgcCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGc                                                                                                   <  1:10727/93‑1 (MQ=255)
    ttGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTg                                                                                                 <  2:2155691/100‑1 (MQ=255)
    ttGCCCCGCCCGCCCACCAGTTCAATCACCCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGACGACATTTttggttg                                                                                                     <  1:2460011/96‑1 (MQ=255)
     tGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTAGCTGc                                                                                                >  1:1515173/1‑100 (MQ=255)
     tGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCAACATAAAGCCCCATGTCAGATGACATTTTTGGTTGGCTg                                                                                                 >  1:390962/1‑99 (MQ=255)
      gCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCa                                                                                               >  2:86644/1‑100 (MQ=255)
      gCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCa                                                                                               >  2:1600993/1‑100 (MQ=255)
        cgcgCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGa                                                                                             >  2:2036102/1‑100 (MQ=255)
         gcgcCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGcagatcgg                                                                                            >  2:10727/1‑93 (MQ=255)
           ctccgccgACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAata                                                                                          <  2:203405/98‑1 (MQ=255)
            cccgccGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAatat                                                                                         <  2:235440/100‑1 (MQ=255)
            cccgccGACCAGTTCAATCAACCAATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAata                                                                                          <  1:1573223/99‑1 (MQ=255)
             ccgccgACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAataag                                                                                        >  1:203403/1‑98 (MQ=255)
                  gACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAat                                                                                           <  1:2534032/92‑1 (MQ=255)
                   aCCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTTTTTGGCTGTAGAATACATTgcg                                                                                   >  1:485301/1‑99 (MQ=255)
                   aCCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTgcgc                                                                                  <  2:1891746/100‑1 (MQ=255)
                   aCCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTAGCTGCAGAATATATTgcgc                                                                                  <  2:2200178/100‑1 (MQ=255)
                     cAGTTCAATCAACCGATCGATCTCCGTTTGGTTTATTTTGCTCATCATAAGGCCCCATGGCAGATGACATTTTTGGTTGGCTGGAGAATATATTGCGCgg                                                                                >  1:2303907/1‑100 (MQ=255)
                      aGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTAGCTGCAGAATATATTGCGCGGa                                                                               >  1:39763/1‑100 (MQ=255)
                       gTTCAATCAACCGATGGATATCCGTTTGGGTAATTTTGCTCATCATAAAGCCCCATGGCAGGTCACATTTTTGGTTGGCTGCAGAATATATTGCGCGGaa                                                                              >  1:2369754/1‑100 (MQ=255)
                       gTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGaa                                                                              >  1:934560/1‑100 (MQ=255)
                       gTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGaa                                                                              >  2:1242975/1‑100 (MQ=255)
                       gTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGaa                                                                              >  2:146854/1‑100 (MQ=255)
                          caatcaaCCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGaaaaa                                                                           >  2:1518074/1‑100 (MQ=255)
                          caatcaaCCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGaaaaa                                                                           >  2:1174763/1‑100 (MQ=255)
                          caatcaaCCGATCGATATCCGTTTGGTTTAATTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGaaaaa                                                                           >  2:1669108/1‑100 (MQ=255)
                            atcaaCCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAAtt                                                                         >  1:145565/1‑100 (MQ=255)
                            atcaaCCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAAtt                                                                         <  1:1296821/100‑1 (MQ=255)
                             tcaaCCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAAttt                                                                        >  2:2377803/1‑100 (MQ=255)
                                aCCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTaaa                                                                     >  1:1873110/1‑100 (MQ=255)
                                   gatcgatATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAAc                                                                  <  2:529667/100‑1 (MQ=255)
                                   gatcgatATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAAGTTaaaaa                                                                   >  2:1677238/1‑99 (MQ=255)
                                   gatcgatATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTAGCTGCAGAATATATTGCGCGGAAAAATTTAAAAAc                                                                  >  2:2514553/1‑100 (MQ=255)
                                    atcgatATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTAGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACg                                                                 <  1:1638854/100‑1 (MQ=255)
                                      cgatATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACggg                                                               <  1:689140/100‑1 (MQ=255)
                                      cgatATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACggg                                                               >  2:513330/1‑100 (MQ=255)
                                      cgatATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACggg                                                               >  1:240595/1‑100 (MQ=255)
                                       gatATCCGTTTGGTTTATTTTGCTCATCATAAAACCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGa                                                              >  1:1553199/1‑100 (MQ=255)
                                           tCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGt                                                          >  1:1968832/1‑100 (MQ=255)
                                           tCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGt                                                          >  1:1305175/1‑100 (MQ=255)
                                            ccGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGtt                                                         >  1:980315/1‑100 (MQ=255)
                                             cGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTc                                                        >  1:1734053/1‑100 (MQ=255)
                                              gTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTcc                                                       >  2:539385/1‑100 (MQ=255)
                                              gTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTcc                                                       <  1:299530/100‑1 (MQ=255)
                                                ttGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCg                                                     <  1:1711906/100‑1 (MQ=255)
                                                ttGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTCCGCGGAAAAATTTAAAAACGGGAACGTTCCCg                                                     <  2:2457842/100‑1 (MQ=255)
                                                 tGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGTTTATCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGa                                                    >  1:307843/1‑100 (MQ=255)
                                                  ggTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTTTTTTTTTTTAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGaa                                                   >  1:406601/1‑100 (MQ=255)
                                                  ggTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTAGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGaa                                                   <  1:584064/100‑1 (MQ=255)
                                                   gTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGATATTTTTTTTTTTTTGTAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGaaa                                                  >  1:1382484/1‑100 (MQ=255)
                                                   gTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGaaa                                                  >  1:1760586/1‑100 (MQ=255)
                                                   gTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGaaa                                                  <  1:819790/100‑1 (MQ=255)
                                                    tttATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTAGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAAc                                                 >  1:417076/1‑100 (MQ=255)
                                                       aTTTTGCTCATCATAAAGCCCCATAGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCa                                              >  1:453066/1‑100 (MQ=255)
                                                         tttGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGc                                            <  2:94792/100‑1 (MQ=255)
                                                           tGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGa                                          <  1:1001852/100‑1 (MQ=255)
                                                           tGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTAGCTGCAGAATATCTTGCGCGGAAAAATTTAAAACCGGGAACGTTCCCGAAACGCAgaga                                          >  2:2248775/1‑100 (MQ=255)
                                                           tGCTCAGCATAAAGCGCCACGGCAGATGAAATTTTTGATTGCCTGCAGAATAAAATGCACGGAAAAATGAAAAAACGGATAACTTCCCGAAACGCAGCg                                           >  1:976012/1‑99 (MQ=255)
                                                            gCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGaa                                         >  2:715448/1‑100 (MQ=255)
                                                             cTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAg                                        <  2:2330706/100‑1 (MQ=255)
                                                             cTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAg                                        >  1:597420/1‑100 (MQ=255)
                                                               catcatAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGAt                                      >  1:1594676/1‑100 (MQ=255)
                                                               catcatAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGAt                                      <  2:2437480/100‑1 (MQ=255)
                                                               catcatAAAGCCCCATGGCAGATGACATTTTTGGTTAGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGAt                                      >  2:1225571/1‑100 (MQ=255)
                                                                atcatAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATc                                     <  1:1181787/100‑1 (MQ=255)
                                                                 tcatAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATca                                    >  2:1264970/1‑100 (MQ=255)
                                                                   atAAAGCCCCATGGCAGATGACATTTTTGGTTAGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATcaca                                  <  1:1153497/100‑1 (MQ=255)
                                                                   atAAAGCCCCATGCCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATcaca                                  >  2:1760909/1‑100 (MQ=255)
                                                                    tAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACaa                                 >  1:676622/1‑100 (MQ=255)
                                                                    tAAAGCCCCATGGCAGATGACATTTTTGGTTAGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACaa                                 >  1:744483/1‑100 (MQ=255)
                                                                     aaaGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAAt                                >  2:2508902/1‑100 (MQ=255)
                                                                      aaGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAAtt                               >  2:1645018/1‑100 (MQ=255)
                                                                       aGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAAtt                               <  2:482530/99‑1 (MQ=255)
                                                                         ccccATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTAt                            <  1:1959671/100‑1 (MQ=255)
                                                                          cccATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATc                           >  2:1039663/1‑100 (MQ=255)
                                                                          cccATGGCAGATGACATTTTTGGTTAGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATc                           <  2:725128/100‑1 (MQ=255)
                                                                           ccATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCg                          <  2:1361237/100‑1 (MQ=255)
                                                                            cATGGCAGATGACATTTTTGGTTAGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGt                         <  2:1382492/100‑1 (MQ=255)
                                                                              tGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGt                         <  1:592649/98‑1 (MQ=255)
                                                                                gCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAg                     <  1:1163280/100‑1 (MQ=255)
                                                                                gCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAg                     <  1:415107/100‑1 (MQ=255)
                                                                                 cAGATGACATTCTTGGTTGGTTGCAGAATATATTGAGCGTAAACATTTAAAAACGGGAACTATCCAGAAAAGCAGCGGACATCACAATTTATCg                          >  1:1258476/1‑94 (MQ=255)
                                                                                     tGACATTTTTGTTTGGCTGCAAAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAAc                <  1:694118/100‑1 (MQ=255)
                                                                                     tGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAAc                <  1:1421821/100‑1 (MQ=255)
                                                                                      gACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACg               >  2:747632/1‑100 (MQ=255)
                                                                                       aCATTTTTGATTTGCTACAGAATATATTGCTCAGATAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAATATCACAATTTATCGTTCAGGAAACGa              <  1:1272639/100‑1 (MQ=255)
                                                                                          tttttGGTTGGCTACAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCa           <  2:409093/100‑1 (MQ=255)
                                                                                           ttttGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAg          <  2:2235673/100‑1 (MQ=255)
                                                                                            tttGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAgg         <  1:2322699/100‑1 (MQ=255)
                                                                                            tttGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAgg         >  1:889501/1‑100 (MQ=255)
                                                                                             ttggttggCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACTCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAGGa        <  2:1181127/100‑1 (MQ=255)
                                                                                             ttggttggCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAGGa        <  1:2457045/100‑1 (MQ=255)
                                                                                                nttggCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAGGACAg     >  2:356487/2‑100 (MQ=255)
                                                                                                gttagCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAGGACAg     >  1:657679/1‑100 (MQ=255)
                                                                                                gttagCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAGGACAg     >  2:1110016/1‑100 (MQ=255)
                                                                                                  tggCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAGGACAggg   >  2:934841/1‑100 (MQ=255)
                                                                                                   ggCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAGGACAGGGt  >  2:325061/1‑100 (MQ=255)
                                                                                                   ggCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAGGACAGGGt  >  1:1561587/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTATCGTTCAGGAAACGATCAGGACAGGGT  >  NC_000913/4466112‑4466310

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