Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,795,228 G→C 8.7% intergenic (‑76/+25) btuC ← / ← ihfA vitamin B12 ABC transporter permease/integration host factor (IHF), DNA‑binding protein, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,795,2280GC 8.7% 286.3 / 13.9 90intergenic (‑76/+25)btuC/ihfAvitamin B12 ABC transporter permease/integration host factor (IHF), DNA‑binding protein, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (45/37);  new base C (4/4);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

AGCATCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTT  >  NC_000913/1795082‑1795377
                                                                                                                                                  |                                                                                                                                                     
agcaTCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGATGTTGGCGGGCAAGTGTCAGCATCGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGc                                                                                                                                                    <  1:284334/150‑1 (MQ=255)
 gcaTCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCg                                                                                                                                                   >  1:970069/1‑150 (MQ=255)
 gcaTCAAAACTGACAGGCATAATAACAAGAGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATTGAATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCagagac                                                                                                                                                   >  1:1463602/1‑148 (MQ=255)
  caTCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCgg                                                                                                                                                  <  2:544927/150‑1 (MQ=255)
  caTCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCgg                                                                                                                                                  >  1:1641590/1‑150 (MQ=255)
   aTCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCATGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGc                                                                                                                                                 <  1:1046269/150‑1 (MQ=255)
   aTCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGc                                                                                                                                                 >  2:72095/1‑150 (MQ=255)
    tCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGcc                                                                                                                                                <  1:1492018/150‑1 (MQ=255)
    tCAAAACTGACAGGCATAATAACCAGCGAATATTTTCTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGATCGGcc                                                                                                                                                >  2:116769/1‑150 (MQ=255)
     cAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCt                                                                                                                                               >  1:538670/1‑150 (MQ=255)
     cAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCt                                                                                                                                               >  2:928667/1‑150 (MQ=255)
     cAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCt                                                                                                                                               >  2:1458545/1‑150 (MQ=255)
       aaaCTGACAGGCATACTACCCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCttt                                                                                                                                             <  2:207041/150‑1 (MQ=255)
       aaaCTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttt                                                                                                                                             <  2:1002430/150‑1 (MQ=255)
       aaaCTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttt                                                                                                                                             <  2:965272/150‑1 (MQ=255)
       aaaCTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCttt                                                                                                                                             >  2:165270/1‑150 (MQ=255)
           tGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttttttc                                                                                                                                         <  1:1489932/150‑3 (MQ=255)
           tGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttttttc                                                                                                                                         <  2:1151410/150‑3 (MQ=255)
                  cATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATc                                                                                                                                  >  1:189144/1‑150 (MQ=255)
                    taataaCCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAg                                                                                                                                >  2:501561/1‑150 (MQ=255)
                          ccAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACt                                                                                                                          <  2:1495932/150‑1 (MQ=255)
                            aGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCg                                                                                                                        >  2:973499/1‑150 (MQ=255)
                            aGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCg                                                                                                                        >  2:559930/1‑150 (MQ=255)
                             gCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGt                                                                                                                       >  2:213459/1‑150 (MQ=255)
                              cGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTc                                                                                                                      <  2:49953/150‑1 (MQ=255)
                                  tatTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTg                                                                                                                  <  1:786858/150‑1 (MQ=255)
                                  tatTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTATTTAGTTAGATCAGATTACTCGTCTTTg                                                                                                                  >  1:904446/1‑150 (MQ=255)
                                   atTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTgg                                                                                                                 <  2:1458079/150‑1 (MQ=255)
                                    ttttGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTTTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTggg                                                                                                                <  1:855422/150‑1 (MQ=255)
                                    ttttGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTggg                                                                                                                <  1:1458545/150‑1 (MQ=255)
                                     tttGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGc                                                                                                               >  1:549652/1‑150 (MQ=255)
                                     tttGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACGCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGc                                                                                                               <  1:949661/150‑1 (MQ=255)
                                       tGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGa                                                                                                             >  1:1038892/1‑150 (MQ=255)
                                         tCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAg                                                                                                           <  2:970069/150‑1 (MQ=255)
                                           gCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCg                                                                                                         >  1:405512/1‑150 (MQ=255)
                                            cTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGt                                                                                                        <  2:574092/150‑1 (MQ=255)
                                                       gggCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCgg                                                                                             <  1:525033/150‑1 (MQ=255)
                                                       gggCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCgg                                                                                             <  1:1599954/150‑1 (MQ=255)
                                                               gtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTaa                                                                                     >  1:572735/1‑150 (MQ=255)
                                                               gtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTaa                                                                                     >  1:902654/1‑150 (MQ=255)
                                                                 cAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGAGTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACt                                                                                   <  1:937486/150‑1 (MQ=255)
                                                                  aGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACtt                                                                                  >  2:1567783/1‑150 (MQ=255)
                                                                      tGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGc                                                                              >  2:1700175/1‑150 (MQ=255)
                                                                       ggTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGcc                                                                             >  1:1585684/1‑150 (MQ=255)
                                                                        gTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGccc                                                                            >  2:586545/1‑150 (MQ=255)
                                                                          aTCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCgg                                                                          >  1:1268585/1‑150 (MQ=255)
                                                                           tCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGAATCTTCGACAGTGAAAAGAAAAAAGGCCGCCGAGCGGCCGTTTTAGTTAGATCAGATTACTCGTCTTTGGTCGAAGCTTTTTCGACCCGGCTTTTTAACTTCTGCCCggg                                                                         <  1:116769/150‑1 (MQ=255)
                                                                             cGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGATAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTACGACCCGGCTTTTTAACTTCTGCCCGGGTc                                                                       >  1:1540745/1‑150 (MQ=255)
                                                                             cGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGACGCAGAGCGTACTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTc                                                                       >  2:1319154/1‑150 (MQ=255)
                                                                                  tgctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGAAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAg                                                                  <  2:97678/150‑1 (MQ=255)
                                                                                    ctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCggg                                                                         >  1:124017/1‑141 (MQ=255)
                                                                                    ctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCggg                                                                         <  2:124017/141‑1 (MQ=255)
                                                                                        aGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTcacc                                                            <  1:72095/150‑1 (MQ=255)
                                                                                         gtgtCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTcacca                                                           >  1:1306829/1‑150 (MQ=255)
                                                                                                   gTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGTGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGc                                                 <  1:501561/150‑1 (MQ=255)
                                                                                                     tGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTg                                               >  1:1633241/1‑150 (MQ=255)
                                                                                                     tGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTg                                               >  1:980236/1‑150 (MQ=255)
                                                                                                      gATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGCGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTATCTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGt                                              >  2:124121/1‑150 (MQ=255)
                                                                                                            aCGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATggg                                        <  1:973499/150‑1 (MQ=255)
                                                                                                            aCGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATggg                                        <  2:1038892/150‑1 (MQ=255)
                                                                                                            aCGGTGACTCATCGACAGTGAAAAGAAAAAAGGACGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGGCACCACGCGCCGTGCTGTAATggg                                        <  1:1706436/150‑1 (MQ=255)
                                                                                                              ggTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGaa                                      >  2:1587188/1‑150 (MQ=255)
                                                                                                                  aCTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATc                                  >  1:1064567/1‑150 (MQ=255)
                                                                                                                    tctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGACGCGTTTTCGACCCGGGTTTTTAACTTCAGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCt                                >  2:81302/1‑150 (MQ=255)
                                                                                                                     ctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTc                               >  1:277237/1‑150 (MQ=255)
                                                                                                                      ttCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCg                              <  2:1269308/150‑1 (MQ=255)
                                                                                                                      ttCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCg                              >  2:516307/1‑150 (MQ=255)
                                                                                                                              tGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGtttt                      >  1:1073155/1‑150 (MQ=255)
                                                                                                                               gaaaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGtttta                     >  1:1433976/1‑149 (MQ=255)
                                                                                                                                 taactaaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCgg                   >  2:1094547/6‑150 (MQ=255)
                                                                                                                                 aaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCgg                   >  1:870792/1‑150 (MQ=255)
                                                                                                                                 aaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTAGCCCGTTTTCgg                   >  2:1296276/1‑150 (MQ=255)
                                                                                                                                  aactaaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCggg                  >  1:831283/5‑150 (MQ=255)
                                                                                                                                   agaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGt                 <  1:616246/150‑1 (MQ=255)
                                                                                                                                   agaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGt                 >  2:1317250/1‑150 (MQ=255)
                                                                                                                                    gaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGtt                >  1:562240/1‑150 (MQ=255)
                                                                                                                                     taaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGCTTACTAGTCTTTGGGCGAAGCGTTTTCGCCCCGGCTTTTTAACTTCTGCCCGGGTCTGAGGGTCACCACGCGCAGTGCTGTAATGGGAATATACTCGAGAGTTTTCGGGTTa               >  1:1260129/2‑150 (MQ=255)
                                                                                                                                     aaaaaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTa               <  1:928667/150‑1 (MQ=255)
                                                                                                                                      aaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTAc              >  1:53791/1‑150 (MQ=255)
                                                                                                                                      aaaaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTAc              <  2:1300834/150‑1 (MQ=255)
                                                                                                                                       aaaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACg             <  2:1064567/150‑1 (MQ=255)
                                                                                                                                        aaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGt            <  1:1700175/150‑1 (MQ=255)
                                                                                                                                         aaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTc           <  1:545102/150‑1 (MQ=255)
                                                                                                                                          aGGCCGCTGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGTGCGAAGCGTTTTCTACCCTGCTTTTTAACTTCTGCCCGGGTCTGAAGGGCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGGTACGTcc          <  2:1463602/150‑1 (MQ=255)
                                                                                                                                            gCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAGTATCCTCGCCCGTTTTCGGGTTACGTCCCg        <  2:431235/150‑1 (MQ=255)
                                                                                                                                            gCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCg        >  1:1421804/1‑150 (MQ=255)
                                                                                                                                            gCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCg        <  2:1306829/150‑1 (MQ=255)
                                                                                                                                             ccGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCgg       >  1:3650/1‑150 (MQ=255)
                                                                                                                                             ccGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCgg       >  2:975459/1‑150 (MQ=255)
                                                                                                                                              cGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCggg      <  2:980236/150‑1 (MQ=255)
                                                                                                                                              cGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGGTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCggg      <  2:1315137/150‑1 (MQ=255)
                                                                                                                                                  gagCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGtt  <  2:1585684/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGCATCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTT  >  NC_000913/1795082‑1795377

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: