Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,570 G→A 7.6% 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 7.6% 265.8 / 9.7 94intergenic (‑472/+126)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (47/37);  new base A (3/4);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

GGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGA  >  NC_000913/4153422‑4153717
                                                                                                                                                    |                                                                                                                                                   
ggcgCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACg                                                                                                                                                     >  4:294837/1‑149 (MQ=39)
 gcgcGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGCTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACgg                                                                                                                                                    >  5:338220/1‑149 (MQ=255)
     gAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                                                                                                                <  3:294837/149‑1 (MQ=255)
            ttATCCGGCCTACAGTTCGGCATCGGTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGGATCGCATGCGACATTATTGCGTCACTGGAGTTTGGCAACAGGGCCGGCTGCGGCGCGAGCGCCGTATCCGGCCGACGGGTGGGCAGCgg                                                                                                                                         >  2:136370/1‑148 (MQ=11)
               tCCGTCCTACAGTAGGGCAGCGTATGATTCACTGTAGGTCGGATAAGTTGCTCAAGTATCGCTTCCGACATTATTGCGGCACTGGAGTTTGGCTACAGTGCCGGATGCGGCGCGAGAGCCTTATCCGGCCTTCTGTTGGGCATCGttt                                                                                                                                       <  6:143110/148‑1 (MQ=14)
                         aGTTTGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCTTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGAGAGAGGGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGt                                                                                                                             <  1:92464/148‑1 (MQ=17)
                         aGTTGGGCATCGTTTGAGTCACTGTAGGTACGAAAAGATGCGCAAGTATCGAATCCTAAATTATTGCGGCACTGGAGTTTGGCAACAGTGACGGATTCGAAGCGAGCGCCTTATCCGGCCTAGGGGTGGGCATCGTTTGAGTGACTGt                                                                                                                             >  2:92464/1‑148 (MQ=14)
                                cATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGCGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGa                                                                                                                     >  6:134257/1‑149 (MQ=255)
                                aatagTTTGAGTCTCTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGa                                                                                                                     <  5:345186/145‑1 (MQ=255)
                                 aTCGTTTGCTTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGAt                                                                                                                    <  6:130176/149‑1 (MQ=255)
                                        gAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGt                                                                                                                                                   <  6:303222/111‑1 (MQ=40)
                                        gAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGt                                                                                                                                                   >  5:303222/1‑111 (MQ=40)
                                        gAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgc                                                                                                             >  5:308657/1‑149 (MQ=255)
                                          gTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCGGCGTGTGAGTGGCTGTAGGTCGGAGAAGATgcgc                                                                                                           >  4:72844/1‑149 (MQ=14)
                                            cccTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCaa                                                                                                         <  3:111302/147‑1 (MQ=38)
                                              cTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAg                                                                                                                           >  6:33903/1‑129 (MQ=38)
                                              cTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATACGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAg                                                                                                        >  7:115058/1‑148 (MQ=38)
                                               tGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgca                                                                                                      >  7:268247/1‑149 (MQ=38)
                                                   ggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAg                                                                                                                           <  1:326768/124‑1 (MQ=34)
                                                   ggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                    <  6:168923/147‑1 (MQ=37)
                                                   ggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                    >  5:168923/1‑147 (MQ=37)
                                                    gtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTGTGAGTCGCTGTAGGTCGGATAAGATGCGCAAgcatcgca                                                                                                 >  1:36604/1‑149 (MQ=255)
                                                     tcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCaa                                                                                                         <  5:294345/140‑1 (MQ=34)
                                                     tcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGCCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCaa                                                                                                         >  6:294345/1‑140 (MQ=21)
                                                        gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                                                                                                                >  2:251726/1‑98 (MQ=33)
                                                        gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                                                                                                                <  1:251726/98‑1 (MQ=33)
                                                        gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCCCTGTAGGTCGGAGAAGATGCGCACGCCTCGCCTcc                                                                                              >  6:264255/1‑148 (MQ=11)
                                                          tAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGAc                                                                                           <  6:308657/149‑1 (MQ=37)
                                                             gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTGTCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACatt                                                                                        >  6:156517/1‑149 (MQ=37)
                                                                    aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                       <  7:391501/95‑1 (MQ=33)
                                                                    aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                       >  8:391501/1‑95 (MQ=33)
                                                                    aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGa                                                                                                                                     <  1:199201/97‑1 (MQ=33)
                                                                    aaGTATCGCATCCAACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGa                                                                                                                                     >  2:199201/1‑97 (MQ=18)
                                                                            cATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTAACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAg                                                                         >  2:138134/1‑149 (MQ=37)
                                                                             aTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGTATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTgg                                                                           <  4:343063/146‑1 (MQ=37)
                                                                             aTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGTATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTgg                                                                           >  3:343063/1‑146 (MQ=37)
                                                                             aTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGt                                                                        >  8:264037/1‑149 (MQ=37)
                                                                             aTCCGACATTATTGCGGCACTGGAGTTTGACAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCa                                                                               <  2:35514/142‑1 (MQ=37)
                                                                             aTCCGACATTATTGCGGCACTGGAGTTTGACAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCa                                                                               >  1:35514/1‑142 (MQ=37)
                                                                              tCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTGTCATCGTTTGAGTAACTGTAGGTCGGATAAGATGCGCAAGCCTCGCAGCAGACATTATTGCGGGACTGGAGtt                                                                       >  6:21321/1‑149 (MQ=11)
                                                                                cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGttt                                                                      >  4:164004/1‑148 (MQ=37)
                                                                                cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGttt                                                                      <  3:164004/148‑1 (MQ=37)
                                                                                cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGAGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTg                                                                     >  2:3837/1‑149 (MQ=37)
                                                                                 gACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTgg                                                                    >  8:218371/1‑149 (MQ=37)
                                                                                 gACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCATAAAAATCACTGAAGGACGGATAAGAAGCGCAAGCATCGCATTCCACATAATTGAGGCACTAGAGATTg                                                                     >  3:18799/1‑148 (MQ=11)
                                                                                   cATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGa                                                                                            <  7:168287/123‑1 (MQ=38)
                                                                                   cATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGa                                                                                            >  8:168287/1‑123 (MQ=38)
                                                                                    attattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                 <  1:3837/149‑1 (MQ=37)
                                                                                    attattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                 >  4:56416/1‑149 (MQ=37)
                                                                                    attattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTAATGCGGCACTTGAGTTTGGCaa                                                                 >  3:304365/1‑149 (MQ=255)
                                                                                     ttattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTAACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                 >  8:81258/1‑148 (MQ=37)
                                                                                                cTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATg                                                     >  6:223628/1‑149 (MQ=39)
                                                                                                 tGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATgc                                                    >  6:243838/1‑149 (MQ=39)
                                                                                                 tGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATcc                                                                                              >  3:19289/1‑107 (MQ=17)
                                                                                                 tGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATcc                                                                                              <  4:19289/107‑1 (MQ=17)
                                                                                                   gAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGa                                                                                                                                     <  3:22029/66‑1 (MQ=35)
                                                                                                      tttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCg                                                                                                                                          >  4:79608/1‑58 (MQ=35)
                                                                                                      tttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCg                                                                                                                                          <  3:79608/58‑1 (MQ=35)
                                                                                                      tttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGgcgc                                               >  7:322555/1‑149 (MQ=39)
                                                                                                       ttGGCAACAGTTCCGGTTGCGGCGCTAGCGCCTTATCCGGGATACTGTTGGGCATTTTTTGATTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGgcgc                                               <  2:198301/148‑1 (MQ=11)
                                                                                                         ggCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGa                                             <  8:115058/148‑1 (MQ=39)
                                                                                                         ggCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGa                                             <  5:156517/148‑1 (MQ=39)
                                                                                                         ggCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAg                                            >  8:220411/1‑149 (MQ=39)
                                                                                                               aGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCCAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCt                                       >  1:207394/1‑148 (MQ=255)
                                                                                                                 tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATTAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                    >  3:253800/1‑149 (MQ=255)
                                                                                                                 tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                    <  8:268247/149‑1 (MQ=39)
                                                                                                                    cGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCg                                 <  5:134257/149‑1 (MQ=39)
                                                                                                                    cGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCg                                 <  7:218371/149‑1 (MQ=39)
                                                                                                                      gATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGc                               <  8:322555/149‑1 (MQ=39)
                                                                                                                       aTGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcgcatc                                                                                               <  5:55186/84‑1 (MQ=18)
                                                                                                                       aTGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcgcatc                                                                                               >  6:55186/1‑84 (MQ=18)
                                                                                                                         gcggcgCGCCCGCCTTATCCGGCCTACGGTGGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGTGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                            <  1:136370/149‑1 (MQ=255)
                                                                                                                         gcggcgCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                    <  8:23644/141‑1 (MQ=39)
                                                                                                                         gcggcgCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                    >  7:23644/1‑141 (MQ=39)
                                                                                                                           ggcgCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACat                                                                                         <  1:284682/86‑1 (MQ=255)
                                                                                                                           ggcgCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACat                                                                                         >  2:284682/1‑86 (MQ=18)
                                                                                                                              gcgAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATTCTCAAGAATCGCATACGCCATTATTTCTGCACTGGAGTTTGGAAACAGCGCCGGATGAGGCGAGAGCGCCTTATCCGGCCTACGttt                       >  7:271189/1‑149 (MQ=14)
                                                                                                                               agAGCGCCTTATCCGGCCTACGGTTGGCAATCGTTTAAGTCACTGTAAGTCGGATAAGAAGCCCAAGCATCGCATCCGACAAAATTTCTGAACTGGATTTTTGCAACACCGACAGATGAGTCACTACCGCCTTATCCGGCCAACCttt                       >  8:142325/2‑148 (MQ=255)
                                                                                                                                  gcgcCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAAt                   >  2:153866/1‑149 (MQ=255)
                                                                                                                                  gcgcCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcg                                                                                                   <  2:292512/69‑1 (MQ=25)
                                                                                                                                    gcCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGc                 <  7:81258/149‑1 (MQ=255)
                                                                                                                                     ccTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGc                 >  5:87136/1‑148 (MQ=255)
                                                                                                                                      cTTATCCGGCCTTCGGTTGGGCATCGTTTGCCGCACTGTAGGTCGGATAATATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCt                             <  4:32199/135‑1 (MQ=14)
                                                                                                                                      cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                    >  3:189711/1‑64 (MQ=255)
                                                                                                                                      cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                    <  4:189711/64‑1 (MQ=255)
                                                                                                                                      cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGc                                                                   >  5:37562/1‑97 (MQ=40)
                                                                                                                                      cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGc                                                                   <  6:37562/97‑1 (MQ=40)
                                                                                                                                      cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCt                             >  3:32199/1‑135 (MQ=40)
                                                                                                                                      cTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                 <  8:394725/99‑1 (MQ=18)
                                                                                                                                         aTCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGc             >  1:233768/1‑148 (MQ=255)
                                                                                                                                               ccTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCCTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATg      >  7:134510/1‑149 (MQ=255)
                                                                                                                                               ccTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATg      <  3:72844/149‑1 (MQ=255)
                                                                                                                                               acTACTGTTGGGCCTGGTTTGCGTCACTGTCGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATg      <  5:307004/148‑1 (MQ=17)
                                                                                                                                                cTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATgg     <  7:264037/149‑1 (MQ=255)
                                                                                                                                                   cGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGt    <  5:223628/147‑1 (MQ=255)
                                                                                                                                                   cGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGa  <  7:142325/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGA  >  NC_000913/4153422‑4153717

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: