Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,746,162 G→T L23L (CTG→CTT yiaO → 2,3‑diketo‑L‑gulonate‑binding periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,746,1620GT99.3% 502.8 / ‑4.6 149L23L (CTG→CTTyiaO2,3‑diketo‑L‑gulonate‑binding periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (73/75);  total (74/75)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 156 total aligned reads displayed.

GAAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTGGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGGAGAGAACC  >  NC_000913/3746066‑3746258
                                                                                                |                                                                                                
gAAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTtgc                                                                                                >  1:3190950/1‑99 (MQ=255)
 aaaTGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGAAATATGAATAGCTGACCTGGCTGCATTCAGCACATATTCTCTtgcg                                                                                               >  1:3018480/1‑99 (MQ=255)
  aaTGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTtgcggc                                                                                             >  2:930342/1‑100 (MQ=255)
  aaTGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTtgcggc                                                                                             <  2:1043690/100‑1 (MQ=255)
     ggATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGcaca                                                                                          >  2:2488696/1‑100 (MQ=255)
        ttaattaaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATc                                                                                       <  2:1227865/100‑1 (MQ=255)
          aattaaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCtt                                                                                     <  1:2796998/100‑1 (MQ=255)
           attaaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTTAGCACATCTTCTCTTGCGGCACAATCttt                                                                                    >  1:1666860/1‑100 (MQ=255)
           attaaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTTAGAACATCTTCTCTTGCGGCACAATCtt                                                                                     >  1:273919/1‑99 (MQ=255)
             taaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCATCACATCTTCTCTTGCGGCACAATCTTTAc                                                                                  >  2:692557/1‑100 (MQ=255)
             taaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTAc                                                                                  <  2:419742/100‑1 (MQ=255)
              aaaaGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAGTCTTTACg                                                                                 >  1:804872/1‑100 (MQ=255)
              aaaaGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACg                                                                                 >  2:3150146/1‑100 (MQ=255)
                 aGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGttt                                                                              >  1:1716082/1‑100 (MQ=255)
                  ggAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTc                                                                             >  2:1123163/1‑100 (MQ=255)
                  ggAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTc                                                                             >  2:1730324/1‑100 (MQ=255)
                  ggAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTc                                                                             >  1:307502/1‑100 (MQ=255)
                  ggAAAATATTATGAAATTACGCTCTGTAACCTACGCAATATTCATTGATGGCCTGGCTGCATTCAGCACATATTCTCTTGCGGCCCAATCTTTACGtt                                                                               >  1:3143591/1‑98 (MQ=255)
                    aaaaTATTATGAAATTACGCTCTGTATCCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCgg                                                                           <  1:1510904/100‑1 (MQ=255)
                    aaaaTATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCgg                                                                           <  2:923279/100‑1 (MQ=255)
                    aaaaTATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCgg                                                                           <  2:3014440/100‑1 (MQ=255)
                     aaaTATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGt                                                                          <  2:1716105/100‑1 (MQ=255)
                       atatTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTa                                                                        <  2:75202/100‑1 (MQ=255)
                       atatTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTa                                                                        >  1:1085898/1‑100 (MQ=255)
                       atatTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTa                                                                        <  2:2236809/100‑1 (MQ=255)
                       atatTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTa                                                                        >  2:885068/1‑100 (MQ=255)
                       atatTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTa                                                                        <  1:28141/100‑1 (MQ=255)
                        tattatGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTAt                                                                       <  2:2525604/100‑1 (MQ=255)
                         attatGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATg                                                                      <  2:2119371/100‑1 (MQ=255)
                           tatGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGaa                                                                    >  2:2670497/1‑100 (MQ=255)
                           tatGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGaa                                                                    >  1:2139018/1‑100 (MQ=255)
                            atGAAATTACGCTCTGTAACCTACGCATTAGTCCTTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTGGGTTATGaaa                                                                   >  1:2285673/1‑100 (MQ=255)
                             tgcacttACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTAACCACATCTTCTCTTGCGGCACAATATTTACGTTTCGGTTATGAAAc                                                                  <  2:3305444/95‑1 (MQ=255)
                             tGAAATTACGCTCTGTAACCTACGCATTATTCATTGTTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAAc                                                                  >  1:1502227/1‑100 (MQ=255)
                             tGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAAc                                                                  >  2:2801853/1‑100 (MQ=255)
                               aaaTTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACAt                                                                <  2:1992518/100‑1 (MQ=255)
                               aaaTTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACAt                                                                >  2:748051/1‑100 (MQ=255)
                               aaaTTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACAt                                                                >  1:1455029/1‑100 (MQ=255)
                               aaaTTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGCCACAATCTTTACGTTTCGGTTATGAAAc                                                                  >  1:569315/1‑98 (MQ=255)
                                aaTTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATc                                                               <  1:2646997/100‑1 (MQ=255)
                                     cGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACaaa                                                          <  2:1499226/100‑1 (MQ=255)
                                     cGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACaaa                                                          <  1:536549/100‑1 (MQ=255)
                                       ctctGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAAcc                                                        >  1:2773735/1‑100 (MQ=255)
                                            tAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTc                                                   >  2:1982958/1‑100 (MQ=255)
                                            tAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTc                                                   >  2:3420102/1‑100 (MQ=255)
                                             aaCCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCg                                                  >  1:1458039/1‑100 (MQ=255)
                                              aCCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGc                                                 >  2:2199797/1‑100 (MQ=255)
                                               ccTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCa                                                <  1:744654/100‑1 (MQ=255)
                                               ccTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCa                                                >  2:1954617/1‑100 (MQ=255)
                                               ccTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCa                                                <  2:1753151/100‑1 (MQ=255)
                                               ccTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCa                                                >  2:1693913/1‑100 (MQ=255)
                                                 tACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAAc                                              >  1:330753/1‑100 (MQ=255)
                                                 tACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAAc                                              <  2:275708/100‑1 (MQ=255)
                                                  aCGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACa                                             >  2:2637444/1‑100 (MQ=255)
                                                    gCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACata                                           <  1:2985688/100‑1 (MQ=255)
                                                    gCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACata                                           >  2:3079838/1‑100 (MQ=255)
                                                    gCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACata                                           <  2:50648/100‑1 (MQ=255)
                                                      attattCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATAtt                                         >  1:807709/1‑100 (MQ=255)
                                                       ttattCAGTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTg                                        <  1:2331100/100‑1 (MQ=255)
                                                          ttCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcgg                                     <  1:2372751/100‑1 (MQ=255)
                                                          ttCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcgg                                     <  2:275409/100‑1 (MQ=255)
                                                           tCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcggc                                    <  2:133734/100‑1 (MQ=255)
                                                            cATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcggcg                                   >  2:2792606/1‑100 (MQ=255)
                                                            cATTGCTGGCCTGGCTACATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcggcg                                   <  1:440346/100‑1 (MQ=255)
                                                             aTTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGa                                  >  2:1602152/1‑100 (MQ=255)
                                                             aTTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGa                                  <  2:364049/100‑1 (MQ=255)
                                                                gCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGaaaa                               >  2:1289929/1‑100 (MQ=255)
                                                                  tggcctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGaaaaaa                             <  2:1930914/100‑1 (MQ=255)
                                                                     cctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTc                          >  1:2970732/1‑100 (MQ=255)
                                                                     cctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTc                          <  2:619480/100‑1 (MQ=255)
                                                                     cctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTc                          <  1:609075/100‑1 (MQ=255)
                                                                     cctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTc                          >  1:361665/1‑100 (MQ=255)
                                                                      ctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCa                         <  2:1593393/100‑1 (MQ=255)
                                                                       tggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCAAAAAAAAAATTCaa                        <  1:747325/100‑1 (MQ=255)
                                                                       tggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGATTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCac                        <  1:1557427/100‑2 (MQ=255)
                                                                        ggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAAt                       <  2:28382/100‑1 (MQ=255)
                                                                        ggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAAt                       <  2:1776692/100‑1 (MQ=255)
                                                                          cTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGa                     >  2:2124108/1‑100 (MQ=255)
                                                                          cTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGa                     >  2:3295129/1‑100 (MQ=255)
                                                                               ttCGGCACATCTTCTGTTGCGGCACATTCTTTACGTTTCGGTTATGAATCTTCACAAACCGACTCGCAACATATTGTGGCGAAAAAATTCAATGATTTAt                >  1:2251868/1‑100 (MQ=255)
                                                                               ttCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTAt                >  1:1498768/1‑100 (MQ=255)
                                                                                tCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATTAAACATCACAAACCGACTCGCAACATATTTCGGCGAAAAAATTCAATGATTTAtt               >  1:1770896/1‑100 (MQ=255)
                                                                                 cAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTg              <  1:113275/100‑1 (MQ=255)
                                                                                 cAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTg              >  2:1224777/1‑100 (MQ=255)
                                                                                 cAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTg              >  2:359322/1‑100 (MQ=255)
                                                                                  aGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGc             >  2:2416179/1‑100 (MQ=255)
                                                                                  aGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGc             <  2:674642/100‑1 (MQ=255)
                                                                                  aGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGc             >  2:1418127/1‑100 (MQ=255)
                                                                                  aGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGc             <  2:3065191/100‑1 (MQ=255)
                                                                                     acaTCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCAAAAAAATTCAATGATTTATTGCAgg          >  1:2174517/1‑100 (MQ=255)
                                                                                      caTCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGga         >  1:594728/1‑100 (MQ=255)
                                                                                      caTCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGga         <  2:1085912/100‑1 (MQ=255)
                                                                                      caTCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGga         >  1:3236193/1‑100 (MQ=255)
                                                                                         cttctCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGgagag      <  1:571697/100‑1 (MQ=255)
                                                                                         cttctCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGgagag      <  2:3227073/100‑1 (MQ=255)
                                                                                           tctcTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGGAGAGaa    <  2:3408421/100‑1 (MQ=255)
                                                                                            ctctTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGGAGAGAAc   <  2:2187630/100‑1 (MQ=255)
                                                                                            ctctTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGGAGAGAAc   <  2:2392948/100‑1 (MQ=255)
                                                                                             tctTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGGAGAGAAcc  >  2:1664029/1‑100 (MQ=255)
                                                                                             tctTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGGAGAGAAcc  >  1:395666/1‑100 (MQ=255)
                                                                                                |                                                                                                
GAAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTGGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGGAGAGAACC  >  NC_000913/3746066‑3746258

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