Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,570 G→A 6.7% intergenic (‑472/+126) 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,5700GA 6.7% 215.6 / 4.8 76intergenic (‑472/+126)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (43/26);  new base A (2/3);  total (46/30)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-01

GGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGT  >  NC_000913/4153422‑4153715
                                                                                                                                                    |                                                                                                                                                 
ggcgCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACg                                                                                                                                                   >  5:382455/1‑149 (MQ=39)
        cgcCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCa                                                                                                                                           >  1:240934/1‑149 (MQ=255)
              ccgatctcTACAGTTGGGCATCGTTTGAGTCACTGTCGGGCGGATAAGATGCGCAAGTATCGCATCCGAAATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                     <  1:175563/142‑1 (MQ=255)
              aTCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGt                                                                                                                                       <  8:74404/147‑1 (MQ=255)
                  ggCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGt                                                                                                                                 <  2:240934/149‑1 (MQ=255)
                      tACAGTTGGGCATCGTTTGAGTCACGATCAGACGGATAAGATGCGCAAGTATCGCATCCGAAATTATTGCGGCACTGGAGTTTGGCAAAATTGACGGATGCGGCGCGAGCGCCGTATCCGGCCTACGGTTGGGCAGCGTTTGAgctcg                                                                                                                              >  2:175563/1‑144 (MQ=17)
                                 aTCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGt                                                                                                                           <  8:282008/140‑1 (MQ=255)
                                 aTCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGt                                                                                                                           >  7:282008/1‑140 (MQ=255)
                                            cgcTGTCGGTCGGATCAGATGCGCAAGGATCGCATCCGACATTTTTGCTGCACTGGAGTGTTGCAAAAGTTCAGGATGCGGCGCGAGCGCCTTATCCGTCCTACTGTTGGGAATATTTTGAGTCACTGTAGGTCGGATAAGATGCGCaa                                                                                                       <  2:265701/147‑1 (MQ=11)
                                                   ggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCGCTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcg                                                                                                 >  7:436355/1‑148 (MQ=37)
                                                   ggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTg                                                                                                                                               <  6:426315/102‑1 (MQ=32)
                                                   ggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTg                                                                                                                                               >  5:426315/1‑102 (MQ=32)
                                                   ggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                     <  8:94405/112‑1 (MQ=32)
                                                   ggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                     >  7:94405/1‑112 (MQ=32)
                                                    gtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                                                                                                              <  1:36268/102‑1 (MQ=32)
                                                    gtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgc                                                                                                     >  5:11355/1‑143 (MQ=35)
                                                    gtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgc                                                                                                     <  6:11355/143‑1 (MQ=35)
                                                      cggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcgcatc                                                                                             >  3:43885/1‑149 (MQ=37)
                                                         aTAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgc                                                                                                     <  8:237106/138‑1 (MQ=25)
                                                         aTAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgc                                                                                                     >  7:237106/1‑138 (MQ=25)
                                                             gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGa                                                                                                                                   >  2:350329/1‑104 (MQ=32)
                                                             gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGa                                                                                                                                   <  1:350329/104‑1 (MQ=32)
                                                                  gcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTc                                                                                                                      <  4:35730/112‑1 (MQ=35)
                                                                  gcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACCGTTGGGCATCGTTTGAGTCACTGTAGGTc                                                                                                                      >  3:35730/1‑112 (MQ=17)
                                                                   cAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGATTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCg                                                                                <  8:436355/149‑1 (MQ=37)
                                                                               ccGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTGGCATCGTTTTAGACACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCTACATAATAGCGGCACAGGAGttt                                                                    >  2:458537/1‑149 (MQ=11)
                                                                                cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                  >  6:326682/1‑118 (MQ=39)
                                                                                cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                  <  5:326682/118‑1 (MQ=39)
                                                                                cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCAc                                                                            >  5:111965/1‑140 (MQ=37)
                                                                                cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCAc                                                                            <  6:111965/140‑1 (MQ=37)
                                                                                 gACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGAGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATCCGCAAGCATAGCATCCGACATTATTGCTGCACTGGAGATTgg                                                                  >  3:194103/1‑149 (MQ=11)
                                                                                    attattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                     >  1:280410/1‑79 (MQ=34)
                                                                                    attattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                     <  2:280410/79‑1 (MQ=34)
                                                                                        ttGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                  <  3:402646/110‑1 (MQ=40)
                                                                                        ttGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                  >  4:402646/1‑110 (MQ=40)
                                                                                          gCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAgcgc                                                         >  4:154055/1‑149 (MQ=39)
                                                                                              cACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGa                                                     >  3:413780/1‑149 (MQ=39)
                                                                                                cTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGATCGCCTTATCCTGCCTAAGGTTCGACATCGTATGAGACACTATAGCTCGGATAACAACCTCAAAAATCAAATCCTACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATg                                                   >  4:301423/1‑149 (MQ=11)
                                                                                                 tGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATgc                                                  >  5:162031/1‑149 (MQ=39)
                                                                                                   gAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                  >  6:34462/1‑99 (MQ=40)
                                                                                                   gAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                  <  5:34462/99‑1 (MQ=40)
                                                                                                    tgTTTGGCATCAGTGCCGGATGCGGCGCGAGCGCCTTAGCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGTATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATgcggc                                               <  3:301423/148‑1 (MQ=255)
                                                                                                      tttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGgcgc                                             >  1:391352/1‑149 (MQ=39)
                                                                                                            aaCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATcc                                                                                            <  4:234025/96‑1 (MQ=18)
                                                                                                             aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGcc                                      >  1:343105/1‑149 (MQ=39)
                                                                                                              cAGTGCCGGATGCGGCGCGTGCGCCTTATCCGGCCTACGGTTGGGAATCGTTTGAGTCACTGTAGGTCGGATAAg                                                                                                               <  5:390946/75‑1 (MQ=14)
                                                                                                              cAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAg                                                                                                               >  4:340609/1‑75 (MQ=39)
                                                                                                              cAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAg                                                                                                               <  3:340609/75‑1 (MQ=39)
                                                                                                              cAGAGCCGGAAGCGGCGCGAGCGCCTTATCAGGCATACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAg                                                                                                               >  6:390946/1‑75 (MQ=14)
                                                                                                                 tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTGGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                  >  1:249420/1‑149 (MQ=37)
                                                                                                                      gataaggcgCTCGCGCCGCATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGc                             >  8:473885/6‑149 (MQ=14)
                                                                                                                           ggcgCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACg                        >  5:194493/1‑149 (MQ=39)
                                                                                                                             cgcgAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGtt                      >  1:8856/1‑149 (MQ=40)
                                                                                                                              gcgAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGttt                     <  4:43885/149‑1 (MQ=255)
                                                                                                                                gAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTa                   <  4:413780/149‑1 (MQ=255)
                                                                                                                                 aGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCCCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGt                       <  5:170969/144‑1 (MQ=255)
                                                                                                                                 aGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                  >  4:264781/1‑133 (MQ=40)
                                                                                                                                 aGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                  <  3:264781/133‑1 (MQ=40)
                                                                                                                                 aGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGAGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGt                       >  6:170969/1‑144 (MQ=255)
                                                                                                                                   cgcCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGAGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATg                >  8:12632/1‑149 (MQ=255)
                                                                                                                                    gcCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGTCGCGAGCGCCTTATCCGTCGTCCGTTTTAATga               >  3:253256/1‑148 (MQ=255)
                                                                                                                                     ccTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACt                                                                                                                             <  5:51146/38‑1 (MQ=40)
                                                                                                                                     ccTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACt                                                                                                                             >  6:51146/1‑38 (MQ=40)
                                                                                                                                     ccTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGAc                                                                                         <  8:356056/74‑1 (MQ=21)
                                                                                                                                     ccTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGAc                                                                                         >  7:356056/1‑74 (MQ=21)
                                                                                                                                     ccTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACa                                                                                        <  1:57654/75‑1 (MQ=21)
                                                                                                                                     ccTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACa                                                                                        >  2:57654/1‑75 (MQ=21)
                                                                                                                                      cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAgcgc                                       >  8:133417/1‑123 (MQ=40)
                                                                                                                                      cTTATCCGGCCTACGGTTGGGCATCGCTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGc               >  7:414347/1‑147 (MQ=255)
                                                                                                                                       ttATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAg            >  4:229533/1‑149 (MQ=255)
                                                                                                                                       ttATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAg            >  4:398104/1‑149 (MQ=255)
                                                                                                                                            cGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGACAGCaaa        >  5:456094/1‑148 (MQ=255)
                                                                                                                                             ggCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCaaaaa      >  1:289895/1‑149 (MQ=255)
                                                                                                                                              gCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAAt     >  2:81935/1‑149 (MQ=255)
                                                                                                                                               ccTACGTTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCTCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATg    >  4:356549/1‑149 (MQ=255)
                                                                                                                                                 tACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGt  >  3:247848/1‑149 (MQ=255)
                                                                                                                                                   cGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGt  <  2:8856/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGT  >  NC_000913/4153422‑4153715

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: