Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,118,221 G→A Q62* (CAA→TAA)  wcaK ← colanic acid biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,118,2210GA97.9% 334.5 / ‑2.1 96Q62* (CAA→TAA) wcaKcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (44/50);  total (44/52)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTGGATTGAG  >  NC_000913/2118124‑2118317
                                                                                                 |                                                                                                
gTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACa                                                                                                <  1:2039274/100‑1 (MQ=255)
 tGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTTTTGTGTTGTTTCATTTACAg                                                                                               <  2:1831443/100‑1 (MQ=255)
 tGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAg                                                                                               >  1:520844/1‑100 (MQ=255)
  gACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAgg                                                                                              <  2:2659040/100‑1 (MQ=255)
   aCGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGa                                                                                             <  1:1510822/100‑1 (MQ=255)
      cgTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAAcn                                                                                          <  1:2934642/100‑2 (MQ=255)
      cgTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACa                                                                                          >  1:1665048/1‑100 (MQ=255)
       gTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAg                                                                                         >  2:1556166/1‑100 (MQ=255)
        tGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGc                                                                                        <  1:2317351/100‑1 (MQ=255)
         gagaGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCg                                                                                       >  2:2524904/1‑100 (MQ=255)
             gCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATc                                                                                   >  1:2643032/1‑100 (MQ=255)
                atgccTGATGCTGGGAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGcc                                                                                <  2:1273462/97‑1 (MQ=255)
                aTACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGcc                                                                                <  1:1391948/100‑1 (MQ=255)
                aTACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGcc                                                                                <  2:224995/100‑1 (MQ=255)
                aTACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGcc                                                                                <  2:2817266/100‑1 (MQ=255)
                aTACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCCGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGcc                                                                                <  1:2110327/100‑1 (MQ=255)
                 tACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGccc                                                                               >  1:1289466/1‑100 (MQ=255)
                  aCCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCa                                                                              >  2:2097814/1‑100 (MQ=255)
                   ccTGATGGTGGTAGCGGCGACGGAGGACTTTTTTAACGCCCCCGACACCGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGcctt                                                                              <  2:762023/99‑3 (MQ=255)
                   ccTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGGGCTGTTGTGATGTTTCATTTACAGGAACAGCGGATCGCCCAt                                                                             >  2:2094612/1‑100 (MQ=255)
                   ccTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCAt                                                                             >  1:2052108/1‑100 (MQ=255)
                    cTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCAtt                                                                            <  1:274539/100‑1 (MQ=255)
                    cTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCAtt                                                                            >  2:349551/1‑100 (MQ=255)
                     tGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTa                                                                           <  1:2741033/100‑1 (MQ=255)
                      gATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTAc                                                                          <  2:2937578/100‑1 (MQ=255)
                      gATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTAc                                                                          >  2:1598702/1‑100 (MQ=255)
                       aTGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTAcc                                                                         <  2:2218563/100‑1 (MQ=255)
                       aTGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTAcc                                                                         <  1:1059041/100‑1 (MQ=255)
                             gTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCgg                                                                   <  2:34676/100‑1 (MQ=255)
                               aGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGtt                                                                 <  2:1650528/100‑1 (MQ=255)
                               aGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGtt                                                                 <  2:1289460/100‑1 (MQ=255)
                                 cggcgACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGa                                                               >  1:1084148/1‑100 (MQ=255)
                                  ggcgACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                                              <  1:2543387/100‑1 (MQ=255)
                                  ggcgACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                                              <  2:761451/100‑1 (MQ=255)
                                  ggcgACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                                              <  1:51463/100‑1 (MQ=255)
                                  ggcgACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                                              <  1:2961181/100‑1 (MQ=255)
                                   gcgACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGagc                                                             <  1:1647664/100‑1 (MQ=255)
                                       cGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGcc                                                         >  1:734885/1‑100 (MQ=255)
                                       cGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGcc                                                         <  1:446906/100‑1 (MQ=255)
                                          aggaCTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAgg                                                      <  2:1254026/100‑1 (MQ=255)
                                          aggaCTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAgg                                                      <  2:1621824/100‑1 (MQ=255)
                                           ggaCTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGCACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGga                                                     >  1:607437/1‑100 (MQ=255)
                                           ggaCTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGga                                                     <  2:2643017/100‑1 (MQ=255)
                                            gaCTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGgaa                                                    >  1:855034/1‑100 (MQ=255)
                                              cTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGgaaga                                                  >  1:1527483/1‑100 (MQ=255)
                                              cTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGgaag                                                   >  2:581147/1‑99 (MQ=255)
                                                tttttAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAAc                                                >  1:1653537/1‑100 (MQ=255)
                                                tttttAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAAc                                                <  2:1694711/100‑1 (MQ=255)
                                                tttttAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAAc                                                >  2:165070/1‑100 (MQ=255)
                                                tttttAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAAc                                                <  2:1120556/100‑1 (MQ=255)
                                                tttttAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAAc                                                <  2:1120577/100‑1 (MQ=255)
                                                 ttttAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACt                                               <  1:1212018/100‑1 (MQ=255)
                                                 ttttAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACt                                               >  2:1385830/1‑100 (MQ=255)
                                                    tAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGAc                                            <  2:1001935/100‑1 (MQ=255)
                                                      aCGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCg                                          <  1:1186387/100‑1 (MQ=255)
                                                      aCGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCg                                          >  1:2979437/1‑100 (MQ=255)
                                                          gcCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATa                                      <  1:2524919/100‑1 (MQ=255)
                                                            ccGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATa                                      <  1:597419/98‑1 (MQ=255)
                                                            ccGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGc                                    <  2:2662230/100‑1 (MQ=255)
                                                            ccGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGc                                    <  2:2883431/100‑1 (MQ=255)
                                                             cGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCg                                   <  2:1380749/100‑1 (MQ=255)
                                                               aCAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGTCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGc                                 <  2:1607803/100‑1 (MQ=255)
                                                               aCAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGc                                 <  2:2302437/100‑1 (MQ=255)
                                                                  aCGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtca                              >  1:2568527/1‑100 (MQ=255)
                                                                  aCGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtca                              >  2:139099/1‑100 (MQ=255)
                                                                  aCGCCCGCCGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGTGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtcat                             >  2:440416/1‑100 (MQ=255)
                                                                    gcccgccGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtcatc                            <  1:1385838/100‑1 (MQ=255)
                                                                    gcccgccGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtcatc                            <  2:179116/100‑1 (MQ=255)
                                                                     cccgccGCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtcatca                           >  1:1466773/1‑100 (MQ=255)
                                                                      ccgccgCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATcac                          <  1:2268196/100‑1 (MQ=255)
                                                                      ccgccgCGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGAAAGCGGCTCATcac                          >  1:1537590/1‑100 (MQ=255)
                                                                       cgccgcGCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATcaca                         >  1:2424672/1‑100 (MQ=255)
                                                                          cgcgCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATcc                      <  1:349553/100‑1 (MQ=255)
                                                                          cgcgCTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATcc                      <  2:379639/100‑1 (MQ=255)
                                                                             gcTGTTGTGTTGTTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACt                   <  2:1215010/100‑1 (MQ=255)
                                                                                  tgtgttgtTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGc              <  2:115888/100‑1 (MQ=255)
                                                                                  tgtgttgtTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGc              >  1:809119/1‑100 (MQ=255)
                                                                                   gtgttgtTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCg             >  1:1060817/1‑100 (MQ=255)
                                                                                   gtgttgtTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCg             >  1:484250/1‑100 (MQ=255)
                                                                                    tgttgtTTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCgt            >  1:2331050/1‑100 (MQ=255)
                                                                                     gttgttTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCgtg           <  2:1528743/100‑1 (MQ=255)
                                                                                     gttgttTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCgtg           >  2:2476257/1‑100 (MQ=255)
                                                                                       tgttTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGtgtg         >  1:1622433/1‑100 (MQ=255)
                                                                                        gttTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTgg        >  1:243289/1‑100 (MQ=255)
                                                                                        gttTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTgg        >  2:794767/1‑100 (MQ=255)
                                                                                        gttTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTgg        >  1:1546160/1‑100 (MQ=255)
                                                                                        gttTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTgg        >  2:561336/1‑100 (MQ=255)
                                                                                        gttTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTgg        >  1:686487/1‑100 (MQ=255)
                                                                                        gttTCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTgg        >  1:987522/1‑100 (MQ=255)
                                                                                        gttTAATTTACAGGAACAGCGGATCGCCCATTACCGGGCTGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTgg        >  2:2790864/1‑100 (MQ=255)
                                                                                           tCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTGGAtt     >  1:2901556/1‑100 (MQ=255)
                                                                                           tCATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTGGAtt     <  1:1312158/100‑1 (MQ=255)
                                                                                            cATTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTGGATTg    <  1:165072/100‑1 (MQ=255)
                                                                                             aTTTACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTGGATTga   >  1:1354093/1‑100 (MQ=255)
                                                                                              tttACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTGGATTgag  >  2:1429372/1‑100 (MQ=255)
                                                                                              tttACAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCCCATCACATCCACTTCGGCGTGTGGAtcgcg  >  2:2596254/1‑96 (MQ=255)
                                                                                                 |                                                                                                
GTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTGGATTGAG  >  NC_000913/2118124‑2118317

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