Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,795,228 G→C 9.5% 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 9.5% 366.0 / 25.5 116intergenic (‑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 (47/58);  new base C (5/6);  total (52/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

Only 100 of 118 total aligned reads displayed.

GCAGCATCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCC  >  NC_000913/1795080‑1795369
                                                                                                                                                    |                                                                                                                                             
gcagcaTCAAAACTGACAGGCATAATAACCAGCGTATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCATCGACAGTGAAAAGAAAAAAGGCCGCaga                                                                                                                                              >  2:142268/1‑150 (MQ=255)
 cagcgTCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCCGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCagcg                                                                                                                                             >  1:730830/1‑148 (MQ=255)
  agcaTCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAATAGAAAAAAGGCCGCAGAGc                                                                                                                                            >  1:268731/1‑150 (MQ=255)
    caTCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCgg                                                                                                                                          <  1:1565924/150‑1 (MQ=255)
       cAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTAATGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCt                                                                                                                                       >  2:1177271/1‑150 (MQ=255)
         aaaCTGACAGGCATAATACCCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGAATCCGTTCTGCTGATGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttt                                                                                                                                     <  2:650313/150‑1 (MQ=255)
          aaCTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCtttt                                                                                                                                    >  2:1080666/1‑150 (MQ=255)
           aCTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAGAAGGCCGCAGAGCGGCCttttt                                                                                                                                   <  2:1097990/150‑1 (MQ=255)
           aCTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttttt                                                                                                                                   <  2:355671/150‑1 (MQ=255)
           aCTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCttttt                                                                                                                                   <  2:3984/150‑1 (MQ=255)
            cTGACAGGCATAAAAACCAGAGAATATTTTGTCGCTGTTGTTGGCGGGAAAGTGTAAGAATGGTATCCGTTCTGCTGAAGGGTCATGGCGGTGATTTTACTGTGACTCTTCGCCAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTtta                                                                                                                                  <  2:1095173/150‑1 (MQ=255)
             tGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttttttc                                                                                                                                 <  1:13233/150‑3 (MQ=255)
             tGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttag                                                                                                                                 <  1:11007/150‑1 (MQ=255)
             tGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttag                                                                                                                                 >  1:737790/1‑150 (MQ=255)
             tGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTGCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttttttc                                                                                                                                 <  2:929384/150‑3 (MQ=255)
              gACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttttttct                                                                                                                                <  1:1805108/150‑4 (MQ=255)
                cAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagtta                                                                                                                              <  2:39964/150‑1 (MQ=255)
                cAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGAGACTCTTCGACAGTTAAAAGAAAAAAGGCCGAAGAGCGGCCTTTttagtta                                                                                                                              >  2:252716/1‑150 (MQ=255)
                cAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGATGGCGGGCAAGTGTCAGCATGGAATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGTTGACTCTTCGACAGTGAAAAGAAAAAAGTACGCAGACCGGCCTTTttagtta                                                                                                                              >  1:1674692/1‑150 (MQ=255)
                 aGGCATAATACCCCGCGAATATTATGTCGCTGTTGTAGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagttag                                                                                                                             <  1:347671/150‑1 (MQ=255)
                 aGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagttag                                                                                                                             >  2:1531892/1‑150 (MQ=255)
                  ggCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttttttcttttc                                                                                                                            <  2:219803/150‑8 (MQ=255)
                  ggCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGa                                                                                                                            >  1:1402258/1‑150 (MQ=255)
                  ggCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGa                                                                                                                            >  1:365887/1‑150 (MQ=255)
                   gCATAATTACGACCGAATATTTTGTCGCTGTTGTTGGAGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTTATTATACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGAt                                                                                                                           <  2:1633580/150‑1 (MQ=255)
                    cATAATAACCAGCGATTATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAACGGCCGCAGAGCGGCCTTTTTAGTTAGATc                                                                                                                          <  1:445074/150‑1 (MQ=255)
                    cATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATc                                                                                                                          <  2:278247/150‑1 (MQ=255)
                    cATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATc                                                                                                                          <  1:1443249/150‑1 (MQ=255)
                    cATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATc                                                                                                                          >  2:149787/1‑150 (MQ=255)
                     ataataACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGACATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCa                                                                                                                         <  1:80796/150‑1 (MQ=255)
                      taataaGCAGCTTATATTTTGTTGCTGTTGTTGTCGGGTAATTTTTATCCAGGTATCCCTTCTTCTGAAGAGTCACGGCGTTGTATTTACGGGGACTCATCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAg                                                                                                                        <  1:110146/150‑1 (MQ=255)
                        ataaCCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGAt                                                                                                                      <  2:483668/150‑1 (MQ=255)
                        ataaCCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGAt                                                                                                                      <  1:905295/150‑1 (MQ=255)
                           aCCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTTAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTAc                                                                                                                   <  1:1420931/150‑1 (MQ=255)
                              aGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCg                                                                                                                >  2:523261/1‑150 (MQ=255)
                               gCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGt                                                                                                               >  1:42308/1‑150 (MQ=255)
                                      ttttGGCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTggg                                                                                                        <  1:739947/150‑1 (MQ=255)
                                        ttGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCg                                                                                                      <  2:910762/150‑1 (MQ=255)
                                        ttGTCGCTGGTGTCGGCGCGCAAGTGTCCGCACGGCATCCCTGCTGCCGAAGCGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGACACGAAAAAAGGCCGAAGAGCGGCCTTTTTAGTTAGCTCAGATTACTCGTCTTTGGGCg                                                                                                      <  2:267948/150‑1 (MQ=255)
                                             gCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGGCATGGCGTTGAGTTTACGGGGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCg                                                                                                 <  1:251322/150‑1 (MQ=255)
                                                  tgttgGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTc                                                                                            >  1:1582159/1‑150 (MQ=255)
                                                  tgttgGCGGGCAAGTGTCAGCATGGTATCCGTTATGCTGAAGTGTCATGTCGTTGATTTTACGGAGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGCCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTc                                                                                            >  2:1591287/1‑150 (MQ=255)
                                                   gttgGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCg                                                                                           <  1:1569866/150‑1 (MQ=255)
                                                   gttgGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCg                                                                                           <  1:578878/150‑1 (MQ=255)
                                                   gttgGAGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGCCATGGCGTAGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTGTTTTGTTGGATCAGATTACTCGTGTTTGGGCGAAGCGTTTTCg                                                                                           >  2:21253/1‑150 (MQ=255)
                                                    ttgGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCGTTTTTGTTAGATCAGATTACTCGTCTTTGGGCGAAGAGTTTTCGa                                                                                          >  1:1746955/1‑150 (MQ=255)
                                                      ggCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGAcc                                                                                        >  1:954674/1‑150 (MQ=255)
                                                         gggCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCgg                                                                                     <  1:609314/150‑1 (MQ=255)
                                                            cAAGTGTCAGCATGTTATCCGTTCTGCTGAAGTGTCTTGGCGTTGATTTTACGGTTACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCtt                                                                                  <  2:1609774/150‑1 (MQ=255)
                                                                tgtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTa                                                                              <  2:1556482/150‑1 (MQ=255)
                                                                 gtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTAGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTaa                                                                             <  1:905368/150‑1 (MQ=255)
                                                                    aGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACtt                                                                          >  1:775921/1‑150 (MQ=255)
                                                                      cATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCt                                                                        >  1:1318701/1‑150 (MQ=255)
                                                                      atgggTTTCCGTTCTGGGTATGTGTCTTGGCGTTGATTTTACGGTTACTCTTCTACAGTGAAAATAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCt                                                                        <  2:754348/147‑1 (MQ=255)
                                                                       aTGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTg                                                                       <  1:1718257/150‑1 (MQ=255)
                                                                       aTGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTg                                                                       <  1:839755/150‑1 (MQ=255)
                                                                            aTCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCgg                                                                  >  2:1792022/1‑150 (MQ=255)
                                                                             tCCCTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCggg                                                                 <  2:774970/150‑1 (MQ=255)
                                                                                    tgctgAAGTGTCATGGCGTTGATTTTCCGGTGACTCTTCGACAGTGAAAACAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTGTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAg                                                          <  1:345751/150‑1 (MQ=255)
                                                                                    tgctgAAGTGTCATGGCGTTGATTTTAGGGTGACTCTTCGACAGTGAAATGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAg                                                          <  1:120850/150‑1 (MQ=255)
                                                                                     gctgAAGTGTCATGGCGTTGATTTTCCGGGGACTCTTCGACAGTGAAAACAAAAACGCCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGGCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAgg                                                         <  2:146281/150‑1 (MQ=255)
                                                                                      ctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGTGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGt                                                        >  2:183370/1‑150 (MQ=255)
                                                                                             gtCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCAcg                                                 >  1:1164747/1‑150 (MQ=255)
                                                                                              tCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCAcgc                                                <  1:602401/150‑1 (MQ=255)
                                                                                                       tGATTTTACGGGGACTCATCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTg                                       <  1:1495963/150‑1 (MQ=255)
                                                                                                           tttACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGGGCTGTAAt                                   >  1:641776/1‑150 (MQ=255)
                                                                                                            ttACGGGGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGGCACCACGCGCCGGGCTGTAATg                                  <  2:1113908/150‑1 (MQ=255)
                                                                                                             tACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATgg                                 <  1:1152681/150‑1 (MQ=255)
                                                                                                              aCGGCGACTCTTCGACAGTGAAAAGAAAAAAGCCCCCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCGGTAATggg                                <  2:270498/150‑1 (MQ=255)
                                                                                                               cGGTGCCTCTGCGACAGTGAACAGACAAACGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTCACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGa                               <  1:259960/150‑1 (MQ=255)
                                                                                                               cGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGa                               >  2:420497/1‑150 (MQ=255)
                                                                                                                ggTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGaa                              <  2:1164747/150‑1 (MQ=255)
                                                                                                                ggTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGaa                              <  2:594787/150‑1 (MQ=255)
                                                                                                                ggTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCATGCTGTAATGGGaa                              <  2:1470777/150‑1 (MQ=255)
                                                                                                                   gACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAatat                           >  2:737023/1‑150 (MQ=255)
                                                                                                                     ctctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGTTCACCACGCGCCGTGATGTAATGGGAATATcc                         >  1:143133/1‑150 (MQ=255)
                                                                                                                       ctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTc                       >  2:6619/1‑150 (MQ=255)
                                                                                                                          cGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGcc                    <  2:365887/150‑1 (MQ=255)
                                                                                                                                 gaaaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTc             >  2:594124/1‑150 (MQ=255)
                                                                                                                                   aaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCgg           >  2:67244/1‑150 (MQ=255)
                                                                                                                                   aaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCgg           >  2:188978/1‑150 (MQ=255)
                                                                                                                                   aaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGACGTGCTGTAATGGGAATATCCTCTCCCGTTTTCgg           >  2:287563/1‑150 (MQ=255)
                                                                                                                                   aaagaaaaAACGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTACgg           >  2:551286/1‑150 (MQ=255)
                                                                                                                                    aactaaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACGCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCggg          >  1:265666/5‑150 (MQ=255)
                                                                                                                                     agaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCCCGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGt         >  2:1696029/1‑150 (MQ=255)
                                                                                                                                     agaaaaAAGGCCGCAGAGCGGCCGTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCAGGCTTTTTAACTTCTGACCGGGTCGGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGt         >  2:747811/1‑150 (MQ=255)
                                                                                                                                     actaaaaaGGCCGCTCTGCGGGCTTGTTAGTTAGATCAGATTACTCGGCGTTGGGCGCAGCGGTTTCGACCCGGCGATGTAACTTCTGCCCGGGTCGGAAGGTCACCACGCGCCGTGCGGTAATGCGTATATACTCGCCCGGTTGCgggg         >  1:925338/4‑149 (MQ=255)
                                                                                                                                     actaaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGt         >  2:1657136/4‑150 (MQ=255)
                                                                                                                                      gaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGtt        <  2:42308/150‑1 (MQ=255)
                                                                                                                                      gaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGtt        >  2:383468/1‑150 (MQ=255)
                                                                                                                                      gaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGtt        >  1:1055345/1‑150 (MQ=255)
                                                                                                                                       aaaaaaGGCCGCAGCGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGATTATCCTCGCCCGTTTTCGGGTTa       >  1:533716/1‑150 (MQ=255)
                                                                                                                                       aaaaaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGGGCCGTGCTGTAATGGGAATATGCTCGCCCGTTTTCGGGGTa       >  1:1383819/1‑150 (MQ=255)
                                                                                                                                       aaaaaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTa       >  2:1775500/1‑150 (MQ=255)
                                                                                                                                        aaaaaGGCCGCTCTGGGTCATATTTAGTTAGATCAGATTACTCGTCTTTGGGCGCAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGACTCAAGGTCACCACGCCGCGTGCTGGAATGCGAGTATCCTCGCACTTTTTCGGGTTAc      >  2:110680/1‑150 (MQ=255)
                                                                                                                                        aaaaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGCATCGTTTTCGACCCGGCTTGTTAACTTCTGCCCGGGTGTGAAGGTCACCACGCGCCGTGCGGTATTGGGTATCTCGTCGCCGGTGTTCGGGttgc      >  1:299641/1‑148 (MQ=255)
                                                                                                                                        aaaaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTAc      >  1:1524525/1‑150 (MQ=255)
                                                                                                                                         aaaaGCCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGGTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACg     <  2:308079/150‑1 (MQ=255)
                                                                                                                                          aaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGACCGGGTCTGAAGGTCACCACTCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGt    >  1:314209/1‑150 (MQ=255)
                                                                                                                                            aGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCACGTTTTCGGGTTACGTcc  >  1:1253063/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GCAGCATCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCC  >  NC_000913/1795080‑1795369

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