Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,881,358 T→A 17.1% D151V (GAC→GTC)  dmlR ← transcriptional activator of dmlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,881,3580TA17.1% 282.1 / 50.0 117D151V (GAC→GTC) dmlRtranscriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base T (56/41);  new base A (10/10);  total (66/51)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 118 total aligned reads displayed.

AATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTC  >  NC_000913/1881260‑1881456
                                                                                                  |                                                                                                  
aaTGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTc                                                                                                   >  1:1664427/1‑100 (MQ=255)
aaTGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTc                                                                                                   >  2:532364/1‑100 (MQ=255)
 aTGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCa                                                                                                  >  2:159104/1‑100 (MQ=255)
  tGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTc                                                                                                   >  1:92949/1‑98 (MQ=255)
  tGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCAt                                                                                                 >  2:2888144/1‑100 (MQ=255)
       cGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCa                                                                                            >  1:2353995/1‑100 (MQ=255)
       cGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGccc                                                                                            <  1:606805/100‑2 (MQ=255)
        ggCGCGCCATGTTGCGCAATAAATTCCGGCGACTCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccc                                                                                           <  2:652638/100‑3 (MQ=255)
        ggCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccc                                                                                           <  2:2101872/100‑3 (MQ=255)
         gcgcgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccct                                                                                          <  2:410059/100‑4 (MQ=255)
          cgcgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGccccta                                                                                         <  2:2635825/100‑5 (MQ=255)
          cgcgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGGCGCCAGTTTGCGGGCGATAAAATTAGGGGCGATGACATCGccc                                                                                            <  1:607037/97‑2 (MQ=255)
            cgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCg                                                                                       <  2:2339152/100‑1 (MQ=255)
              ccATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGaa                                                                                     >  2:1545094/1‑100 (MQ=255)
              ccATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccctaattt                                                                                     <  1:745263/100‑9 (MQ=255)
               cATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAat                                                                                    <  1:2258098/100‑1 (MQ=255)
               cATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccctaattta                                                                                    <  2:2623667/100‑10 (MQ=255)
                aTGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAata                                                                                   >  1:492490/1‑100 (MQ=255)
                aTGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAata                                                                                   >  2:311760/1‑100 (MQ=255)
                 tGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAatat                                                                                  <  2:2602552/100‑1 (MQ=255)
                  gTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTAGCGGGCGAATAAATTAGGGGCGATGTCATCGCCAATGCGAATATc                                                                                 >  1:1752415/1‑100 (MQ=255)
                  gTTGCGCAATAAATTCCGGAGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATc                                                                                 >  2:669395/1‑100 (MQ=255)
                   ttGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCg                                                                                >  1:2871133/1‑100 (MQ=255)
                    tGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCga                                                                               >  1:725931/1‑100 (MQ=255)
                     gcgcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCgag                                                                              <  2:539948/100‑1 (MQ=255)
                      cgcAATAAATTGCGGCGACGCGCAGAGAATATGGTAATTGGTCGTCAGTTTGCTGGCGATGAAATTAGGGGCGATGACGTCGcccctaat                                                                                       <  1:1115071/90‑7 (MQ=255)
                      cgcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTACATTAGAGGCGATGTCATCGCCAATGCGAATATCgaga                                                                             >  1:1959418/1‑100 (MQ=255)
                      cgcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAAATAGGGGCGATGTCATCGCCAATGAGAATATCgaga                                                                             >  1:1727410/1‑100 (MQ=255)
                       gcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCCGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGAt                                                                            >  1:1728092/1‑100 (MQ=255)
                        cAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATc                                                                           >  2:197837/1‑100 (MQ=255)
                        cAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATc                                                                           <  2:696958/100‑1 (MQ=255)
                        cAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATc                                                                           >  2:896761/1‑100 (MQ=255)
                          aTAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCg                                                                          >  1:1740003/1‑99 (MQ=255)
                           tAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGac                                                                        <  1:2145414/100‑1 (MQ=255)
                             aaTTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGacac                                                                      >  1:1307666/1‑100 (MQ=255)
                              aTTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACAcc                                                                     >  1:1937076/1‑100 (MQ=255)
                               ttCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCt                                                                    >  1:1951991/1‑100 (MQ=255)
                                tCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTTATCGCCAATGCGAACATCGAGATCGACACCtt                                                                   <  1:63379/100‑1 (MQ=255)
                                tCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCtt                                                                   <  1:2158388/100‑1 (MQ=255)
                                 ccGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTc                                                                  <  2:1542572/100‑1 (MQ=255)
                                  cGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCa                                                                 >  2:917399/1‑100 (MQ=255)
                                  cGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCa                                                                 >  2:2030087/1‑100 (MQ=255)
                                   aacGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCAt                                                                <  2:623973/98‑1 (MQ=255)
                                    gcggcgCGCAGAGAATACGGTAATTGGTTCTCAGTTTGCGGGCGATTAAATTAGTGTCTATGTTAATGCTAAGGCGAATACCGAGATCGACACCTTCAtt                                                               <  1:1744206/100‑1 (MQ=255)
                                    gCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCAtt                                                               >  2:852462/1‑100 (MQ=255)
                                     cGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGTGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTg                                                              <  1:1027489/100‑1 (MQ=255)
                                      gACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGa                                                             >  2:2474020/1‑100 (MQ=255)
                                       aCGCTCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGTCGATGTCAATGATAATGCGAATACCGAGATCTACACCTTCATTGa                                                             <  1:391304/99‑1 (MQ=255)
                                        cgcgCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGAcc                                                           <  2:1781854/100‑1 (MQ=255)
                                         gcgcAGAGAATACGTGAACTGGTCGCCAGTTTGCTGGTGATTAAATTAGGGGCGATGTTATCGCCAATGCGAATATCGATATCGACACCTTCATTGAc                                                            <  1:2369399/98‑1 (MQ=255)
                                         gcgcAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCa                                                          <  1:933259/100‑1 (MQ=255)
                                         gcgcAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCa                                                          <  2:2501244/100‑1 (MQ=255)
                                         cttccgatcttACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCa                                                          <  2:2404140/90‑1 (MQ=255)
                                           gcAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCaaa                                                        >  2:2102933/1‑100 (MQ=255)
                                              gagaATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCCAGATCGACACCTTCATTGACCAAATCg                                                     >  1:2211568/1‑100 (MQ=255)
                                                gaATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGAc                                                   >  2:2239153/1‑100 (MQ=255)
                                                gaATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGAc                                                   >  2:923181/1‑100 (MQ=255)
                                                gaATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGAc                                                   >  2:93853/1‑100 (MQ=255)
                                                 aaTACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGAcc                                                  <  1:2888133/100‑1 (MQ=255)
                                                  aTACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCa                                                 >  2:2175021/1‑100 (MQ=255)
                                                   tACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCaagatcggaag                                                >  1:2404129/1‑91 (MQ=255)
                                                    aCGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGt                                               >  2:1327585/1‑100 (MQ=255)
                                                      ggTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCg                                             >  2:1071390/1‑100 (MQ=255)
                                                      ggTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATCGACCAAATCGACCAg                                                <  1:2814434/97‑1 (MQ=255)
                                                        tAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGt                                           <  1:720014/100‑1 (MQ=255)
                                                         aaTTGGTCGCCAGTTTGCGGGCGATTAAATTTTTTTTGTTTTTAATGATACGGCGACCACCGAGATCTACACCTTCATTGACCAAATCGACCAGTCGGTc                                          <  1:1951948/100‑1 (MQ=255)
                                                          aTTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCt                                         <  2:1317176/100‑1 (MQ=255)
                                                             ggTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTcc                                      <  2:758470/100‑1 (MQ=255)
                                                              gTCGCCAGTTTGCGGGCGATTAAATTAGGTGTTTTGTCATTGCTACGGCGAACATCGAGATCTACACCTTCATTGACCAAATCGACCAGTCGGTCTTCca                                     <  1:229980/100‑1 (MQ=255)
                                                                cGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCcaca                                   >  1:809547/1‑100 (MQ=255)
                                                                cGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCcaca                                   >  1:601514/1‑100 (MQ=255)
                                                                cGCCAGTTTGCGGACGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCcaca                                   >  1:1395073/1‑100 (MQ=255)
                                                                 gCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACAt                                  <  1:975400/100‑1 (MQ=255)
                                                                    aGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCaa                               <  2:1307668/100‑1 (MQ=255)
                                                                    aGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCaa                               <  1:1473510/100‑1 (MQ=255)
                                                                     gTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCaaa                              >  1:267582/1‑100 (MQ=255)
                                                                      tttGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACCTCaa                               >  1:334116/1‑98 (MQ=255)
                                                                        tGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACg                           <  2:1375373/100‑1 (MQ=255)
                                                                          cGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTa                         >  1:141667/1‑100 (MQ=255)
                                                                            ggCGATTAAATTAGGGGCGTTGTTATCGCCAATGCGACCATCGAGATCTACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAAt                       <  1:618532/100‑1 (MQ=255)
                                                                            ggCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTACACATCAAAACGTAAt                       <  1:1991391/100‑1 (MQ=255)
                                                                               gATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTcc                    >  1:779848/1‑100 (MQ=255)
                                                                                aTTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCa                   <  2:601515/100‑1 (MQ=255)
                                                                                    aaTTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTg               >  2:2721559/1‑100 (MQ=255)
                                                                                     cgcattggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCCAAACGTAATTCCAGtt                >  1:466333/7‑98 (MQ=255)
                                                                                      ttAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCg             <  2:55591/100‑1 (MQ=255)
                                                                                      ttAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCg             <  2:2822744/100‑1 (MQ=255)
                                                                                      gcattggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCg             >  2:337919/6‑100 (MQ=255)
                                                                                      gcattggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCg             >  2:2784664/6‑100 (MQ=255)
                                                                                       tAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCgg            <  2:290108/100‑1 (MQ=255)
                                                                                        attggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGa           >  2:1319665/4‑100 (MQ=255)
                                                                                         ttggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGAt          >  2:1263982/3‑100 (MQ=255)
                                                                                         ggggCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCg             >  1:214977/1‑97 (MQ=255)
                                                                                           ggCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATAc        >  2:649409/1‑100 (MQ=255)
                                                                                           ggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATAc        >  2:611255/1‑100 (MQ=255)
                                                                                           ggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATAc        >  2:1681152/1‑100 (MQ=255)
                                                                                               aTGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCt    >  2:1631689/1‑100 (MQ=255)
                                                                                                tGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCtt   >  2:1723225/1‑100 (MQ=255)
                                                                                                tGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCtt   >  1:2064380/1‑100 (MQ=255)
                                                                                                  tCAGCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTc  <  1:2216720/99‑1 (MQ=255)
                                                                                                  |                                                                                                  
AATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTC  >  NC_000913/1881260‑1881456

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: