Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,756,562:1 +G 15.3% intergenic (‑51/+114) aldB ← / ← yiaY aldehyde dehydrogenase B/L‑threonine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,756,5621.G15.3% 396.0 / 76.2 124intergenic (‑51/+114)aldB/yiaYaldehyde dehydrogenase B/L‑threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (58/47);  new base G (12/7);  total (70/54)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-02

Only 100 of 130 total aligned reads displayed.

AACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGT  >  NC_000913/3756463‑3756661
                                                                                                    |                                                                                                   
aaCCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGaa                                                                                                      <  1:2455526/100‑1 (MQ=255)
 aCCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGaa‑a                                                                                                    <  2:1588629/100‑1 (MQ=255)
 aCCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGaa‑a                                                                                                    <  1:212646/100‑1 (MQ=255)
   cATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATc                                                                                                  <  2:95387/100‑1 (MQ=255)
   cATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATc                                                                                                  >  2:1960435/1‑100 (MQ=255)
     tACTCTCCGGGCTTAATCTGTGCTGAAGTGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTTAGGTATAACAAGGGTAGACGTGACTTGCGAA‑ATCtt                                                                                                <  1:357444/100‑1 (MQ=255)
      aCTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTc                                                                                               <  1:41998/100‑1 (MQ=255)
       cTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCg                                                                                              <  2:1569744/100‑1 (MQ=255)
       cTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCg                                                                                              >  2:597310/1‑100 (MQ=255)
       cTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCg                                                                                              <  2:517933/100‑1 (MQ=255)
       cTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCg                                                                                              >  2:1689349/1‑100 (MQ=255)
         cGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGaagatttcgcc                                                                                             <  2:561622/100‑10 (MQ=255)
         cGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGcc                                                                                            <  2:1553028/100‑1 (MQ=255)
          gCCGGGCTTAATCTGTTCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCa                                                                                           >  1:1014054/1‑100 (MQ=255)
           ccGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGAT‑TTCGCCAg                                                                                          <  2:1379724/100‑1 (MQ=255)
           ccGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑GATTTCGCCAg                                                                                          <  1:742430/100‑1 (MQ=255)
            cGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGAT‑TTCGCCAGt                                                                                         <  1:1247525/100‑1 (MQ=255)
            cGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑GATTTCGCCAGt                                                                                         <  2:159105/100‑1 (MQ=255)
             gggCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGAT‑TTCGCCAGTc                                                                                        <  2:1834181/100‑1 (MQ=255)
             gggCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTc                                                                                        <  1:2579078/100‑1 (MQ=255)
              ggCGTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGAT‑TTCGCCAGTca                                                                                       <  1:2072743/100‑2 (MQ=255)
              ctcTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGCGACTGGCGAAGAT‑TTCGCCAGTca                                                                                       <  1:1525030/98‑2 (MQ=255)
               gCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGAT‑TTCGCCAGTcac                                                                                      <  2:2119933/100‑3 (MQ=255)
               gCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑GATTTCGCCAGTcan                                                                                      <  1:1071373/100‑3 (MQ=255)
                cTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGAT‑TTCGCCAGTcacg                                                                                     <  1:518390/100‑4 (MQ=39)
                cTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGt                                                                                     >  1:1848375/1‑100 (MQ=255)
                cTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGCGACTGGCGAAGAT‑TTCGCCAGTcaag                                                                                     >  2:1525039/1‑97 (MQ=255)
                  tAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTaa                                                                                   <  2:256064/100‑1 (MQ=255)
                   aaTCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAAc                                                                                  <  2:431811/100‑1 (MQ=255)
                   aaTCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAAc                                                                                  >  1:945583/1‑100 (MQ=255)
                         tgCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCttt                                                                            <  1:1732465/100‑1 (MQ=255)
                         tgCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCttt                                                                            >  2:1170812/1‑100 (MQ=255)
                            tGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACg                                                                         >  2:1274368/1‑100 (MQ=255)
                             gAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGa                                                                        >  2:1728767/1‑100 (MQ=255)
                                gggggATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAg                                                                     >  2:781326/1‑100 (MQ=255)
                                 ggtgATTATTAGTCATGATCGTCTCCATGCGGTGTGAGGTATAACAAGGGTAGAAGTGACTAGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGATAg                                                                     <  1:1391000/96‑1 (MQ=255)
                                  gggATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGAGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCt                                                                   >  1:1741173/1‑100 (MQ=255)
                                     attattGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATc                                                                <  1:1845857/100‑1 (MQ=255)
                                       ccgatctCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCAc                                                              <  2:84527/94‑1 (MQ=255)
                                        attGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACg                                                             >  1:2305875/1‑100 (MQ=255)
                                         ttGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGa                                                            >  2:2219927/1‑100 (MQ=255)
                                         ttGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGa                                                            >  2:1598103/1‑100 (MQ=255)
                                          tGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGaa                                                           >  2:176331/1‑100 (MQ=255)
                                           ggTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAAt                                                          >  2:1264266/1‑100 (MQ=255)
                                           ggTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAAt                                                          >  2:1903335/1‑100 (MQ=255)
                                              cATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCAcagatcgg                                                       >  1:84526/1‑93 (MQ=255)
                                                tGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACg                                                     <  2:63426/100‑1 (MQ=255)
                                                tGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACg                                                     >  1:2309044/1‑100 (MQ=255)
                                                   tcgtcgCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGc                                                  <  1:183077/100‑1 (MQ=255)
                                                   tCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGc                                                  <  2:268435/100‑1 (MQ=255)
                                                         ccTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGt                                            <  2:1343523/100‑1 (MQ=255)
                                                           ttGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTa                                          <  1:368707/100‑1 (MQ=255)
                                                              cGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACag                                       >  2:1264512/1‑100 (MQ=255)
                                                              cGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACag                                       >  2:1812104/1‑100 (MQ=255)
                                                               ggTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACaga                                      <  1:457467/100‑1 (MQ=255)
                                                               ggTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACaga                                      <  2:1953025/100‑1 (MQ=255)
                                                                gtgtGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACagag                                     >  1:2092937/1‑100 (MQ=255)
                                                                   tgAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATa                                  >  2:2276964/1‑100 (MQ=255)
                                                                   tgAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATa                                  <  2:2065462/100‑1 (MQ=255)
                                                                    gAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTGACGACAGCGATCACGAATTTACGGGCAAGGGTTCCAGAGATAg                                 >  2:1948214/1‑100 (MQ=255)
                                                                      ggTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGtt                               >  1:1252907/1‑100 (MQ=255)
                                                                      ggTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGtt                               >  2:2240073/1‑100 (MQ=255)
                                                                      ggTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGtt                               >  2:1999980/1‑100 (MQ=255)
                                                                      ggTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGtt                               <  2:1472600/100‑1 (MQ=255)
                                                                       gTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTg                              >  1:1011535/1‑100 (MQ=255)
                                                                        tataACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGa                             >  2:422777/1‑100 (MQ=255)
                                                                        tataACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGa                             <  1:2534689/100‑1 (MQ=255)
                                                                          taACAAGGGTAGACGTGACTGGCGAA‑ATCTTCACCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaa                           >  2:2459859/1‑100 (MQ=255)
                                                                           aaCAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaaa                          >  2:323622/1‑100 (MQ=255)
                                                                           aaCAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaaa                          <  2:2351408/100‑1 (MQ=255)
                                                                           aaCAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaaa                          <  1:1028906/100‑1 (MQ=255)
                                                                            aCAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaaaa                         >  2:1112975/1‑100 (MQ=255)
                                                                            aCAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaaaa                         >  1:157916/1‑100 (MQ=255)
                                                                                gggTAGACGTGACTGGCGAA‑ATCTTAGACAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCg                      <  2:354033/99‑1 (MQ=255)
                                                                                 ggTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGc                    <  1:51714/100‑1 (MQ=255)
                                                                                 gctgttaccgACTGGCGAAGAT‑TTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGc                    >  1:2133688/10‑100 (MQ=37)
                                                                                      taccgACTGGCGAAGAT‑TTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTAGAGAGATAGTTGAAAAAGCGGCTAACa               >  1:663000/5‑100 (MQ=38)
                                                                                       cGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACaa              >  1:1286849/1‑100 (MQ=255)
                                                                                        ccgACTGGCGAAGAT‑TTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAAt             >  2:495424/3‑100 (MQ=255)
                                                                                        ccgACTGGCGAAGAT‑TTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAAt             >  1:1652687/3‑100 (MQ=255)
                                                                                         cgACTGGCGAA‑GATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAAtt            >  2:1989413/2‑100 (MQ=255)
                                                                                          gACTGGCGAAGAT‑TTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAAttt           >  2:2400551/1‑100 (MQ=255)
                                                                                          gACTGGCGAAGAT‑TTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAAttt           >  2:376794/1‑100 (MQ=255)
                                                                                           acgggCGGA‑ATCTCCGCCAGTAGGTAACAGCTATAGAACAGGTATCACGAATTAACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTg          <  1:700374/97‑1 (MQ=255)
                                                                                           aCTGGCGAAGAT‑TTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTg          >  2:2352872/1‑100 (MQ=255)
                                                                                           aCTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTg          <  2:2007336/100‑1 (MQ=255)
                                                                                            cTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGc         >  2:429818/1‑100 (MQ=255)
                                                                                            cTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGc         >  2:430875/1‑100 (MQ=255)
                                                                                            cTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGc         >  2:485111/1‑100 (MQ=255)
                                                                                            cTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGc         >  2:1643704/1‑100 (MQ=255)
                                                                                            cTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGc         >  2:738520/1‑100 (MQ=255)
                                                                                             tGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGcc        <  2:2095012/100‑1 (MQ=255)
                                                                                             tGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGcc        <  2:1520165/100‑1 (MQ=255)
                                                                                             tGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGcc        <  1:1977425/100‑1 (MQ=255)
                                                                                              ggCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCa       >  2:2278008/1‑100 (MQ=255)
                                                                                              ggCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCa       >  2:542199/1‑100 (MQ=255)
                                                                                               gCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAg      >  1:879519/1‑100 (MQ=255)
                                                                                                cGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGc     >  2:46871/1‑100 (MQ=255)
                                                                                                 gAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGcc    >  1:1525833/1‑100 (MQ=255)
                                                                                                   a‑aTCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGt  >  2:608994/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
AACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGT  >  NC_000913/3756463‑3756661

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