Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,746,162 G→T 100% 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,1620GT98.8% 314.6 / ‑3.8 95L23L (CTG→CTTyiaO2,3‑diketo‑L‑gulonate‑binding periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (54/38);  total (54/41)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTGGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGGAGAGAACC  >  NC_000913/3746067‑3746258
                                                                                               |                                                                                                
aaaTGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTtgcgg                                                                                              <  1:1418204/100‑1 (MQ=255)
 aaTGGCTTAATTAAAAGGAAAATATTATGAAATCACGCTCTCTACCATACGCACTATTAATTGCTAGCCCGGCCGCATTCAGCACATCTTCTATgtcgg                                                                                              >  2:501005/1‑99 (MQ=255)
     gATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACaa                                                                                         <  2:1158831/100‑1 (MQ=255)
     gATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACaa                                                                                         >  1:2341977/1‑100 (MQ=255)
     gATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTACTTGCCTGGCGGCATTCAGCACATCTTCTCTTGCGGCACaa                                                                                         >  2:777572/1‑100 (MQ=255)
         aattaaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCtt                                                                                     >  1:1967998/1‑100 (MQ=255)
         aattaaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCt                                                                                      <  2:2217214/99‑1 (MQ=255)
          attaaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCttt                                                                                    <  1:289566/100‑1 (MQ=255)
          attaaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCttt                                                                                    >  1:811559/1‑100 (MQ=255)
          attaaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCttt                                                                                    >  1:2095827/1‑100 (MQ=255)
           ttaaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTa                                                                                   <  1:1282635/100‑1 (MQ=255)
           ttaaAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTa                                                                                   >  2:134512/1‑100 (MQ=255)
             aaaaGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACg                                                                                 >  1:2091934/1‑100 (MQ=255)
             aaaaGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACg                                                                                 >  1:459762/1‑100 (MQ=255)
             aaaaGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACg                                                                                 <  1:442648/100‑1 (MQ=255)
             aaaaGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACg                                                                                 <  1:2209437/100‑1 (MQ=255)
             aaaaGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACg                                                                                 >  1:1415280/1‑100 (MQ=255)
                aGGAAAATATTATGCAATTACGCCCTCTAACCTACTCAGCATTCATTGCTGGCCTTGCTGCATTCAGCACCTCTTCTCCAGCTGCACAATCTTTAc                                                                                  <  2:510781/96‑1 (MQ=255)
                 ggAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGttcc                                                                             >  2:1552754/1‑98 (MQ=255)
                 ggAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTc                                                                             <  1:2053780/100‑1 (MQ=255)
                 ggAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTc                                                                             >  2:1328650/1‑100 (MQ=255)
                 ggAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTc                                                                             >  2:579605/1‑100 (MQ=255)
                  gAAAATATTATGAAATTACGCTCTGTAACCTAGGCATTATTCATTGATGGCCTGGCTGAATTCAGCACATCTTCTCTTGCGGCACAATCTTTACtttta                                                                             >  2:1888113/1‑98 (MQ=255)
                   aaaaTATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCgg                                                                           <  1:457780/100‑1 (MQ=255)
                   aaaaTATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCgg                                                                           <  1:625845/100‑1 (MQ=255)
                    aaaTATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGt                                                                          <  2:171693/100‑1 (MQ=255)
                    aaaTATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTGGCGGCACAATCTTTACGTTTCGGt                                                                          <  1:952524/100‑1 (MQ=255)
                     aaTACTATAAAATTACTCTCTGTAACCTACGCCTTAATCGCTGCTGGCCTGGTCGCATTCACCATATGTTTTCTGGTGGCATAATCTATACGTTTCGCtt                                                                         <  2:607763/100‑1 (MQ=255)
                      atatTATGAAATTACGCTCGGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACGATCTTTACGTTTCGGTTa                                                                        <  1:2462431/100‑1 (MQ=255)
                       tattatGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTCGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTAt                                                                       <  2:740790/100‑1 (MQ=255)
                         ttatGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGa                                                                     >  2:1245803/1‑100 (MQ=255)
                          tatGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGaa                                                                    <  2:1687274/100‑1 (MQ=255)
                          tatGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGaa                                                                    <  2:2291592/100‑1 (MQ=255)
                          tatGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGaa                                                                    >  2:2035745/1‑100 (MQ=255)
                           atGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGaaa                                                                   >  1:850051/1‑100 (MQ=255)
                           atGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGaaa                                                                   >  2:1349892/1‑100 (MQ=255)
                            tGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAAc                                                                  <  1:529081/100‑1 (MQ=255)
                             gAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACa                                                                 <  2:1415269/100‑1 (MQ=255)
                             gAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACa                                                                 <  2:951689/100‑1 (MQ=255)
                                 ttACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATcac                                                             >  2:2457310/1‑100 (MQ=255)
                                       tctGTAACCTACGCATTATTCATTGCTGGCCTGCCTGCATTCAGCACATCTTCTCTTGCGGCACACTCTTTACGGTTCGGTTATGAAACATCACACAcc                                                        >  2:1227303/1‑99 (MQ=255)
                                         tGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGAc                                                     <  1:1349901/100‑1 (MQ=255)
                                           tAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTc                                                   >  1:1085056/1‑100 (MQ=255)
                                              ccTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCa                                                <  1:2121382/100‑1 (MQ=255)
                                              ccTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCa                                                <  1:289126/100‑1 (MQ=255)
                                              ccTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCa                                                >  1:816666/1‑100 (MQ=255)
                                                tACTCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAAc                                              <  1:1754060/100‑1 (MQ=255)
                                                tACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAAc                                              <  1:481145/100‑1 (MQ=255)
                                                 aCGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACa                                             >  1:2006204/1‑100 (MQ=255)
                                                  cGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACat                                            >  1:780568/1‑100 (MQ=255)
                                                   gCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACata                                           >  2:673054/1‑100 (MQ=255)
                                                       tattCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgc                                       >  1:1315584/1‑100 (MQ=255)
                                                        attCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcg                                      <  1:823227/100‑1 (MQ=255)
                                                        aatCATTGCTGGCATGGCTGCATTCAGCACAACTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcg                                      >  1:1728634/3‑100 (MQ=255)
                                                         ttCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcgg                                     <  1:2128758/100‑1 (MQ=255)
                                                          tCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcggc                                    >  2:1819107/1‑100 (MQ=255)
                                                          tCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcggc                                    <  2:462596/100‑1 (MQ=255)
                                                           cATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcggcg                                   >  1:1279487/1‑100 (MQ=255)
                                                           cATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTgcggcg                                   >  1:2328299/1‑100 (MQ=255)
                                                             ttGCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGaa                                 >  1:1610353/1‑100 (MQ=255)
                                                               gCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGaaaa                               >  2:1916520/1‑100 (MQ=255)
                                                               gCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGaaaa                               >  2:386740/1‑100 (MQ=255)
                                                               gCTGGCCTGGCTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGaaaa                               >  1:2035863/1‑100 (MQ=255)
                                                                 tggcctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGaaaaaa                             >  1:1237606/1‑100 (MQ=255)
                                                                   gcctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAAtt                           <  1:803485/100‑1 (MQ=255)
                                                                   gcctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAAtt                           <  1:2163933/100‑1 (MQ=255)
                                                                   gcctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAAtt                           <  1:51491/100‑1 (MQ=255)
                                                                    cctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTc                          <  2:1315575/100‑1 (MQ=255)
                                                                    cctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTc                          >  1:1998335/1‑100 (MQ=255)
                                                                    cctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTc                          <  2:758875/100‑1 (MQ=255)
                                                                    cctggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTCCGGCGaaaaa                              >  2:749156/1‑96 (MQ=255)
                                                                       ggcTGCATTCAGCACCTCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAAt                       <  1:2515687/100‑1 (MQ=255)
                                                                       ggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAAt                       >  2:682290/1‑100 (MQ=255)
                                                                       ggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAAt                       >  1:414380/1‑100 (MQ=255)
                                                                       ggcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAAt                       >  2:1358653/1‑100 (MQ=255)
                                                                        gcTGTATTCATCACATCTTCTCTTGCGGCACAATCTTTACGTTTTGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAAt                       <  2:2506698/99‑1 (MQ=255)
                                                                        gcTGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATg                      >  1:1990021/1‑100 (MQ=255)
                                                                          tGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGAt                    >  2:1920148/1‑100 (MQ=255)
                                                                          tGCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGAt                    >  1:2500755/1‑100 (MQ=255)
                                                                           gCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGAtt                   >  2:627227/1‑100 (MQ=255)
                                                                           gCATTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGAtt                   >  2:1303019/1‑100 (MQ=255)
                                                                             aTTCAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTa                 >  2:208628/1‑100 (MQ=255)
                                                                                cAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTg              >  1:2545768/1‑100 (MQ=255)
                                                                                cAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTg              >  2:1514959/1‑100 (MQ=255)
                                                                                cAGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTg              >  1:785278/1‑100 (MQ=255)
                                                                                 aGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCACAACATATTGCGGCGAAAAAATTCAATGATTTATTGc             <  1:9859/100‑1 (MQ=255)
                                                                                 aGCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTAGGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGc             <  2:517687/100‑1 (MQ=255)
                                                                                  gCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCa            <  2:183778/100‑1 (MQ=255)
                                                                                  gCACATCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCa            >  2:1355804/1‑100 (MQ=255)
                                                                                   cacaTCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAg           <  1:1096568/100‑1 (MQ=255)
                                                                                   cacaTATTCTCTTGCGGCACACTCTTGACGTTTCGGTTATGAAACTTCACAAACCGACTCGCCACACATTGAGACGAGAAAAATCagagaatta                 >  2:2155910/1‑86 (MQ=255)
                                                                                    acaTCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAgg          <  2:1001971/100‑1 (MQ=255)
                                                                                     caTCTTCTCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGga         >  2:1893624/1‑100 (MQ=255)
                                                                                       tcttctCTTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGgaga       >  1:1662239/1‑100 (MQ=255)
                                                                                         ttctCATGCGGCACAATTTTTACGTTTCGGTTATGAAACATCACAAACCGACTCTCAACATATAGCGGCGAAAAAATTCAATGATTTATTGCAGgagag      <  2:2472288/99‑1 (MQ=255)
                                                                                            tctTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGTGGCGAAAAAATTCAATGATTTATTGCAGGAGAGAAcc  >  1:2364795/1‑100 (MQ=255)
                                                                                            tctTGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGGAGAGAAc   >  2:654377/1‑99 (MQ=255)
                                                                                               |                                                                                                
AAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTGGCGGCACAATCTTTACGTTTCGGTTATGAAACATCACAAACCGACTCGCAACATATTGCGGCGAAAAAATTCAATGATTTATTGCAGGAGAGAACC  >  NC_000913/3746067‑3746258

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