Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,211,612 A→T 7.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,6120AT 7.2% 372.8 / 13.6 125intergenic (+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 (68/47);  new base T (3/6);  total (71/54)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 132 total aligned reads displayed.

TATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGG  >  NC_000913/1211463‑1211760
                                                                                                                                                     |                                                                                                                                                    
tATTGAAGTCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACTTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtata                                                                                                                                                      <  2:1648204/150‑1 (MQ=255)
tATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtata                                                                                                                                                      <  1:74339/150‑1 (MQ=255)
  tGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAac                                                                                                                                                    <  2:710922/150‑1 (MQ=255)
   ggAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAaca                                                                                                                                                   <  1:423142/150‑1 (MQ=255)
   ggAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAaca                                                                                                                                                   >  2:1413728/1‑150 (MQ=255)
     aaGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACg                                                                                                                                                 >  1:153328/1‑150 (MQ=255)
       ggCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGcc                                                                                                                                               >  1:1137510/1‑150 (MQ=255)
       ggCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGcc                                                                                                                                               <  1:1322495/150‑1 (MQ=255)
         cgcAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCg                                                                                                                                             >  1:530813/1‑150 (MQ=255)
            aatcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGCTTGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgttt                                                                                                                                          <  2:393205/150‑2 (MQ=25)
             atcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                                                                         <  2:495935/150‑3 (MQ=255)
             atcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttgt                                                                                                                                         >  1:257962/1‑150 (MQ=255)
              tcaatGCCAAGCCCGTAACTTATGACTTCGAACGTCTGAGGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttca                                                                                                                                        <  1:1863250/150‑4 (MQ=25)
                 atGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTttta                                                                                                                                     <  1:575118/150‑1 (MQ=255)
                   gCCATGACGGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGGGAAGCGATCATCGAAGACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttatt                                                                                                                                   <  1:1214880/150‑1 (MQ=255)
                    ccAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattt                                                                                                                                  >  1:1249027/1‑150 (MQ=255)
                     cAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattta                                                                                                                                 >  1:866062/1‑150 (MQ=255)
                     cAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattta                                                                                                                                 >  2:19481/1‑150 (MQ=255)
                       aGACCGTAACTTATGACTTCGAACGTCGGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtg                                                                                                                               <  1:1504353/150‑1 (MQ=255)
                        gACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgt                                                                                                                              >  1:1108926/1‑150 (MQ=255)
                           cGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGa                                                                                                                           <  1:409445/150‑1 (MQ=255)
                                  tATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAATCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt                                                                                                                    <  1:1942024/150‑1 (MQ=255)
                                    tGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTaa                                                                                                                  >  2:450113/1‑150 (MQ=255)
                                       cTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTTAATATTGAAACGGGCGTATAACCCGCCCGTTGTTTTATTTATGTGGATATTATTAATAg                                                                                                               >  2:1289688/1‑150 (MQ=255)
                                       cTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCAATGTGTGAAATATTGAAACGGGCGTATAACACGCCCGTGGTTTTATTTATGTGGATATTATTAATAg                                                                                                               >  1:287684/1‑150 (MQ=255)
                                           gAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCata                                                                                                           >  2:1953348/1‑150 (MQ=255)
                                              cGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCg                                                                                                        <  1:1109238/150‑1 (MQ=255)
                                               gTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGa                                                                                                       <  2:202579/150‑1 (MQ=255)
                                               gTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGa                                                                                                       >  1:1329170/1‑150 (MQ=255)
                                                tCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACTCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAg                                                                                                      >  1:1247469/1‑150 (MQ=255)
                                                 cTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGc                                                                                                     <  1:668025/150‑1 (MQ=255)
                                                  tGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCa                                                                                                    >  2:1494022/1‑150 (MQ=255)
                                                  tGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCa                                                                                                    <  2:334797/150‑1 (MQ=255)
                                                   gATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCat                                                                                                   >  1:1798445/1‑150 (MQ=255)
                                                      ggAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATAGTATTAATAGCATATCGAGCATAtt                                                                                                >  2:301214/1‑150 (MQ=255)
                                                        aaGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTta                                                                                              <  2:1284261/150‑1 (MQ=255)
                                                        aaGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTta                                                                                              >  1:63375/1‑150 (MQ=255)
                                                          ggCGCTAAGCTTCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTtata                                                                                            <  1:1963405/150‑1 (MQ=255)
                                                          ggCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTtata                                                                                            <  2:408100/150‑1 (MQ=255)
                                                            cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATg                                                                                          <  1:1589584/150‑1 (MQ=255)
                                                             gcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGa                                                                                         <  1:846930/150‑1 (MQ=255)
                                                              cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGaa                                                                                        >  2:1101406/1‑150 (MQ=255)
                                                              cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGaa                                                                                        >  1:930576/1‑150 (MQ=255)
                                                              cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGaa                                                                                        >  1:1620890/1‑150 (MQ=255)
                                                                aaGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGc                                                                                      >  2:1094700/1‑150 (MQ=255)
                                                                   ctgctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATAGAGCATATTTATATGAAGCCCa                                                                                   >  1:119012/1‑150 (MQ=255)
                                                                   atgatgAAATGTTCATAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTTAAACGGGCTTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCa                                                                                   <  1:1297852/146‑1 (MQ=255)
                                                                    tgctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCAt                                                                                  >  2:172430/1‑150 (MQ=255)
                                                                        gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATTTTGAAACGGGCGTATAACACGCCAGTTGTTTTATTTATGTGGATATTATTAATAGCTTATCGAGCATATTTATATGAATCCCATTACt                                                                              >  1:331569/1‑150 (MQ=255)
                                                                          aaTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTg                                                                            <  1:1259743/150‑1 (MQ=255)
                                                                            tGTTCAGAGTTTGGTGAAGCGATCATCTAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGTGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAg                                                                          <  1:1785550/150‑1 (MQ=255)
                                                                              ttCAGAGTTTGGTGAAGCGATCATCGAAAACATGTTATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                                                                        <  1:420365/150‑1 (MQ=255)
                                                                              ttCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                                                                        <  1:1185254/150‑1 (MQ=255)
                                                                              ttCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                                                                        <  2:708646/150‑1 (MQ=255)
                                                                                cAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCa                                                                      >  2:635483/1‑150 (MQ=255)
                                                                                cAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTAATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCa                                                                      >  2:48162/1‑150 (MQ=255)
                                                                                 agagTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCat                                                                     <  2:226700/150‑1 (MQ=255)
                                                                                  gagTTTGGTGAAGCCATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCata                                                                    <  2:341761/150‑1 (MQ=255)
                                                                                      ttGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg                                                                <  2:281316/150‑1 (MQ=255)
                                                                                      ttGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg                                                                >  2:357136/1‑150 (MQ=255)
                                                                                       tGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggga                                                               >  1:844446/1‑149 (MQ=255)
                                                                                          tGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCACATATGTGCata                                                            >  1:32034/1‑150 (MQ=255)
                                                                                           gAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGCGCACATATGGGCatat                                                           >  2:29525/1‑150 (MQ=255)
                                                                                             aGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATAATGAAACGGGCGTATAACCCGCCCGTTGTTATATTTATGTGGATATTATTAATAGAATATCGAGCATAATTATATGAAGCCCATAACTAGAGACCAAATGGGCatatat                                                         >  2:12740/1‑148 (MQ=255)
                                                                                                gATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTta                                                      >  2:1676862/1‑150 (MQ=255)
                                                                                                 atcatcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtat                                                     >  1:322518/1‑150 (MQ=255)
                                                                                                  tcatcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtata                                                    >  1:238192/1‑150 (MQ=255)
                                                                                                        aaaaCATGTAATCTCTCCATGTTTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCaa                                              <  2:1092595/150‑1 (MQ=255)
                                                                                                        aaaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTAATAATGCaa                                              >  2:95517/1‑150 (MQ=255)
                                                                                                              tGTAATCTCTCCATGTGTTAAATATTGAAACGGGCTTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACtatta                                        <  2:1964758/150‑1 (MQ=255)
                                                                                                              tGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACtatta                                        <  2:1162318/150‑1 (MQ=255)
                                                                                                                tAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATg                                      >  2:1560562/1‑150 (MQ=255)
                                                                                                                    ctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCAATTACTTGAGCCCATATGGGAATATTTTTATAATGCAACTATTATGTaaa                                  >  2:1269223/1‑150 (MQ=255)
                                                                                                                         cATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAATCCCATTACTTGAGCCCATATGGGCATAATTTTATAATGCAACTATTATGTAAACattt                             >  1:228315/1‑150 (MQ=255)
                                                                                                                         cATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTAAAATGCAACTATTATGTAAACattt                             >  2:959701/1‑150 (MQ=255)
                                                                                                                            gtgtTAAATATTGCAAAGTGCGTATCACACGCGCGTTGTTTTATTTATGTGGATATTATTAAAAGCATATCGAGCATATATATATGAAGCCCATTACTTGAGCCCATAGGGGCATATTTTTATAATGCAACTATTATGTAAACatttatt                          <  2:1056910/150‑1 (MQ=255)
                                                                                                                             tgtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACCTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACatttattt                         >  2:1665171/1‑150 (MQ=255)
                                                                                                                              gtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTg                        <  1:1635739/150‑1 (MQ=255)
                                                                                                                              gtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTg                        >  2:750946/1‑150 (MQ=255)
                                                                                                                              gtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTg                        <  1:1101406/150‑1 (MQ=255)
                                                                                                                                  aaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTTTGTAAACATTTATTTGTTAt                    >  2:1622894/1‑150 (MQ=255)
                                                                                                                                   ataaaacaaCGGGCGGGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTAtt                   >  1:164247/8‑150 (MQ=255)
                                                                                                                                      ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTg                <  2:1181819/150‑1 (MQ=255)
                                                                                                                                      ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTg                <  2:866062/150‑1 (MQ=255)
                                                                                                                                      ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTg                >  2:1468602/1‑150 (MQ=255)
                                                                                                                                      ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTg                >  2:1745384/1‑150 (MQ=255)
                                                                                                                                      ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGc               <  1:306115/150‑1 (MQ=255)
                                                                                                                                      ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGACCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTg                >  1:635177/1‑150 (MQ=255)
                                                                                                                                       tGAAACGGGCGTATAACACTCCCGTTGTTTTATTTATGTGGATAATATTAATATCATATCGATCATATTTATATGTAGCCCATTACTTGTGCCCATATGGGTATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGc               <  2:202219/150‑1 (MQ=255)
                                                                                                                                       tGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGc               >  2:1995934/1‑150 (MQ=255)
                                                                                                                                        gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTTGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCt              <  1:1900115/150‑1 (MQ=255)
                                                                                                                                        gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATTGGCATATTTTTATAATGCAAATTTTATGTAAACATTTATTTGTTATTTTGCt              >  1:1347191/1‑150 (MQ=255)
                                                                                                                                        gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATATTGCAACTATTATGTAAACATTTATTTGTTATTTTGCt              >  1:1353922/1‑150 (MQ=255)
                                                                                                                                        gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCt              >  1:723852/1‑150 (MQ=255)
                                                                                                                                        acaaCGGGCGTGTTATACGCCCGATGTTTAATTTATGTGGATATTATTAATAGCATAACGAGCATATTTATATGAAGCCCATTACTAGAGCCCATATGGGCATATTTTTATAAAGCAACTATTAAGTAAACATTTATTTGTTATTTTGCt              >  1:1757116/3‑150 (MQ=255)
                                                                                                                                          aaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCACATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTTCttt            >  1:491858/1‑150 (MQ=255)
                                                                                                                                            cGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTtct          <  1:750946/150‑1 (MQ=255)
                                                                                                                                                  tataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTgga    <  1:1470063/150‑1 (MQ=255)
                                                                                                                                                    taACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTggagg  <  1:2006102/150‑1 (MQ=255)
                                                                                                                                                    taACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTggagg  <  2:153328/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGG  >  NC_000913/1211463‑1211760

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