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,8940AG98.7% 492.9 / ‑3.2 155T175A (ACG→GCG) dapFdiaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (67/86);  total (67/88)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 161 total aligned reads displayed.

GCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGT  >  NC_000913/3992804‑3992983
                                                                                          |                                                                                         
gccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCg                                                                                         >  2:4513539/1‑93 (MQ=255)
 ccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAgcgc                                                                                        <  2:3609130/93‑1 (MQ=255)
   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                                                      <  2:1211153/93‑1 (MQ=255)
      cAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGt                                                                                   <  1:4501971/93‑1 (MQ=255)
       aGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTc                                                                                  >  1:1993681/1‑93 (MQ=255)
       aGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTc                                                                                  >  1:1500351/1‑93 (MQ=255)
        gACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCTGCGGTAGAAGAGCTTGGTcc                                                                                 >  2:3996821/1‑93 (MQ=255)
        gACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTcc                                                                                 >  1:2820877/1‑93 (MQ=255)
        gACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTcc                                                                                 >  2:583754/1‑93 (MQ=255)
        gACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTcc                                                                                 >  2:1989917/1‑93 (MQ=255)
        aaCAATCTTATGCGGCGTGGTGTCGATGGAAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGCTAGAAACGCTTGGTca                                                                                 <  2:4508566/92‑2 (MQ=255)
         gcgaTCTTGTGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGTTGATGTCGATCCCGCGGCGGTAGACGCGCTTGGTCCt                                                                                <  2:4210099/90‑1 (MQ=255)
         aCAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCt                                                                                >  1:894861/1‑93 (MQ=255)
         aCAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCt                                                                                <  2:2189126/93‑1 (MQ=255)
         aCAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGCCCt                                                                                <  2:791790/93‑1 (MQ=255)
          cAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTg                                                                               >  2:1309216/1‑93 (MQ=255)
             tCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTc                                                                            <  2:4725669/93‑1 (MQ=255)
             tCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTc                                                                            <  2:862375/93‑1 (MQ=255)
             tCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTc                                                                            <  2:4258345/93‑1 (MQ=255)
             tCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTc                                                                            >  1:311427/1‑93 (MQ=255)
              cTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCt                                                                           <  1:1825022/93‑1 (MQ=255)
               ttATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTg                                                                          <  2:4192967/93‑1 (MQ=255)
                   gcggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGaaa                                                                      <  1:4007679/93‑1 (MQ=255)
                   gcggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGaaa                                                                      >  1:3684053/1‑93 (MQ=255)
                   gcggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGAGGTAGAAGCGCTTGGTCCTGTTCTGGaaa                                                                      <  1:378669/93‑1 (MQ=255)
                     ggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGc                                                                    <  1:948037/93‑1 (MQ=255)
                      gcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGcc                                                                   <  1:1027162/93‑1 (MQ=255)
                      gcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGcc                                                                   <  1:3310911/93‑1 (MQ=255)
                       cgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCa                                                                  <  1:2839563/93‑1 (MQ=255)
                       cgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCa                                                                  >  2:3197648/1‑93 (MQ=255)
                       cgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCa                                                                  >  1:2781387/1‑93 (MQ=255)
                        gtggtgTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCAc                                                                 <  2:3185738/93‑1 (MQ=255)
                        gtggtgTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCAc                                                                 <  1:336644/93‑1 (MQ=255)
                         tggtgTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACg                                                                <  1:2113989/93‑1 (MQ=255)
                         gggtgTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCTATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACg                                                                <  2:3049852/92‑1 (MQ=255)
                            ggtCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGc                                                             <  1:2847398/92‑1 (MQ=255)
                              tCGATGGGCAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGt                                                           >  2:199981/1‑93 (MQ=255)
                              gcGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGt                                                           <  2:4071339/92‑1 (MQ=255)
                                gATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGttt                                                         >  2:4145336/1‑93 (MQ=255)
                                gATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGttt                                                         >  2:2191936/1‑93 (MQ=255)
                                gATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCACGGCGGTAGAAGCGCTTGGACCTGTTATGGAAAGCCACGAGCGttt                                                         >  1:3444825/1‑93 (MQ=255)
                                 aTGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGtttt                                                        >  2:1510822/1‑93 (MQ=255)
                                  tGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTCGAAGAGCTTGGTCCTGTTCTGCAAAGCCACGAGCGTTTTc                                                       >  1:1182556/1‑93 (MQ=255)
                                  tGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTc                                                       >  1:4276152/1‑93 (MQ=255)
                                  tGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTc                                                       >  2:2635072/1‑93 (MQ=255)
                                  tGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTGTTc                                                       >  2:258445/1‑93 (MQ=255)
                                   gggAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTcc                                                      >  2:4396357/1‑93 (MQ=255)
                                   gggAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTcc                                                      <  1:1984336/93‑1 (MQ=255)
                                    ggAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGCTCTGGAAAGCCACGAGCGTTTTCCg                                                     >  2:1922960/1‑93 (MQ=255)
                                      agaTCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGa                                                   <  1:4557897/91‑1 (MQ=255)
                                          ccGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgc                                               >  1:1930841/1‑93 (MQ=255)
                                          ccGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgc                                               <  2:2978537/93‑1 (MQ=255)
                                           cGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgcg                                              <  2:2346331/93‑1 (MQ=255)
                                            gCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGGCACGAGCGTTTTCCGGAgcgcgc                                             >  2:3381924/1‑93 (MQ=255)
                                            gCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgcgc                                             >  1:2564391/1‑93 (MQ=255)
                                              aTTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCa                                           >  2:2094759/1‑93 (MQ=255)
                                                tGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAat                                         <  2:3470966/93‑1 (MQ=255)
                                                tGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAat                                         <  1:3079409/93‑1 (MQ=255)
                                                 gCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAata                                        >  1:1839098/1‑93 (MQ=255)
                                                 gCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAata                                        <  2:1591382/93‑1 (MQ=255)
                                                  cGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAatat                                       <  1:1285719/93‑1 (MQ=255)
                                                    tGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCg                                     >  1:1111993/1‑93 (MQ=255)
                                                     gATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCgg                                    <  2:1315911/93‑1 (MQ=255)
                                                     gATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCgg                                    <  1:1786398/93‑1 (MQ=255)
                                                      aTTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGt                                   <  2:919579/93‑1 (MQ=255)
                                                      aTTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGt                                   >  1:2782876/1‑93 (MQ=255)
                                                      aTTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTCGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGt                                   <  1:1588672/93‑1 (MQ=255)
                                                       ttCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGtt                                  <  1:433313/93‑1 (MQ=255)
                                                       ttCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGtt                                  >  1:1452478/1‑93 (MQ=255)
                                                        tCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATCTCGGttt                                 >  2:1058374/1‑93 (MQ=255)
                                                        tCAGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGttt                                 <  2:1758553/93‑1 (MQ=255)
                                                          aGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTa                               <  2:3833615/93‑1 (MQ=255)
                                                           ggTCGATGATGCCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGGTTTCCGGAGCGCGCCAATATCGGTTTTAt                              >  2:2749278/1‑93 (MQ=255)
                                                             tCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGc                            <  1:1879655/93‑1 (MQ=255)
                                                              cGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCa                           <  1:3198507/93‑1 (MQ=255)
                                                               gatgatGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCaa                          >  2:2447526/1‑93 (MQ=255)
                                                                atgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAg                         <  1:2631070/93‑1 (MQ=255)
                                                                  gatgTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTg                       <  1:4037361/93‑1 (MQ=255)
                                                                   atgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTgg                      <  1:3718604/93‑1 (MQ=255)
                                                                    tgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGt                     >  1:1903381/1‑93 (MQ=255)
                                                                      tCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTa                   >  2:4158471/1‑93 (MQ=255)
                                                                       cGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTaa                  >  1:3306045/1‑93 (MQ=255)
                                                                       cGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTaa                  <  2:854629/93‑1 (MQ=255)
                                                                       cGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTaa                  <  2:1690108/93‑1 (MQ=255)
                                                                        gATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAg                 <  1:544534/93‑1 (MQ=255)
                                                                         aTACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAgc                >  2:1818465/1‑93 (MQ=255)
                                                                         aTACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAgc                <  2:1635400/93‑1 (MQ=255)
                                                                         aTACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAgc                >  1:3488540/1‑93 (MQ=255)
                                                                              gcggcgGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAg           >  2:366650/1‑93 (MQ=255)
                                                                               cggcggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGc          >  2:3807602/1‑93 (MQ=255)
                                                                                  cggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCGCGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCata       <  1:1611893/93‑1 (MQ=255)
                                                                                  cggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCata       <  1:770391/93‑1 (MQ=255)
                                                                                    gTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATAtt     <  1:423881/93‑1 (MQ=255)
                                                                                    gTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATAtt     <  1:1258439/93‑1 (MQ=255)
                                                                                      aGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCg   >  2:115298/1‑93 (MQ=255)
                                                                                       gAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGt  <  1:3829508/93‑1 (MQ=255)
                                                                                       gAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGt  <  1:4546261/93‑1 (MQ=255)
                                                                                       gAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGt  <  1:760759/93‑1 (MQ=255)
                                                                                       gAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGt  >  2:1417305/1‑93 (MQ=255)
                                                                                       gAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGt  <  2:157001/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGT  >  NC_000913/3992804‑3992983

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