Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,992,894 A→G T175A (ACG→GCG)  dapF → diaminopimelate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,992,8940AG99.7% 1184.6 / ‑5.0 367T175A (ACG→GCG) dapFdiaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (166/200);  total (166/201)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

Only 100 of 390 total aligned reads displayed.

GCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGT  >  NC_000913/3992804‑3992983
                                                                                          |                                                                                         
 ccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAgcgc                                                                                        <  1:94303/93‑1 (MQ=255)
   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                                                      <  1:4620866/93‑1 (MQ=255)
   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                                                      >  2:3928603/1‑93 (MQ=255)
   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                                                      >  2:5426728/1‑93 (MQ=255)
   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                                                      >  1:3884965/1‑93 (MQ=255)
   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                                                      <  2:471862/93‑1 (MQ=255)
    aGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTg                                                                                     >  1:7634105/1‑93 (MQ=255)
     gCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTgg                                                                                    >  2:9654087/1‑93 (MQ=255)
     gCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTgg                                                                                    >  2:2658776/1‑93 (MQ=255)
     gCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTgg                                                                                    >  2:6332394/1‑93 (MQ=255)
      cAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTggc                                                                                   >  2:1047809/1‑92 (MQ=255)
      cAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGt                                                                                   <  1:7462092/93‑1 (MQ=255)
        gACAATCTTATGCGGCGTGGTGTCGATGGGCAATCCGCAATGCGTGATTCAGGTCGAGGAGGTCGATACCGCGGCGGTAGAAGAGCTTGGtac                                                                                 >  1:4452395/1‑91 (MQ=255)
        gACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTcc                                                                                 >  1:10017818/1‑93 (MQ=255)
           aaTATTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTACTGt                                                                              >  2:5007178/1‑93 (MQ=255)
            aTCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGtt                                                                             >  1:7870590/1‑93 (MQ=255)
             tCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTc                                                                            >  2:740794/1‑93 (MQ=255)
               ttATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTg                                                                          >  2:127166/1‑93 (MQ=255)
                tATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTgg                                                                         >  2:7063822/1‑93 (MQ=255)
                  tGCGGCGTGGTGGCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGaa                                                                       <  1:3411203/93‑1 (MQ=255)
                   gcggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGCaaa                                                                      >  2:2304054/1‑93 (MQ=255)
                     ggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGc                                                                    <  1:2560670/93‑1 (MQ=255)
                      gcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGcc                                                                   >  2:8082675/1‑93 (MQ=255)
                      gcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGcc                                                                   <  2:5257412/93‑1 (MQ=255)
                       cgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCa                                                                  >  2:3398991/1‑93 (MQ=255)
                        gtggtgTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCAc                                                                 <  2:9773959/93‑1 (MQ=255)
                         tggtgTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACg                                                                >  1:9252180/1‑93 (MQ=255)
                          ggtgTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGa                                                               <  1:8443543/93‑1 (MQ=255)
                            tgtCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGc                                                             >  2:2978512/1‑93 (MQ=255)
                             gtCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCg                                                            <  2:1802398/93‑1 (MQ=255)
                             gtCGATGGGAAAGCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCg                                                            <  2:9886373/93‑1 (MQ=255)
                              tCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGt                                                           >  1:4111648/1‑93 (MQ=255)
                               cGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGtt                                                          <  1:6381209/93‑1 (MQ=255)
                                gATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGATAGAAGCGCTTGGACCTGTTCTGGAAAGCCACGAGCGttt                                                         >  2:3028315/1‑93 (MQ=255)
                                  tGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTc                                                       >  1:4076767/1‑93 (MQ=255)
                                   gggAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTcc                                                      >  2:8194838/1‑93 (MQ=255)
                                   gggAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTcc                                                      <  2:3589532/93‑1 (MQ=255)
                                    ggAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCg                                                     >  1:8520584/1‑93 (MQ=255)
                                    ggAAACGCGCATTGCGTGATTCAGGTCGATGATGTCGATGCCGCGGCGGTAGAAGCCGTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCg                                                     <  2:2458173/93‑1 (MQ=255)
                                     gAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCgg                                                    <  2:5714150/93‑1 (MQ=255)
                                       aaTCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAg                                                  >  1:2820845/1‑93 (MQ=255)
                                        aTCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgc                                                 >  2:7714838/1‑93 (MQ=255)
                                         tCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcg                                                >  1:3504766/1‑93 (MQ=255)
                                         tCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcg                                                >  1:3181713/1‑93 (MQ=255)
                                         tCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcg                                                <  2:5324718/93‑1 (MQ=255)
                                          ccGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgc                                               <  2:8716387/93‑1 (MQ=255)
                                          ccGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgc                                               >  1:9481628/1‑93 (MQ=255)
                                            gCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgcgc                                             <  2:908713/93‑1 (MQ=255)
                                            gCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgcgc                                             >  2:1886723/1‑93 (MQ=255)
                                            gCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgcgc                                             <  1:5394862/93‑1 (MQ=255)
                                              aTTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCa                                           >  1:507947/1‑93 (MQ=255)
                                                tGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAat                                         >  2:6941685/1‑93 (MQ=255)
                                                 gCGTGATTCAGGTCGATGATGTCGATCCCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAata                                        <  2:8603423/93‑1 (MQ=255)
                                                 gCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAata                                        <  1:5936080/93‑1 (MQ=255)
                                                 gCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAata                                        <  1:2433413/93‑1 (MQ=255)
                                                  cGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAatat                                       >  1:9857595/1‑93 (MQ=255)
                                                  cGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAatat                                       >  2:1583901/1‑93 (MQ=255)
                                                  cGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAatat                                       >  2:4422454/1‑93 (MQ=255)
                                                  cGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAatat                                       <  2:3470186/93‑1 (MQ=255)
                                                   gTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATc                                      >  1:5748553/1‑93 (MQ=255)
                                                   gTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATc                                      <  2:6618732/93‑1 (MQ=255)
                                                      aTTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGt                                   >  2:5993827/1‑93 (MQ=255)
                                                      aTTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGt                                   >  1:1149371/1‑93 (MQ=255)
                                                      aTTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGt                                   <  1:7144598/93‑1 (MQ=255)
                                                       ttCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGtt                                  >  2:8533919/1‑93 (MQ=255)
                                                       ttCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGtt                                  >  2:1721673/1‑93 (MQ=255)
                                                        tCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGttt                                 <  2:7629284/93‑1 (MQ=255)
                                                        tCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGttt                                 <  2:555928/93‑1 (MQ=255)
                                                          aGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTa                               <  1:2055009/93‑1 (MQ=255)
                                                          aGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTa                               <  1:1274700/93‑1 (MQ=255)
                                                          aGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTa                               <  1:6913587/93‑1 (MQ=255)
                                                          aGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTa                               <  2:4559431/93‑1 (MQ=255)
                                                            gTCGATGATGTCGATACCGCGGCGGTAGAAGCGATTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATg                             <  1:9731334/93‑1 (MQ=255)
                                                              cGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCa                           <  2:3722527/93‑1 (MQ=255)
                                                               gatgatGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCaa                          <  1:6839755/93‑1 (MQ=255)
                                                               gatgatGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCaa                          <  1:9997583/93‑1 (MQ=255)
                                                                atgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAg                         >  2:1384910/1‑93 (MQ=255)
                                                                atgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAg                         <  1:3262393/93‑1 (MQ=255)
                                                                  tatgTTGATACCGCGGCTGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATTCAAGTg                       <  1:5565376/92‑1 (MQ=255)
                                                                   atgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTgg                      >  1:8419883/1‑93 (MQ=255)
                                                                    tgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGt                     >  1:4953346/1‑93 (MQ=255)
                                                                      tCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTa                   >  1:6592210/1‑93 (MQ=255)
                                                                       cGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTaa                  <  1:5173947/93‑1 (MQ=255)
                                                                       cGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTaa                  >  1:1499673/1‑93 (MQ=255)
                                                                         atgccGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAgc                <  1:6209899/90‑1 (MQ=255)
                                                                          tACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAgcg               <  1:4274676/93‑1 (MQ=255)
                                                                           aCCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAgcgc              >  2:4037191/1‑93 (MQ=255)
                                                                            ccGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGcgcg             <  2:9486694/93‑1 (MQ=255)
                                                                            ccGCGGCGGGAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGAGCGCCAATATCGGGTTTATGCAAGTGGTTAAGcgcg             >  2:5888217/1‑93 (MQ=255)
                                                                             cgcgGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGa            <  2:6424695/93‑1 (MQ=255)
                                                                              gcggcgGTAGAAGCGCTTGGTCCTGTTCTGGAAATCCACGAGCGTTTTCCTGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAg           <  1:8910141/93‑1 (MQ=255)
                                                                                ggcggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCa         <  1:2634425/93‑1 (MQ=255)
                                                                                  cggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCata       <  1:1722200/93‑1 (MQ=255)
                                                                                  cggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCata       <  2:8278314/93‑1 (MQ=255)
                                                                                  cggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCata       <  1:7350389/93‑1 (MQ=255)
                                                                                  cggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCata       <  2:4895219/93‑1 (MQ=255)
                                                                                    gTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATAtt     >  1:8235684/1‑93 (MQ=255)
                                                                                     tAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTc    <  1:4859209/93‑1 (MQ=255)
                                                                                     tAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTc    >  2:3272685/1‑93 (MQ=255)
                                                                                      aGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCg   >  2:2829819/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGT  >  NC_000913/3992804‑3992983

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