Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,616 C→T 5.2% intergenic (‑518/+80) ppc ← / ← argE phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,153,6160CT 5.2% 231.5 / 3.5 82intergenic (‑518/+80)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (47/28);  new base T (2/2);  total (51/31)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

ACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTC  >  NC_000913/4153467‑4153764
                                                                                                                                                     |                                                                                                                                                    
aCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgc                                                                                                                                                      >  2:844517/1‑150 (MQ=39)
  tGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGACGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcat                                                                                                                                                    >  1:734112/1‑150 (MQ=255)
   gtcggtcgGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                                                                   <  2:1378295/150‑1 (MQ=39)
   gtaggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                                                                   <  2:821577/147‑1 (MQ=21)
   gtaggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                                                                   <  2:720576/147‑1 (MQ=21)
   gtaggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                                                                   <  2:720590/147‑1 (MQ=21)
                gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACatta                                                                                                                                      >  1:800932/1‑150 (MQ=37)
                 aTGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACattat                                                                                                                                     >  1:66397/1‑150 (MQ=37)
                 aTGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACattat                                                                                                                                     >  1:675345/1‑150 (MQ=37)
                  tGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTATGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCAGACaatatt                                                                                                                                    >  2:1247350/1‑150 (MQ=255)
                   gcgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTg                                                                                                                                   >  1:447381/1‑150 (MQ=37)
                    cgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCTCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGc                                                                                                                                  <  2:413474/150‑1 (MQ=37)
                    cgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGc                                                                                                                                  <  1:1308692/150‑1 (MQ=37)
                      cAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCgg                                                                                                                                >  1:421655/1‑150 (MQ=37)
                      cAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTGTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCGTCCTACATAATTGAgg                                                                                                                                >  1:421846/1‑150 (MQ=11)
                        aGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCAACGGGTGAGTCACTGGAGGTCGGAGAAGATGCGCAAGCAGCGCATCCGCCAGTATAGCGGCa                                                                                                                              >  2:722391/1‑150 (MQ=11)
                          tATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCAcg                                                                                                                            >  2:1377611/1‑149 (MQ=37)
                            tCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTgg                                                                                                                          >  1:874387/1‑150 (MQ=37)
                             cGCATCCGTCAGTATTGCGGCTCTGTAGTTGGGCTACCGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCTTTTGACTCCCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGa                                                                                                                         <  2:1345154/150‑1 (MQ=11)
                             cGCATCCGACATTATTGCGGCTCTGGAGTTTGGCAACAGTGCCGGTTGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGACTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGa                                                                                                                         <  2:1345894/150‑1 (MQ=255)
                                 tCCGACATTATTGCTGCACTGGATTTTGGCAACAGTGCCGGTTGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGttt                                                                                                                     <  1:1211016/150‑1 (MQ=255)
                                 tCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTTCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACTGTTGGGCATAGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGttt                                                                                                                     <  1:234995/150‑1 (MQ=255)
                                 tCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCAACGTTTGAGTCACTGTAGGTCGGATAAGATGAGCAAGCATCGCATCAGACATTATTGCGGCACTGGAGttt                                                                                                                     >  2:676498/1‑150 (MQ=255)
                                  ccGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACTGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTg                                                                                                                    <  2:968235/150‑1 (MQ=21)
                                  ccGACATGCTTGCGGCACTGGAGTTTGTCCACCGTGCCGGTTGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTg                                                                                                                    <  2:968236/150‑1 (MQ=11)
                                   cGACATTATTGCGGCACTGGAGTTTGTCTAGAGTGCAGGATGCGGCGCGAGCGCCTTATCCGGCCTACTGTTGGGCATCGATTGATTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTgg                                                                                                                   <  2:839027/150‑1 (MQ=11)
                                   cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGAATCGCATCCGACATTATTGCGGCACTGGAGTTTgg                                                                                                                   >  2:774880/1‑150 (MQ=21)
                                    gACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGc                                                                                                                  <  1:433309/150‑1 (MQ=37)
                                    gACATGATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGGGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGc                                                                                                                  <  2:769520/150‑1 (MQ=255)
                                      cATTATTGCTGCACTGGATTTTGGCTACTGTGCAGGATGCGGCGAGTGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGCGTCTGTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                                                                <  2:683277/150‑1 (MQ=11)
                                      cATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACTGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                                                                <  1:729568/150‑1 (MQ=21)
                                       attattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAAc                                                                                                               >  1:155140/1‑150 (MQ=37)
                                       attattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAAc                                                                                                               >  1:121330/1‑150 (MQ=37)
                                         tattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAg                                                                                                             >  1:448754/1‑150 (MQ=37)
                                             gCGGCACTGGAGTTTGGCAACAGTGCCGTATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGcc                                                                                                         <  1:358103/150‑1 (MQ=255)
                                             gCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGcc                                                                                                         >  2:1330518/1‑150 (MQ=39)
                                                     ggAGTTTGGCAACCGTGCCGGATGCGGCGCGAGCGCCTGAGCCGGCCTACGGTGGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATgcgg                                                                                                 <  2:71075/150‑1 (MQ=255)
                                                       aGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTGTGAGTCACTGTAGGTCGGAGAAGATGCGCACGCTGCGCATCCGACATTATTGGGGCACGGGAGTTTGGCAACAGCGCCGGATgcggcg                                                                                               >  2:1377082/1‑150 (MQ=14)
                                                        gTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcgca                                                                                                                                                >  1:47157/1‑100 (MQ=40)
                                                        gTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcgca                                                                                                                                                <  2:47157/100‑1 (MQ=40)
                                                          ttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCTGTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTGTGGCAACAGCGCCGGATGCGGCGCGa                                                                                            >  2:328062/1‑150 (MQ=255)
                                                             gCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAgcg                                                                                         >  1:1168253/1‑150 (MQ=39)
                                                                   gagactgATGCGGCGCGCGCCCCTCATCCGGCCTACGGGTGGGCATCGTTCGATACTCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                                                                   <  2:351404/144‑1 (MQ=14)
                                                                    tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATc                                                                                  >  2:696242/1‑150 (MQ=39)
                                                                     gCCGGATGCGGCGCGAGCGCCATATACGGCCAACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAACCAACGCATCCGACATAAAAGCGGCACAGGAGTTTGGCAAAAGCGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                 >  1:860881/1‑150 (MQ=14)
                                                                       cgcttGCGGCGCGAGCGCGAAATCCGGCCCACGGTGGGGCAGCGACCGCGCCCGTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGGGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                               <  1:509444/146‑1 (MQ=11)
                                                                         gATGCGGCGCGAGCGCCTTATCCGGCCTACTGTTGGGAATCGTTTGAGTCAGTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                             <  1:281328/150‑1 (MQ=255)
                                                                         gATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                             <  1:281299/150‑1 (MQ=39)
                                                                         gATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                             >  2:459326/1‑150 (MQ=39)
                                                                         gATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGGTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                             >  2:429482/1‑150 (MQ=255)
                                                                              ggcgCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGt                                                                        >  1:1385906/1‑150 (MQ=25)
                                                                               gcgcGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGGATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGtt                                                                       >  2:1284825/1‑150 (MQ=38)
                                                                               gcgcGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGtt                                                                       >  1:730358/1‑150 (MQ=40)
                                                                                cgcgAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGttt                                                                      >  2:36057/1‑150 (MQ=255)
                                                                                 gcgAGCGCCTTATCCGGCCTACTGTTGGGCATCGTTTGAGTCCCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGtttt                                                                     <  2:1359867/150‑1 (MQ=255)
                                                                                   gAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTaa                                                                   >  2:1148010/1‑150 (MQ=255)
                                                                                    cgcgcCTTATCCGTCCTAATCGTGGGCAGATTTAGATTCTCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAAt                                                                  <  2:1335225/149‑1 (MQ=14)
                                                                                      cgcCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGc                                                                >  1:1031826/1‑150 (MQ=255)
                                                                                       gcCTTATCCGGCCTACGGTTGGACATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGcc                                                               >  1:853340/1‑150 (MQ=255)
                                                                                        ccTTATCCGGCCTACGGGTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCa                                                              <  1:716314/150‑1 (MQ=255)
                                                                                         cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAg                                                             >  2:35299/1‑150 (MQ=255)
                                                                                         cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTAATCCGGCCTACGTTTTAATGCCAg                                                             >  2:1294507/1‑150 (MQ=255)
                                                                                          ttATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGc                                                            >  1:871631/1‑150 (MQ=255)
                                                                                                  ccatcGGGTGGGAATAGTGTGAGTCTCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCTCTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATgg                                                    <  1:819379/146‑1 (MQ=17)
                                                                                                   cTACTGTTGGGCATCGTTTGCGTCGGTGTAGGTCGGATAAGATGCGCAAGGATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGt                                                   <  2:990261/150‑1 (MQ=17)
                                                                                                        gTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAtt                                              >  2:765229/1‑150 (MQ=255)
                                                                                                           gggCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTAcc                                           >  1:775234/1‑150 (MQ=255)
                                                                                                             gCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAAATGTGAATTACCTg                                         >  1:531088/1‑150 (MQ=255)
                                                                                                              cATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTgg                                        >  1:781099/1‑150 (MQ=255)
                                                                                                                   tttGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTAt                                   <  1:1168102/150‑1 (MQ=255)
                                                                                                                       aGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGt                               >  2:1149464/1‑150 (MQ=255)
                                                                                                                             tGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCgg                         >  1:1195820/1‑150 (MQ=255)
                                                                                                                             tGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCgg                         >  1:1210745/1‑150 (MQ=255)
                                                                                                                                 ggTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTgg                     <  1:1279969/150‑1 (MQ=255)
                                                                                                                                  gTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCGAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTggg                    >  1:518194/1‑150 (MQ=255)
                                                                                                                                     ggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCtt                 <  2:682952/150‑1 (MQ=255)
                                                                                                                                           gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTggg                    >  2:606477/1‑141 (MQ=255)
                                                                                                                                           gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTggg                    <  1:606477/141‑1 (MQ=255)
                                                                                                                                            aTGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAAc          >  1:1092045/1‑150 (MQ=255)
                                                                                                                                            aTGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAAc          >  1:1092044/1‑150 (MQ=255)
                                                                                                                                            aTGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAAc          >  1:1092042/1‑150 (MQ=255)
                                                                                                                                              gcgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCg        >  1:549659/1‑150 (MQ=255)
                                                                                                                                               cgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt       >  1:386995/1‑150 (MQ=255)
                                                                                                                                               cgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt       <  1:849219/150‑1 (MQ=255)
                                                                                                                                               cgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt       >  2:501556/1‑150 (MQ=255)
                                                                                                                                               cgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt       >  2:501553/1‑150 (MQ=255)
                                                                                                                                               cgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt       <  1:288541/150‑1 (MQ=255)
                                                                                                                                               cgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt       >  1:195508/1‑150 (MQ=255)
                                                                                                                                               cgcAAGCATCGCATCCGACAGTCTTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt       <  1:537456/150‑1 (MQ=255)
                                                                                                                                                gcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgtg      >  1:1167887/1‑150 (MQ=255)
                                                                                                                                                    gcatcgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGGTTc  >  1:1067020/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTC  >  NC_000913/4153467‑4153764

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: