Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,140,718 T→A 15.6% I502F (ATC→TTC)  frwA ← putative PTS enzyme: Hpr, enzyme I and II components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,140,7180TA15.6% 385.3 / 58.7 155I502F (ATC→TTC) frwAputative PTS enzyme: Hpr, enzyme I and II components
Reads supporting (aligned to +/- strand):  ref base T (66/64);  new base A (11/13);  total (78/77)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-01

Only 100 of 157 total aligned reads displayed.

CTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATA  >  NC_000913/4140620‑4140817
                                                                                                  |                                                                                                   
ctcATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATc                                                                                                    <  2:1603886/100‑1 (MQ=255)
 tcatcaAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCa                                                                                                   <  2:828762/100‑1 (MQ=255)
 tcatcaAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCa                                                                                                   <  1:163731/100‑1 (MQ=255)
  catcaAACGGAATGTGTTCGTTACTTAGTTTCTGTTTTGCTTCCTCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAAGAGATCATCGGGAACa                                                                                                   <  2:1999302/99‑1 (MQ=255)
  catcaAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCAt                                                                                                  >  1:1139721/1‑100 (MQ=255)
  catcaAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAACAt                                                                                                  >  2:993632/1‑100 (MQ=255)
   atcaAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAtaa                                                                                                   >  1:215022/1‑96 (MQ=255)
   atcaAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATg                                                                                                 <  1:2404622/100‑1 (MQ=255)
   atcaAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATg                                                                                                 <  1:1385342/100‑1 (MQ=255)
    tcaAACGGTATGTGTTCGTTACGTAGTTTCTGTTTGGCTTCCGCCAGTTTTTCTTTCCCCCATAATATCTCTTCCATTGATGAGATCATCGGGATCATGa                                                                                                <  2:571825/100‑1 (MQ=255)
      aaaCGGAATGTGTTCGTAACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTACTTTCACCCATAAGATCTCTTCCATTGAGGAGCTCATCGGGATCATgct                                                                                               >  2:446164/1‑97 (MQ=255)
        aCGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGAtttt                                                                                            <  1:776322/100‑1 (MQ=255)
         cGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTc                                                                                           <  2:1759396/100‑1 (MQ=255)
            actgtgTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAgg                                                                                        <  1:2282454/98‑1 (MQ=255)
             aTGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAACATGATTTTCAGGc                                                                                       <  1:1801895/100‑1 (MQ=255)
              tgtgTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCt                                                                                      >  1:1982130/1‑100 (MQ=255)
               gtgtTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTg                                                                                     >  2:2398117/1‑100 (MQ=255)
                   tCGTTACGTAGTTTCTGGTTGGCTGCCGCCAGCTTTGCTTTAACCCATAAGATCTCTTCCATTGAGGAGTTCATCTTGATCATGATTTTCAGGCTGCCg                                                                                  <  2:745626/99‑1 (MQ=255)
                    cGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTg                                                                                >  1:2161063/1‑100 (MQ=255)
                      ttACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTg                                                                                <  1:940396/98‑1 (MQ=255)
                        aCGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTTTTTTTTTTTTTAGCAGTATTCGATTGAGGAGATCATCGAGATCCTGATTTTCAGGCTGc                                                                                    >  2:495260/1‑92 (MQ=255)
                        aCGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGc                                                                             >  1:622835/1‑99 (MQ=255)
                        aCGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCg                                                                            <  1:1115638/100‑1 (MQ=255)
                         cGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATTATCGGTAACATGATTTTCAGGCTGCCGTGAGCg                                                                            >  2:760677/1‑99 (MQ=255)
                           tAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAACATGATTTTCAGGCTGCCGTGAGCGGAg                                                                         >  1:1018637/1‑100 (MQ=255)
                           tAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAACATGATTTTCAGGCTGCCGTGAGCGGAg                                                                         <  1:1271527/100‑1 (MQ=255)
                             gTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGgc                                                                       <  2:672566/100‑1 (MQ=255)
                               tGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGgcgc                                                                     >  2:1190228/1‑100 (MQ=255)
                                 cTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCgg                                                                   >  1:1560814/1‑100 (MQ=255)
                                 cTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAACATGATTTTCAGGCTGCCGTGAGCGGATGCGCgg                                                                   >  1:1861950/1‑100 (MQ=255)
                                 cTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAACATGATTTTCAGGCTGCCGTGAGCGGAGGCGCgg                                                                   >  2:752855/1‑100 (MQ=255)
                                  tGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCgga                                                                  <  2:1147864/100‑1 (MQ=255)
                                  tGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCgga                                                                  <  1:747592/100‑1 (MQ=255)
                                  tGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAACATGATTTTCAGGCTGCCGTGAGCGGAGGCGCgga                                                                  <  1:1017278/100‑1 (MQ=255)
                                   gTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCggag                                                                 <  1:394464/100‑1 (MQ=255)
                                   gTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCggag                                                                 <  1:718518/100‑1 (MQ=255)
                                    tttGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCggag                                                                 <  1:724298/99‑1 (MQ=255)
                                      ttgCTTCCGCCAGTTTTTCTTTAACCCATAAGATCTCTTCCATTCAGGAGATCATCGGGATCATGATTTTCAGGCTGCGGTGCGGGGAGAAGCGGAGGAt                                                              >  2:1884690/3‑100 (MQ=255)
                                       ggCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATc                                                             >  1:2188846/1‑100 (MQ=255)
                                       ggCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATc                                                             <  2:1937840/100‑1 (MQ=255)
                                       ggCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAACATGATTTTCATGCTTCCTTGAGCGGAGGCGCGGAGGATc                                                             <  2:241446/100‑1 (MQ=255)
                                         cTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCAATTGAGGAGATCATCGGGATCATGATTTTCAGTCTGTCGTGAGCGGAGGCGCGGAGGATCGa                                                           >  2:39573/1‑100 (MQ=255)
                                         cTTCCGACAGTTTTTCATTAAACAATAAAATCTATTCCATTGAGTAACTCATCGGGCTCATGATTTTCATACTGCAGTGGGGGGAGGCGCGGAGGATCg                                                            >  2:791948/1‑99 (MQ=255)
                                          ttCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGAc                                                          <  2:980135/100‑1 (MQ=255)
                                          ttCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGAc                                                          <  1:873502/100‑1 (MQ=255)
                                           tCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGAcc                                                         <  1:2080680/100‑1 (MQ=255)
                                           tCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGAcc                                                         <  1:555215/100‑1 (MQ=255)
                                           tCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGAcc                                                         <  2:1991320/100‑1 (MQ=255)
                                            ccGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCg                                                        <  1:2077704/100‑1 (MQ=255)
                                             cGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGt                                                       <  2:2364986/100‑1 (MQ=255)
                                             cGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGt                                                       >  1:1175070/1‑100 (MQ=255)
                                               ccAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAg                                                     <  2:927872/100‑1 (MQ=255)
                                                 aGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGTGGAGGATCGACCGTAGc                                                    <  2:323872/99‑1 (MQ=255)
                                                  gTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTg                                                  <  2:551945/100‑1 (MQ=255)
                                                  gTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTg                                                  >  2:1528168/1‑100 (MQ=255)
                                                     tttctttcACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGt                                               <  2:2232873/100‑1 (MQ=255)
                                                       tctttcACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTgg                                             >  1:754676/1‑100 (MQ=255)
                                                       tctttcACCCATAAGATCTCTTCCAGTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTgg                                             <  2:66708/100‑1 (MQ=255)
                                                        ctttcACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGt                                            >  2:1698813/1‑100 (MQ=255)
                                                        ctttcACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGt                                            >  2:1313216/1‑100 (MQ=255)
                                                         tttcACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTa                                           >  1:2402871/1‑100 (MQ=255)
                                                         tttcACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTa                                           >  2:476729/1‑100 (MQ=255)
                                                         tttcACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTCCGTGGTa                                           >  1:681476/1‑100 (MQ=255)
                                                         tttcACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAACATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTa                                           <  1:1473852/100‑1 (MQ=255)
                                                          ttcACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTaa                                          <  2:335296/100‑1 (MQ=255)
                                                            cACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAaac                                        >  1:490299/1‑100 (MQ=255)
                                                             aCCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAACATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAaaca                                       <  2:221747/100‑1 (MQ=255)
                                                              cccATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCCGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAaacaa                                      >  2:463579/1‑100 (MQ=255)
                                                               ccATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAaacaac                                     <  1:476732/100‑1 (MQ=255)
                                                                cATAAGATCTCTTCCATTAAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGAGAAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAacg                                    >  2:399156/1‑100 (MQ=255)
                                                                    aGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGAcgc                                <  1:1209596/100‑1 (MQ=255)
                                                                    aGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGAcgc                                >  2:1007184/1‑100 (MQ=255)
                                                                      aTCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGt                              <  1:1854269/100‑1 (MQ=255)
                                                                        ctctTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTAc                            >  1:963888/1‑100 (MQ=255)
                                                                        ctctTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTAc                            <  1:2143678/100‑1 (MQ=255)
                                                                        ctctTCCATTGAGGAAAACATCGGGATCATGATTTTCAGGATGACGTGAACGGAGGCGCGGAGGATCGACCGTAGCTGCGTAGTAAACAACGACGCGTa                             >  2:453932/1‑99 (MQ=255)
                                                                         tctTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTAct                           >  2:1518643/1‑100 (MQ=255)
                                                                         tctTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTAc                            >  2:1064951/1‑99 (MQ=255)
                                                                          ctTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTActc                          >  2:1832088/1‑100 (MQ=255)
                                                                            tCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCtt                        >  1:774994/1‑100 (MQ=255)
                                                                             ccATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACGGTAGCTGCGTGGTAAACAACGACGCGTACTCTTc                       >  2:1255577/1‑100 (MQ=255)
                                                                             ccATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTc                       >  2:976607/1‑100 (MQ=255)
                                                                              cATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTc                       <  2:47246/99‑1 (MQ=255)
                                                                              cATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCa                      <  1:1839468/100‑1 (MQ=255)
                                                                              cATTGAGGAGATCAGCGGGATCATGATTTTCAGTCTGCCGTGAGCGGAGTCGCGGAGGATCGACCGTAGCTTCGTGGTAAACACCGACGCGTACTCTTCa                      >  2:1395242/1‑100 (MQ=255)
                                                                               aTTGAGGAGATCATCGGGAACATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCAt                     <  2:1860274/100‑1 (MQ=255)
                                                                                tttaggagATCATCGGGATCATGAGTTTCAGGCTGCCGTGGGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAATCAACGACGCGTACTCTTCAt                     <  2:2257432/96‑1 (MQ=255)
                                                                                 tGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATaa                   <  1:592155/100‑1 (MQ=255)
                                                                                  gaggagATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATaaa                  >  1:2335867/1‑100 (MQ=255)
                                                                                   aggagATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAAt                 <  1:660053/100‑1 (MQ=255)
                                                                                   aggagATCATCGGGATCATGATGTTTAGGCTGTCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTTGTAATTTTTTACGCGTACTCATCATa                    <  2:887873/97‑1 (MQ=255)
                                                                                    ggagATCATCGGGAACATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATa                >  2:2234514/1‑100 (MQ=255)
                                                                                     gagaTCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATAc               <  1:2046029/100‑1 (MQ=255)
                                                                                      agaTCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACg              >  2:1888332/1‑100 (MQ=255)
                                                                                        atcatcGGGAGCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCa            >  1:934287/1‑100 (MQ=255)
                                                                                         tcatcGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCAc           >  1:2175971/1‑100 (MQ=255)
                                                                                             cGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGgcg       >  1:1802010/1‑100 (MQ=255)
                                                                                             cGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGgcg       >  1:958442/1‑100 (MQ=255)
                                                                                               ggATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGGGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGcgcg     >  2:590235/1‑100 (MQ=255)
                                                                                                 aTCATGATTTTCAGGCTGCCGTGAGCGGAGTCGCGGAGGATCGACTGTAGCTGCGTGGTAAACAACGACGTGTACTCTTCATAAATACGCACGGCGCGa    <  2:304908/99‑1 (MQ=255)
                                                                                                  |                                                                                                   
CTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATA  >  NC_000913/4140620‑4140817

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