Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,211,614 C→T 15.4% intergenic (+37/+66) 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,6140CT15.4% 397.1 / 52.9 156intergenic (+37/+66)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 C (64/67);  new base T (14/10);  total (79/77)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

Only 100 of 164 total aligned reads displayed.

ATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATG  >  NC_000913/1211515‑1211713
                                                                                                   |                                                                                                   
 tGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAaca                                                                                                    >  2:1352303/1‑100 (MQ=255)
    aaGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGc                                                                                                 >  2:979577/1‑100 (MQ=255)
      ggCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGccc                                                                                               >  2:620002/1‑100 (MQ=255)
       gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCg                                                                                              >  1:1633518/1‑100 (MQ=255)
        cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgt                                                                                             >  1:2864860/1‑100 (MQ=255)
        cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgt                                                                                             >  1:3389724/1‑100 (MQ=255)
        cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgt                                                                                             >  1:395580/1‑100 (MQ=255)
        cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAGCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgt                                                                                             <  2:331326/100‑1 (MQ=25)
         gcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGTGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAATCGGTCGTGTTATACGCCCgtt                                                                                            <  2:1506137/100‑1 (MQ=21)
         gcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtt                                                                                            <  1:489000/100‑1 (MQ=25)
          cTAAGCTGCTGAAATTTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                           <  2:890217/100‑1 (MQ=255)
          cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                           <  1:758814/100‑1 (MQ=255)
          cTAAGCTGCGGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                           <  2:2017942/100‑1 (MQ=255)
           tAATCTGCTGAAATGTTCAGAGTTTAGCGGAGTGCTAATCGAAAACATGTAATCACTCCATGTGTTAGATATTGAAACGTGGGTACAACACGTCCg                                                                                              <  2:753303/96‑1 (MQ=16)
           tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                          <  2:1649372/100‑3 (MQ=25)
           tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                          <  2:2271028/100‑3 (MQ=25)
           tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                          <  1:1080517/100‑3 (MQ=25)
           tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                          <  1:1924148/100‑3 (MQ=25)
           tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATAGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                          <  2:2900757/100‑3 (MQ=25)
             aGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttt                                                                                        <  1:1447114/100‑1 (MQ=255)
              gctgctGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGtttt                                                                                       >  2:2430238/1‑100 (MQ=255)
              gctgctGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGtttt                                                                                       <  1:1328286/100‑1 (MQ=255)
              gctgctGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGtttt                                                                                       <  1:2540074/100‑1 (MQ=255)
                 gctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttatt                                                                                    <  2:1339062/100‑1 (MQ=255)
                  ctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattt                                                                                   >  1:2269996/1‑100 (MQ=255)
                    gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                                 >  1:2329210/1‑100 (MQ=255)
                    gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                                 >  2:2075855/1‑100 (MQ=255)
                    gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                                 >  1:2054636/1‑100 (MQ=255)
                     aaaTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtg                                                                                <  2:917693/100‑1 (MQ=255)
                      aaTGTTCAGAGTTAGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGCATAACACGCCCGATGTTTTATTTAtgt                                                                               >  1:873541/1‑100 (MQ=255)
                       aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                              <  1:3045869/100‑1 (MQ=255)
                       aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                              <  1:1022542/100‑1 (MQ=255)
                       aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                              <  1:3528089/100‑1 (MQ=255)
                       aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAAAATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                              >  2:530696/1‑100 (MQ=255)
                            cAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGatat                                                                         >  2:936567/1‑100 (MQ=255)
                                 tttGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATAttatta                                                                    >  1:2757240/1‑100 (MQ=255)
                                  ttGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTaa                                                                   >  1:1270817/1‑100 (MQ=255)
                                    ggTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAATTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATa                                                                 <  2:2333498/100‑1 (MQ=255)
                                     gTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAg                                                                >  1:2353628/1‑100 (MQ=255)
                                           cGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATc                                                          <  2:1997988/100‑1 (MQ=255)
                                            gATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCg                                                         <  2:2028350/100‑1 (MQ=255)
                                               catcGAAAACATGTAATCTCTCCATTTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGc                                                      <  2:1525537/100‑1 (MQ=255)
                                               catcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGc                                                      >  1:2059334/1‑100 (MQ=255)
                                                atcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCa                                                     >  2:1737853/1‑100 (MQ=255)
                                                  atctaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCata                                                   <  2:1453385/96‑1 (MQ=255)
                                                    aaaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATAtt                                                 <  2:2531189/100‑1 (MQ=255)
                                                    aaaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATAtt                                                 >  1:1559592/1‑100 (MQ=255)
                                                      aaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCataagat                                               >  1:1453393/1‑96 (MQ=255)
                                                         aTGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTatat                                            >  2:2665072/1‑100 (MQ=255)
                                                         aTGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTatat                                            >  2:3385889/1‑100 (MQ=255)
                                                         aTGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTatat                                            >  2:61300/1‑100 (MQ=255)
                                                         aTGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTatat                                            <  2:1755017/100‑1 (MQ=255)
                                                           gacaTCTCTCCATGTGTTAAATATTGAAACGGGCTTATACCACGCCCGTTGTTTTATTTATGTGGATATGATTAACAGCATATCGGGCATATTTATATGa                                          <  2:672988/97‑1 (MQ=255)
                                                           gTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGa                                          >  1:1772923/1‑100 (MQ=255)
                                                              aTCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGc                                       <  1:2431949/100‑1 (MQ=255)
                                                               tctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGcc                                      >  1:3478196/1‑100 (MQ=255)
                                                               tctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGcc                                      >  1:3385364/1‑100 (MQ=255)
                                                                ctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGccc                                     >  1:2196543/1‑100 (MQ=255)
                                                                 tctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCa                                    <  1:2990748/100‑1 (MQ=255)
                                                                  ctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCAt                                   >  1:1735108/1‑100 (MQ=255)
                                                                      aTGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACt                               <  1:1221103/100‑1 (MQ=255)
                                                                       tgtgTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACtt                              <  1:2476029/100‑1 (MQ=255)
                                                                        gtgtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTg                             <  2:1670221/100‑1 (MQ=255)
                                                                          gtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAg                           >  1:311512/1‑100 (MQ=255)
                                                                            tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                         <  2:15989/100‑1 (MQ=255)
                                                                            tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                         <  2:2479073/100‑1 (MQ=255)
                                                                            tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                         >  1:2940662/1‑100 (MQ=255)
                                                                             aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGccc                        <  2:811463/100‑1 (MQ=255)
                                                                              aataaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCa                       >  1:3635429/9‑100 (MQ=37)
                                                                                taaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCata                     >  1:1661659/7‑100 (MQ=38)
                                                                                taaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCata                     >  1:564345/7‑100 (MQ=38)
                                                                                 atTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCTATTACTTGAGCCCatat                    >  2:547935/1‑100 (MQ=255)
                                                                                  ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCatat                    >  2:245961/1‑99 (MQ=255)
                                                                                  aaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATg                   >  1:114563/5‑100 (MQ=38)
                                                                                   aacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCatat                    >  1:2778799/4‑98 (MQ=38)
                                                                                    gAAACGGGCGTATAACACGTCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg                 <  2:3260234/100‑1 (MQ=255)
                                                                                    gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCAATTACTTTAGCCCATATg                   >  2:1559374/1‑98 (MQ=255)
                                                                                    acaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg                 >  2:3545222/3‑100 (MQ=255)
                                                                                     caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc                >  1:417177/2‑100 (MQ=255)
                                                                                     caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc                >  2:105053/2‑100 (MQ=255)
                                                                                     caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc                >  2:453590/2‑100 (MQ=255)
                                                                                     caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc                >  2:2507559/2‑100 (MQ=255)
                                                                                     aaaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc                <  1:981738/100‑1 (MQ=255)
                                                                                      aaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa               >  2:30224/1‑100 (MQ=255)
                                                                                      aaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa               <  1:3032095/100‑1 (MQ=255)
                                                                                      aaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa               <  1:1331167/100‑1 (MQ=255)
                                                                                      aaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa               <  1:936574/100‑1 (MQ=255)
                                                                                          ggCGTATAACACGCCCGTTGTTTTATTTATGTGGATGTTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAtt           >  2:266324/1‑100 (MQ=255)
                                                                                           gCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttt          >  1:1382286/1‑100 (MQ=255)
                                                                                             gTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttttt        <  2:2593485/100‑1 (MQ=255)
                                                                                             gTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttttt        <  2:473972/100‑1 (MQ=255)
                                                                                             gTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttttt        <  2:585558/100‑1 (MQ=255)
                                                                                               ataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtat      >  2:1886470/1‑100 (MQ=255)
                                                                                                taACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtata     >  1:580215/1‑100 (MQ=255)
                                                                                                 aaCACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCTCATTACTTGAGCCCATATGGGCATATTTTTATaa    >  2:1742221/1‑100 (MQ=255)
                                                                                                  acacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGGATATTTATATGAAGCACATAACTTGAGCCCAAATGGGCATATTTTtatca    >  2:2815128/1‑97 (MQ=255)
                                                                                                  acacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAAt   >  1:3293402/1‑100 (MQ=255)
                                                                                                   cacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGCGCATATTTATATGAAGCGCATTACTTGAGCCCATATGGGCATATTTTTATAATg  >  1:2498342/1‑100 (MQ=255)
                                                                                                   cacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATg  >  1:869753/1‑100 (MQ=255)
                                                                                                   cacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGACCATATGGGCATATTTTTATAATg  >  1:2535065/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATG  >  NC_000913/1211515‑1211713

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: