Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,211,612 A→T 17.7% 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,6120AT17.7% 461.6 / 86.5 191intergenic (+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 (82/74);  new base T (20/14);  total (103/88)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

Only 100 of 203 total aligned reads displayed.

TGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATA  >  NC_000913/1211513‑1211710
                                                                                                   |                                                                                                  
 gATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATaa                                                                                                   >  1:2062237/1‑100 (MQ=255)
  aTGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAac                                                                                                  <  1:1455954/100‑1 (MQ=255)
    ggAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAacac                                                                                                >  2:3690007/1‑100 (MQ=255)
     gAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACg                                                                                               >  2:192529/1‑100 (MQ=255)
       aGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGcc                                                                                             <  2:1911805/100‑1 (MQ=255)
       aGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGcc                                                                                             <  2:2654631/100‑1 (MQ=255)
       aGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGcc                                                                                             <  1:1861456/100‑1 (MQ=255)
         gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAATCGATCATCGAAAACATGGAATCACCCCCCCTGTTAAATATTGAAACGGGCGTATAACACGcc                                                                                             <  1:1371064/98‑1 (MQ=21)
          cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgt                                                                                          <  2:1105159/100‑1 (MQ=25)
          cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgt                                                                                          >  2:1473257/1‑100 (MQ=255)
           gcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtt                                                                                         <  1:3409310/100‑1 (MQ=25)
             tAAGCTGCTGAAATGTTCAGAGTTTGTTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGACACGGGCGTGTTATACGCCCgtttc                                                                                       <  1:205134/100‑3 (MQ=25)
             tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                       <  2:1769484/100‑3 (MQ=25)
             tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttgt                                                                                       <  2:971195/100‑1 (MQ=255)
             tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttgt                                                                                       <  2:1173975/100‑1 (MQ=255)
               aGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttcaa                                                                                     <  2:687768/100‑5 (MQ=25)
               aGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttt                                                                                     <  2:2128865/100‑1 (MQ=255)
                 ctgctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTttta                                                                                   <  1:1754017/100‑1 (MQ=255)
                  tgctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttat                                                                                  >  2:288356/1‑100 (MQ=255)
                   gctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttatt                                                                                 <  2:126835/100‑1 (MQ=255)
                    ctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattt                                                                                <  2:840675/100‑1 (MQ=255)
                    ctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattt                                                                                >  1:296692/1‑100 (MQ=255)
                     tgAAATGTTCAGAGTTTGGTGCAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAAGGGGCGTATATCACGCCCGTTGTtttattt                                                                                >  1:3021941/1‑99 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                              >  1:2453957/1‑100 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                              >  1:3992623/1‑100 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                              >  1:2572054/1‑100 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                              >  1:729427/1‑100 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                              >  1:2658683/1‑100 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                              >  1:839372/1‑100 (MQ=255)
                           gTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACTGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                           <  1:4382410/98‑1 (MQ=255)
                            ttCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGat                                                                        >  1:4111583/1‑100 (MQ=255)
                             tCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGata                                                                       <  2:2839089/100‑1 (MQ=255)
                             tCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGata                                                                       >  2:240287/1‑100 (MQ=255)
                               agagTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGAtatt                                                                     >  2:1498717/1‑100 (MQ=255)
                                gagTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACCCGCCCGTTGTTTTATTTATGTGGAtatta                                                                    >  1:1222097/1‑100 (MQ=255)
                                  gTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATTTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt                                                                  <  2:195906/100‑1 (MQ=255)
                                      ggTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATa                                                              >  2:580319/1‑100 (MQ=255)
                                       gTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAg                                                             >  2:2671152/1‑100 (MQ=255)
                                       gTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAg                                                             >  1:1672486/1‑100 (MQ=255)
                                        tGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGc                                                            <  1:4191396/100‑1 (MQ=255)
                                        tGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGc                                                            <  1:544929/100‑1 (MQ=255)
                                         gAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGc                                                            <  2:3510631/99‑1 (MQ=255)
                                         gAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCa                                                           >  2:2178164/1‑100 (MQ=255)
                                         gAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCa                                                           <  2:3431810/100‑1 (MQ=255)
                                           aGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCata                                                         <  2:3470392/100‑1 (MQ=255)
                                             cGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCatatn                                                       <  2:1929566/100‑2 (MQ=255)
                                              gATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCg                                                      >  2:595949/1‑100 (MQ=255)
                                                tcatcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAg                                                    <  1:2064067/100‑1 (MQ=255)
                                                tcatcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAg                                                    <  2:2554694/100‑1 (MQ=255)
                                                 catcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTTTTTATGTGGATATTATTAATAGCATATCGAGc                                                   <  2:4172306/100‑1 (MQ=255)
                                                  atcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGc                                                   <  1:1312459/99‑1 (MQ=255)
                                                  atcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCa                                                  <  2:2284285/100‑1 (MQ=255)
                                                  atcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCa                                                  <  2:541024/100‑1 (MQ=255)
                                                     gAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGCGCatat                                               >  1:1125220/1‑100 (MQ=255)
                                                      aaaaCATGTAATCTCTCCATGGGTTGAATATTTAAACGGGCGTATAACACGACAGTTGTATTATTTAGGTGGATATTATTAATAGCATATCGaagatct                                               >  1:4284136/1‑93 (MQ=255)
                                                          cATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTtata                                          <  1:1051106/100‑1 (MQ=255)
                                                          cATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTtata                                          <  1:3683927/100‑1 (MQ=255)
                                                           atttAATCTCTCCATGTGTTAAATATTTAAACTGGCTTATAACACTCCCGTTGTTTTATTTATTTTGATATTATTAATAGCATATCGAGCATATTTatat                                         <  1:2532533/97‑1 (MQ=255)
                                                            tGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATg                                        <  2:79726/100‑1 (MQ=255)
                                                            tGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATg                                        >  2:3056233/1‑100 (MQ=255)
                                                             gTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGa                                       >  2:3793944/1‑100 (MQ=255)
                                                              tAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGaa                                      >  2:3814413/1‑100 (MQ=255)
                                                               aaTCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAg                                     >  2:464620/1‑100 (MQ=255)
                                                                 tctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTGTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGc                                    >  2:3632691/1‑99 (MQ=255)
                                                                    ctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCAt                                <  1:958635/100‑1 (MQ=255)
                                                                       cATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTAc                             <  2:4380080/100‑1 (MQ=255)
                                                                       cATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTAc                             <  1:3807598/100‑1 (MQ=255)
                                                                           tgtTAAATATTGAAATGTGTTTATAATAAGTATGTTGTTTCATTTATGTTGATATTATTAATATCATATCGATTATATTTATATGAAGCCCATTACt                            <  1:2807309/97‑1 (MQ=255)
                                                                            tttAAATAGTGAAACGGGCGTGTAACACGCCCGTTGGTTTATTAATGTGGATATTATTAATAGCATATCTAGCATTTTTATATGATGCCCAGTAGTTGa                         >  1:2622811/2‑99 (MQ=255)
                                                                            gtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAg                        <  1:1792275/100‑1 (MQ=255)
                                                                            gtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTGTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGa                         >  1:2339230/1‑99 (MQ=255)
                                                                            gtTAAATATTGAAACAGGCGTATAACACGCTCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAg                        <  2:4334003/100‑1 (MQ=255)
                                                                             ttAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGc                       <  2:3272879/100‑1 (MQ=255)
                                                                              tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGCAGCCCATTACTTGAGcc                      >  2:4239203/1‑100 (MQ=255)
                                                                              tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                      >  2:2038504/1‑100 (MQ=255)
                                                                              tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                      >  1:3245828/1‑100 (MQ=255)
                                                                              tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATAGAGCATATTTATATGAAGCCCATTACTTGa                         >  1:3144962/1‑97 (MQ=255)
                                                                              tAAATANTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTTTTAATAGCATATCGAGCGTATTTATATGAAGCCNATTACTTGAGc                       >  1:2565080/1‑99 (MQ=255)
                                                                               aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGccc                     >  2:654791/1‑100 (MQ=255)
                                                                                aaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCa                    >  2:1980099/1‑100 (MQ=255)
                                                                                 ataaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCat                   >  2:2376903/8‑100 (MQ=37)
                                                                                  tatTGAATCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCata                  <  1:3011165/100‑1 (MQ=255)
                                                                                  taaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCata                  >  2:3403144/7‑100 (MQ=38)
                                                                                    ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATg                <  1:1671814/100‑1 (MQ=255)
                                                                                    aaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATg                >  2:2223837/5‑100 (MQ=38)
                                                                                      gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg              <  1:439767/100‑1 (MQ=255)
                                                                                      acaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCCTATg                >  1:2876340/3‑98 (MQ=38)
                                                                                      acaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg              >  2:2755130/3‑100 (MQ=255)
                                                                                       caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATgg               >  2:1046116/2‑98 (MQ=255)
                                                                                       caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc             >  2:496622/2‑100 (MQ=255)
                                                                                       caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc             >  2:2722154/2‑100 (MQ=255)
                                                                                       aaaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc             >  2:918176/1‑100 (MQ=255)
                                                                                        aaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa            >  2:1364912/1‑100 (MQ=255)
                                                                                          cGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCata          <  1:3811226/100‑1 (MQ=255)
                                                                                          cGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCata          <  1:3152561/100‑1 (MQ=255)
                                                                                               gTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttttt     >  1:4169965/1‑100 (MQ=255)
                                                                                                tataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTta    <  2:3992650/100‑1 (MQ=255)
                                                                                                tataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTta    <  2:1346753/100‑1 (MQ=255)
                                                                                                 ataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtat   <  2:3003211/100‑1 (MQ=255)
                                                                                                  taACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtata  >  2:1380564/1‑100 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑