Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,881,128 C→A E228* (GAA→TAA)  dmlR ← transcriptional activator of dmlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,881,1280CA98.4% 399.8 / ‑6.5 128E228* (GAA→TAA) dmlRtranscriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (72/53);  total (74/54)
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 130 total aligned reads displayed.

ACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTT  >  NC_000913/1881031‑1881224
                                                                                                 |                                                                                                
aCTTGCACTAAATGACCACTGGCAATGTTTTCGCGAACATCCCACCATGATCGTAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTtatc                                                                                                <  1:1454423/100‑1 (MQ=255)
 cTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATcc                                                                                               >  1:2264807/1‑100 (MQ=255)
    gCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCTCATCGAGACACCACTGATGGACAATTTATCCATg                                                                                            <  1:1453357/100‑1 (MQ=255)
    gCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATg                                                                                            <  2:3467905/100‑1 (MQ=255)
     cACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGa                                                                                           >  1:1538672/1‑100 (MQ=255)
       cTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCGCTGCCCATCGAGACACCACTGATGGACAATTTATCCa                                                                                              >  1:677840/1‑95 (MQ=255)
        tAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTa                                                                                        <  1:25159/100‑1 (MQ=255)
        tAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTa                                                                                        <  1:456458/100‑1 (MQ=255)
        tAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTa                                                                                        >  2:692563/1‑100 (MQ=255)
           aTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgac                                                                                     <  2:2883181/100‑1 (MQ=255)
            tGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgacg                                                                                    <  2:2598852/100‑1 (MQ=255)
             gACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgacga                                                                                   >  2:765996/1‑100 (MQ=255)
              aCCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgacgac                                                                                  <  2:2750685/100‑1 (MQ=255)
               ctaCTGGCTATGTTTTCGCTTACATCCCACCAGGAACGCAGAGCAATTCCCTGCCGATCGTGACACCATTGGTGGTCACTTTATCCATGATTAg                                                                                       >  1:2869337/3‑94 (MQ=255)
               ccACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACa                                                                                 <  1:174174/100‑1 (MQ=255)
               ccACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACa                                                                                 >  1:241838/1‑100 (MQ=255)
               ccACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATATATCCATGATTAgacgac                                                                                  >  1:2221786/1‑99 (MQ=255)
                  cTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCACCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACATTg                                                                              >  1:2950807/1‑100 (MQ=255)
                    ggCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTggg                                                                            >  2:2437262/1‑100 (MQ=255)
                     gCAATGTTTTCGCTAACATCCCACCCGCAACGCAGAGCAATCCCCTGCCCACCGAGACACCACTGACGTCCAATTTATCCATGATTAGACGACAGt                                                                               >  1:900607/1‑96 (MQ=255)
                     gCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGc                                                                           >  1:900577/1‑100 (MQ=255)
                      cAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGcc                                                                          >  2:2735525/1‑100 (MQ=255)
                      cAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGcc                                                                          <  2:1881351/100‑1 (MQ=255)
                         tGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCTTTGATTAGACGACACTGGGc                                                                           >  1:1033252/1‑96 (MQ=255)
                          gTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTg                                                                      <  2:3117037/100‑1 (MQ=255)
                          gTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTg                                                                      <  2:3112541/100‑1 (MQ=255)
                          gTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTg                                                                      <  2:2711387/100‑1 (MQ=255)
                           ttttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGt                                                                       >  1:3464140/1‑98 (MQ=255)
                            tttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGAc                                                                    <  1:1263981/100‑1 (MQ=255)
                            tttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGAc                                                                    <  1:3270760/100‑1 (MQ=255)
                              tCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACtt                                                                  <  1:351201/100‑1 (MQ=255)
                              tCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACtt                                                                  <  2:2518229/100‑1 (MQ=255)
                               cGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACttt                                                                 >  1:2065600/1‑100 (MQ=255)
                               cGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACttt                                                                 >  2:3088880/1‑100 (MQ=255)
                               cGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACtt                                                                  >  2:1495554/1‑99 (MQ=255)
                               cGCTAACATCCCACCAGGAACGCAAAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAATGGGCAGGAGACtt                                                                  >  1:3146328/1‑99 (MQ=255)
                                gCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCCCGATTAGACGACAGTGGGACGGGGACttt                                                                 >  1:360732/1‑99 (MQ=255)
                                 cTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTaa                                                               <  2:837221/100‑1 (MQ=255)
                                   aaCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATc                                                             <  2:924013/100‑1 (MQ=255)
                                     cATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGCTGGACAATCTATCCATGATTAGACGACAGCGGGCCGGTGACTTTAATCGc                                                           >  1:1224386/1‑100 (MQ=255)
                                       tCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCAt                                                         <  1:901132/100‑1 (MQ=255)
                                        cccACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATg                                                        >  2:1034229/1‑100 (MQ=255)
                                        cccACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATg                                                        <  2:757617/100‑1 (MQ=255)
                                           accaGGAACGCAGAGCAATCCCCTGCACATCGAGACACCACCGGTGGACAATTTATCCAGGCCTAGACGACCGTAGGCCGACGACTTGaacgcgataggg                                                     >  1:1190447/1‑90 (MQ=255)
                                            ccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATggggg                                                    <  1:1755956/100‑1 (MQ=255)
                                            ccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATggggg                                                    >  2:1895670/1‑100 (MQ=255)
                                             caGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGc                                                   >  2:2178930/1‑100 (MQ=255)
                                              aGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTACTCACATggg                                                      >  1:198300/1‑96 (MQ=255)
                                               ggAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAACTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCct                                                 >  2:2181130/1‑100 (MQ=255)
                                                 aaCGCGAAGCAAACCCCTGCCCAGCGAGACACCACTGATGGACAATTTCGCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCcttc                                               <  2:2377440/100‑1 (MQ=255)
                                                 aaCGCAGAGCAATCCCCTGCCCATCGAGACAGCACTAGCGGCGAGTATATCCAGGATTAGCCGCCAATTGGCCGGTGACTCTCACCGCATGCAGTCcttc                                               >  1:1466467/1‑100 (MQ=255)
                                                 aaCGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGGTGGTGACTTTAATCACATGGGGGGctt                                                >  1:487431/1‑99 (MQ=255)
                                                  aCGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATgggggagatcg                                              >  1:2506191/1‑94 (MQ=255)
                                                  aCGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCcttct                                              >  2:958668/1‑100 (MQ=255)
                                                   cGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCcttctt                                             >  2:3334441/1‑100 (MQ=255)
                                                     cAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTa                                           >  2:981227/1‑100 (MQ=255)
                                                       gagCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTa                                           >  1:113210/1‑98 (MQ=255)
                                                       gagCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTAtt                                         >  2:1277507/1‑100 (MQ=255)
                                                        tgCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                        <  2:1178510/99‑1 (MQ=255)
                                                        agCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                        <  1:2422524/100‑1 (MQ=255)
                                                        agCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                        <  2:1564391/100‑1 (MQ=255)
                                                        agCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                        <  1:2822901/100‑1 (MQ=255)
                                                         gaaaTCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGc                                       >  2:2344630/3‑100 (MQ=255)
                                                         gCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTga                                       >  1:1178495/1‑99 (MQ=255)
                                                         gCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCATCTTTATTg                                        >  1:422598/1‑99 (MQ=255)
                                                           aaTCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGt                                     <  2:2980459/100‑1 (MQ=255)
                                                             tCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTaa                                   <  1:2851609/100‑1 (MQ=255)
                                                              ccccTGCCCATCGAGATACCACTGATGGACAAGTTATGTATGATTAGACGACAGTGGGCCGGTTACTTTAATTGAATGGAGGACTTCTTTATGGCg                                      >  1:2254447/1‑96 (MQ=255)
                                                              ccccTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAAt                                  >  1:2143790/1‑100 (MQ=255)
                                                              ccccTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGAGCCTTCTTTATTGCGTAAt                                  >  1:68281/1‑100 (MQ=255)
                                                                ccTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTg                                <  1:962702/100‑1 (MQ=255)
                                                                 cTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGc                               >  2:1627605/1‑100 (MQ=255)
                                                                  tGCCCATCGAGACACCACTGATGGACAATTTGTCCCTGATTAGACGACAGTGGGCGGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGc                               >  1:838879/1‑99 (MQ=255)
                                                                    cccATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAg                            <  2:2143821/100‑1 (MQ=255)
                                                                      cATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAg                            >  1:1414575/1‑98 (MQ=255)
                                                                          gagaCACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCa                      <  2:2823755/100‑1 (MQ=255)
                                                                          gagaCACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCATAATTGCCAGACACCa                      <  1:2685364/100‑1 (MQ=255)
                                                                           agaCACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCaa                     >  2:706549/1‑100 (MQ=255)
                                                                            gaCACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCaaa                    >  2:1789802/1‑100 (MQ=255)
                                                                              caccacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAg                  >  1:520150/1‑100 (MQ=255)
                                                                              caccacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAg                  <  1:1470226/100‑1 (MQ=255)
                                                                                ccacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGga                >  2:2986663/1‑100 (MQ=255)
                                                                                 cacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGgat               <  1:1475442/100‑1 (MQ=255)
                                                                                   cTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGTCAGACACCAAAAGgatga             <  2:241841/100‑1 (MQ=255)
                                                                                   cTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGgatga             >  2:685089/1‑100 (MQ=255)
                                                                                     gATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATc           <  1:331094/100‑1 (MQ=255)
                                                                                      aTGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcg          >  1:2233869/1‑100 (MQ=255)
                                                                                       tGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgc         <  1:3035127/100‑1 (MQ=255)
                                                                                        ggACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgcg        >  2:1868009/1‑100 (MQ=255)
                                                                                         gACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCgcgc       >  2:521191/1‑100 (MQ=255)
                                                                                         gACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCgcgc       >  1:3276864/1‑100 (MQ=255)
                                                                                         gACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCgcgc       >  1:464297/1‑100 (MQ=255)
                                                                                           cAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctc     >  1:1922293/1‑100 (MQ=255)
                                                                                           cAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctc     >  2:1385388/1‑100 (MQ=255)
                                                                                           cAATTTATCCACGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctc     >  1:2968860/1‑100 (MQ=255)
                                                                                            aaTTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctct    >  2:1337826/1‑100 (MQ=255)
                                                                                             aTTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCtt   <  2:3470897/100‑1 (MQ=255)
                                                                                             aTTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCtt   >  2:2918024/1‑100 (MQ=255)
                                                                                              tttATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCttt  <  1:1808372/100‑1 (MQ=255)
                                                                                              tttATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCttt  <  1:1248800/100‑1 (MQ=255)
                                                                                                 |                                                                                                
ACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTT  >  NC_000913/1881031‑1881224

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