Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,211,612 A→T 15.2% intergenic (+35/+68) icdC → / ← iraM pseudogene, isocitrate dehydrogenase C‑terminal gene fragment; idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion/RpoS stabilzer during Mg starvation, anti‑RssB factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,211,6120AT15.2% 458.7 / 65.2 171intergenic (+35/+68)icdC/iraMpseudogene, isocitrate dehydrogenase C‑terminal gene fragment; idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion/RpoS stabilzer during Mg starvation, anti‑RssB factor
Reads supporting (aligned to +/- strand):  ref base A (83/61);  new base T (12/14);  total (95/76)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 186 total aligned reads displayed.

TGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATA  >  NC_000913/1211513‑1211710
                                                                                                   |                                                                                                  
tGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATTTCTCCATGTGTTAAATATTGAAACGGGCGtata                                                                                                    <  2:3618339/100‑1 (MQ=255)
tGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtata                                                                                                    <  2:1827159/100‑1 (MQ=255)
 gATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATaa                                                                                                   >  1:969735/1‑100 (MQ=255)
  aTGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATaa                                                                                                   >  1:2814189/1‑99 (MQ=255)
   tGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAaca                                                                                                 <  1:3759864/100‑1 (MQ=255)
    ggAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAacac                                                                                                >  1:1509201/1‑100 (MQ=255)
    ggAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAacac                                                                                                >  1:507967/1‑100 (MQ=255)
      aaGGCGCTAAGCTGCTCAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCACGTCTTAAAAATTGCAACGGGCGTATAACACGc                                                                                              <  2:2300222/100‑1 (MQ=25)
       aGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtatagatcgga                                                                                             >  2:1681943/1‑93 (MQ=255)
       aGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGcc                                                                                             <  2:186450/100‑1 (MQ=255)
        ggCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGccc                                                                                            >  1:425355/1‑100 (MQ=255)
            cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgctg                                                                                        >  2:1164506/1‑100 (MQ=255)
             tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                       <  1:794580/100‑3 (MQ=25)
                gctgctGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttcaat                                                                                    <  1:1871497/100‑6 (MQ=25)
                gctgctGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGtttt                                                                                    >  1:1407348/1‑100 (MQ=255)
                gctgctGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGtttt                                                                                    <  2:2049622/100‑1 (MQ=255)
                 ctgctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTttta                                                                                   >  1:1553693/1‑100 (MQ=255)
                  tgctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAGAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttat                                                                                  <  2:576356/100‑1 (MQ=255)
                    ctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattt                                                                                <  1:3854944/100‑1 (MQ=255)
                     tgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattta                                                                               >  2:3583891/1‑100 (MQ=255)
                         aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGGGGAGGGTTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                           <  1:1111819/100‑1 (MQ=37)
                         aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGAGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                           <  2:84746/100‑1 (MQ=255)
                          tGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTgg                                                                          >  2:1241775/1‑100 (MQ=255)
                           gTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGa                                                                         >  2:3791915/1‑100 (MQ=255)
                            ttCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGat                                                                        >  1:3029729/1‑100 (MQ=255)
                             tCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGata                                                                       <  2:891354/100‑1 (MQ=255)
                             tCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGat                                                                        >  1:1944538/1‑99 (MQ=255)
                              cAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGatat                                                                      <  1:2642551/100‑1 (MQ=255)
                               agagTTTGGTGAAGCGCTCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGACGTATAACACGCCCGTTGTTTTATTTATGTGGat                                                                        >  2:243856/1‑97 (MQ=255)
                                gagTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGAtatta                                                                    >  1:2933819/1‑100 (MQ=255)
                                 agTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGAtattat                                                                   >  2:2831574/1‑100 (MQ=255)
                                 agTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGAtattat                                                                   <  2:3034882/100‑1 (MQ=255)
                                  gTTTTGTGAAGCGATCATCGAAAACATTTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt                                                                  <  2:3225536/100‑1 (MQ=255)
                                  gTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt                                                                  <  1:3491400/100‑1 (MQ=255)
                                  gTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCAATGTGTTACATATTGAAACGGGCGAATAACACGCCCGTTATTTTATTTATGTGGATcttatt                                                                  >  2:3380998/1‑100 (MQ=38)
                                    ttGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTaa                                                                <  1:1167011/100‑1 (MQ=255)
                                      ggTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATa                                                              >  2:3147938/1‑100 (MQ=255)
                                      ggTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATa                                                              >  1:720313/1‑100 (MQ=255)
                                       gTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAg                                                             <  1:2604857/100‑1 (MQ=255)
                                       gTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAg                                                             >  1:1579528/1‑100 (MQ=255)
                                        tGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAACCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGc                                                            <  1:196428/100‑1 (MQ=255)
                                         gAAGCGATCATCTAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGc                                                            <  2:1248192/99‑1 (MQ=255)
                                         gAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCa                                                           <  2:3747780/100‑1 (MQ=255)
                                           aGCAATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGCCATATAACACGCCCGTTATTTTATTTATGTGGATATTATTAATAGCa                                                           <  2:1437980/98‑1 (MQ=255)
                                               atcatcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGa                                                     >  2:2019232/1‑100 (MQ=255)
                                               atcatcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGa                                                     >  2:1319243/1‑100 (MQ=255)
                                                  atcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCa                                                  <  2:3471649/100‑1 (MQ=255)
                                                  atcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCa                                                  >  1:3836594/1‑100 (MQ=255)
                                                   tcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCat                                                 >  1:2965649/1‑100 (MQ=255)
                                                    cGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGc                                                   <  2:513988/97‑1 (MQ=255)
                                                    cGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCata                                                >  1:1258403/1‑100 (MQ=255)
                                                     gAAAACATGTAATCTCTCCACGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCatat                                               >  2:3109351/1‑100 (MQ=255)
                                                        aaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATCGCATATCGAGCATATTta                                            <  2:3294460/100‑1 (MQ=255)
                                                        aaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTta                                            <  1:3294084/100‑1 (MQ=255)
                                                          cATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTtata                                          <  1:2015972/100‑1 (MQ=255)
                                                            tGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATg                                        >  1:836338/1‑100 (MQ=255)
                                                            tGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATg                                        <  1:1672616/100‑1 (MQ=255)
                                                               aaTCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAg                                     >  1:3492154/1‑100 (MQ=255)
                                                               aaTCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAg                                     >  1:2443180/1‑100 (MQ=255)
                                                                aTCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGc                                    <  1:2853902/100‑1 (MQ=255)
                                                                aTCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGc                                    >  1:663409/1‑100 (MQ=255)
                                                                 tctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGcc                                   <  2:939775/100‑1 (MQ=255)
                                                                 tctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGcc                                   <  1:3187296/100‑1 (MQ=255)
                                                                  ctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGcc                                   >  2:2754675/1‑99 (MQ=255)
                                                                      ccATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGGTATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTa                              >  2:384764/1‑100 (MQ=255)
                                                                      ccATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTa                              >  2:737620/1‑100 (MQ=255)
                                                                       cATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCTAGCATATTTATATGAAGCACATTAc                             >  2:436866/1‑100 (MQ=255)
                                                                       cATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTAc                             >  1:2305764/1‑100 (MQ=255)
                                                                       cATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTAc                             >  1:1377401/1‑100 (MQ=255)
                                                                             ttAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGc                       >  2:1356523/1‑100 (MQ=255)
                                                                              tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                      >  1:578666/1‑100 (MQ=255)
                                                                              tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                      <  1:1793369/100‑1 (MQ=255)
                                                                              tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                      >  1:3911574/1‑100 (MQ=255)
                                                                                aataaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCa                    >  2:3895734/9‑100 (MQ=37)
                                                                                 ataaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCat                   >  1:4865/8‑100 (MQ=37)
                                                                                  tatTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCACATTACTTGAGCCCata                  >  2:485587/1‑100 (MQ=255)
                                                                                   atTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCatat                 >  2:1104364/1‑100 (MQ=255)
                                                                                   aaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCatat                 >  2:267431/6‑100 (MQ=38)
                                                                                   aaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCatat                 >  1:3870063/6‑100 (MQ=38)
                                                                                    ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATg                >  2:1211760/1‑100 (MQ=255)
                                                                                     tGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATgg               >  2:1698627/1‑100 (MQ=255)
                                                                                     aacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATgg               >  1:369631/4‑100 (MQ=38)
                                                                                     aacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATgg               >  1:1097226/4‑100 (MQ=38)
                                                                                      gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg              <  1:1846719/100‑1 (MQ=255)
                                                                                      gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg              <  1:3260112/100‑1 (MQ=255)
                                                                                      gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg              <  2:948612/100‑1 (MQ=255)
                                                                                      gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATgg               >  2:578061/1‑99 (MQ=255)
                                                                                      gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATgg               >  2:650040/1‑99 (MQ=255)
                                                                                      acaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg              >  2:2620773/3‑100 (MQ=255)
                                                                                       caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc             >  1:2806332/2‑100 (MQ=255)
                                                                                        aaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa            >  1:2312565/1‑100 (MQ=255)
                                                                                        aaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa            >  1:633590/1‑100 (MQ=255)
                                                                                         aCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCat           >  2:3694829/1‑100 (MQ=255)
                                                                                         aCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAACATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCat           >  1:1711062/1‑100 (MQ=255)
                                                                                          cGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCata          >  1:3669784/1‑100 (MQ=255)
                                                                                              cGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAtttt      >  1:1660204/1‑100 (MQ=255)
                                                                                                tataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTta    <  1:1384261/100‑1 (MQ=255)
                                                                                                tataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTta    >  2:2792289/1‑100 (MQ=255)
                                                                                                tataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTta    <  1:2241510/100‑1 (MQ=255)
                                                                                                  taACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtata  <  1:336381/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATA  >  NC_000913/1211513‑1211710

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: