Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,795,229 A→T 9.4% intergenic (‑77/+24) 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,2290AT 9.4% 292.3 / 17.1 116intergenic (‑77/+24)btuC/ihfAvitamin B12 ABC transporter permease/integration host factor (IHF), DNA‑binding protein, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (46/58);  new base T (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 = 6.74e-01

Only 100 of 118 total aligned reads displayed.

GCAGCATCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCG  >  NC_000913/1795080‑1795375
                                                                                                                                                     |                                                                                                                                                  
  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)
             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:365887/1‑150 (MQ=255)
                  ggCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGa                                                                                                                                  >  1:1402258/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                                                                                                                                >  2:149787/1‑150 (MQ=255)
                    cATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATc                                                                                                                                <  1:1443249/150‑1 (MQ=255)
                     ataataACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGACATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCa                                                                                                                               <  1:80796/150‑1 (MQ=255)
                      taataaGCAGCTTATATTTTGTTGCTGTTGTTGTCGGGTAATTTTTATCCAGGTATCCCTTCTTCTGAAGAGTCACGGCGTTGTATTTACGGGGACTCATCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAg                                                                                                                              <  1:110146/150‑1 (MQ=255)
                        ataaCCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGAt                                                                                                                            <  2:483668/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)
                                     atTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTgg                                                                                                               >  2:1234227/1‑150 (MQ=255)
                                      ttttGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTTATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTggg                                                                                                              <  2:1081696/150‑1 (MQ=255)
                                      ttttGGCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTggg                                                                                                              <  1:739947/150‑1 (MQ=255)
                                        ttGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCg                                                                                                            <  2:910762/150‑1 (MQ=255)
                                        ttGTCGCTGGTGTCGGCGCGCAAGTGTCCGCACGGCATCCCTGCTGCCGAAGCGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGACACGAAAAAAGGCCGAAGAGCGGCCTTTTTAGTTAGCTCAGATTACTCGTCTTTGGGCg                                                                                                            <  2:267948/150‑1 (MQ=255)
                                         tGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCGGATTACTCGTCTTTGGGCGa                                                                                                           <  1:874739/150‑1 (MQ=255)
                                             gCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGGCATGGCGTTGAGTTTACGGGGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCg                                                                                                       <  1:251322/150‑1 (MQ=255)
                                                  tgttgGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTc                                                                                                  >  1:1582159/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:839755/150‑1 (MQ=255)
                                                                       aTGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTg                                                                             <  1:1718257/150‑1 (MQ=255)
                                                                            aTCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCgg                                                                        >  2:1792022/1‑150 (MQ=255)
                                                                            aTCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCgg                                                                        >  2:1699717/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:594787/150‑1 (MQ=255)
                                                                                                                ggTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGaa                                    <  2:1164747/150‑1 (MQ=255)
                                                                                                                ggTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCATGCTGTAATGGGaa                                    <  2:1470777/150‑1 (MQ=255)
                                                                                                                 gTGACTCATCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGTCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGGTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAat                                   <  2:1582507/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              >  1:1055345/1‑150 (MQ=255)
                                                                                                                                      gaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGtt              >  2:383468/1‑150 (MQ=255)
                                                                                                                                      gaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGtt              <  2:42308/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                  
GCAGCATCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCG  >  NC_000913/1795080‑1795375

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: