Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,211,612 A→T 15.6% 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.6% 407.4 / 61.0 154intergenic (+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 (61/69);  new base T (14/10);  total (75/79)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 164 total aligned reads displayed.

TGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAA  >  NC_000913/1211513‑1211711
                                                                                                   |                                                                                                   
tGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtata                                                                                                     >  2:3096231/1‑100 (MQ=255)
 gATGGAAGGCGCTAAGCTGCTGAAATGTTCAGGGTTTGGTGAAGCGATCATCGAAAACATGTAATCTTTCCATGTGTTAATTATTGAAACGGGTGTATaa                                                                                                    >  2:2534681/1‑100 (MQ=21)
 gATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATaa                                                                                                    >  1:2692533/1‑100 (MQ=255)
  aTGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAac                                                                                                   <  2:172704/100‑1 (MQ=255)
   tGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAaca                                                                                                  >  2:1352303/1‑100 (MQ=255)
    ggAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAacac                                                                                                 <  1:2519178/100‑1 (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:395580/1‑100 (MQ=255)
          cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAGCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgt                                                                                           <  2:331326/100‑1 (MQ=25)
           gcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtt                                                                                          <  1:489000/100‑1 (MQ=25)
            cTAAGCTGCTGAAATTTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                         <  2:890217/100‑1 (MQ=255)
            cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCCTCGAAAACATGTAATCTTTTCATGTGTTAAATATTGAAACGGGCGTATAAAACGcc                                                                                              >  2:2243459/1‑95 (MQ=25)
            cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                         <  2:2016920/100‑1 (MQ=255)
            cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                         <  2:1576700/100‑1 (MQ=255)
            cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                         <  1:758814/100‑1 (MQ=255)
             tAATCTGCTGAAATGTTCAGAGTTTAGCGGAGTGCTAATCGAAAACATGTAATCACTCCATGTGTTAGATATTGAAACGTGGGTACAACACGTCCg                                                                                            <  2:753303/96‑1 (MQ=16)
             tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                        <  2:1649372/100‑3 (MQ=25)
             tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                        <  1:1080517/100‑3 (MQ=25)
             tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                        <  1:1924148/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)
                   gctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttatt                                                                                  <  2:1339062/100‑1 (MQ=255)
                    ctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattt                                                                                 >  1:2269996/1‑100 (MQ=255)
                    ctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattt                                                                                 <  1:2929438/100‑1 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                               >  1:2329210/1‑100 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                               >  1:2054636/1‑100 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                               >  2:2075855/1‑100 (MQ=255)
                       aaaTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtg                                                                              <  2:917693/100‑1 (MQ=255)
                        aaTGTTCAGAGTTAGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGCATAACACGCCCGATGTTTTATTTAtgt                                                                             >  1:873541/1‑100 (MQ=255)
                         aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                            <  1:2616901/100‑1 (MQ=255)
                         aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                            <  1:1022542/100‑1 (MQ=255)
                         aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAAAATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                            >  2:530696/1‑100 (MQ=255)
                           gTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGa                                                                          >  1:3041740/1‑100 (MQ=255)
                              cAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGatat                                                                       >  2:936567/1‑100 (MQ=255)
                               agagTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGAtatt                                                                      >  1:611091/1‑100 (MQ=255)
                                  gTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt                                                                   >  1:3028640/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)
                                         gAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCa                                                            >  1:94088/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)
                                                    atctaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCata                                                 <  2:1454157/96‑1 (MQ=255)
                                                    atctaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCata                                                 <  2:1453385/96‑1 (MQ=255)
                                                      aaaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATAtt                                               >  1:1559592/1‑100 (MQ=255)
                                                        aaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCataagat                                             >  1:1453393/1‑96 (MQ=255)
                                                           aTGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCTAGCATATTTatat                                          >  2:2515282/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)
                                                               aaTCTCTCCATGTGTTAAATATTTAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAg                                      <  2:2617630/100‑1 (MQ=255)
                                                                aTCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGc                                     <  1:2431949/100‑1 (MQ=255)
                                                                 tctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGcc                                    >  1:3478196/1‑100 (MQ=255)
                                                                 tctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGcc                                    >  1:3385364/1‑100 (MQ=255)
                                                                 tctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGcc                                    >  1:481116/1‑100 (MQ=255)
                                                                  ctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGccc                                   >  1:2196543/1‑100 (MQ=255)
                                                                   tctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCa                                  <  1:2990748/100‑1 (MQ=255)
                                                                    ctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCAt                                 >  2:2453420/1‑100 (MQ=255)
                                                                    ctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCAt                                 >  1:1735108/1‑100 (MQ=255)
                                                                       cATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTAc                              <  2:3440349/100‑1 (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:2479073/100‑1 (MQ=255)
                                                                              tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                       >  1:3607723/1‑100 (MQ=255)
                                                                               aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGccc                      <  2:811463/100‑1 (MQ=255)
                                                                                aataaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCa                     >  1:3635429/9‑100 (MQ=37)
                                                                                  tatTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCGTATTTATATGAAGCCCATTACTTGAGCCCata                   <  1:3258293/100‑1 (MQ=255)
                                                                                  taaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCata                   >  1:564345/7‑100 (MQ=38)
                                                                                  taaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCata                   >  1:1661659/7‑100 (MQ=38)
                                                                                   atTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCTATTACTTGAGCCCatat                  >  2:547935/1‑100 (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              >  2:2507559/2‑100 (MQ=255)
                                                                                       caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc              >  2:105053/2‑100 (MQ=255)
                                                                                       caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc              >  2:453590/2‑100 (MQ=255)
                                                                                        aaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa             >  2:2883214/1‑100 (MQ=255)
                                                                                        aaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa             <  1:1331167/100‑1 (MQ=255)
                                                                                        aaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa             >  2:30224/1‑100 (MQ=255)
                                                                                        aaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa             <  1:936574/100‑1 (MQ=255)
                                                                                         aCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCat            >  2:3561202/1‑100 (MQ=255)
                                                                                            ggCGTATAACACGCCCGTTGTTTTATTTATGTGGATGTTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAtt         >  2:266324/1‑100 (MQ=255)
                                                                                             gCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttt        >  1:1382286/1‑100 (MQ=255)
                                                                                              cGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATGTGAAGCCCATTACTTGAGCCCATATGGGCATAtttt       >  1:3469252/1‑100 (MQ=255)
                                                                                               gTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttttt      <  2:473972/100‑1 (MQ=255)
                                                                                               gTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttttt      <  2:585558/100‑1 (MQ=255)
                                                                                                 ataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtat    >  2:2675605/1‑100 (MQ=255)
                                                                                                 ataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtat    >  2:1886470/1‑100 (MQ=255)
                                                                                                   aaCACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCTCATTACTTGAGCCCATATGGGCATATTTTTATaa  >  2:1742221/1‑100 (MQ=255)
                                                                                                   aaCACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATaa  >  1:2809505/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAA  >  NC_000913/1211513‑1211711

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: