Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,570 G→A 5.3% 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 5.3% 216.4 / 2.2 77intergenic (‑472/+126)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (40/32);  new base A (2/2);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

GGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAA  >  NC_000913/4153422‑4153718
                                                                                                                                                    |                                                                                                                                                    
ggcgCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACg                                                                                                                                                      <  2:233707/149‑1 (MQ=39)
 gcgcGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGGCGGATAAGATGCGCCGGTATCGCATCCGACATTATTGAGGGACTCGAGTTTGGCAACAGGGACGGATGCGGCGCGAGCGCCTTTTCAGGCGTACgg                                                                                                                                                     >  7:90119/1‑149 (MQ=11)
    cgAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTg                                                                                                                                                  >  3:79415/1‑149 (MQ=255)
     gAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                                                                                                                 <  4:79415/149‑1 (MQ=255)
          ccTTATCCTGACTACAGTTGGGCATGGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATc                                                                                                                                            <  4:103535/149‑1 (MQ=255)
            ttATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGt                                                                                                                                          >  8:60633/1‑149 (MQ=255)
              aTCCGGCCTACAGTTGTGCATCGTTTTAGTCACTGTAGTTCGTATAAGATGAGCAAGTATCTCATCAGAAATTATTGCGGAACTGGATTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTGGGATctggt                                                                                                                                        >  6:222257/1‑145 (MQ=11)
               tCCTGCATACAGTTGGGCATCGTTTGAGTCTCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                        <  7:104625/148‑1 (MQ=255)
                 ggtccTACAGTTGGGAAGCCTCTGACTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGa                                                                                                                                      <  6:261548/145‑1 (MQ=14)
                   gCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGAATCGGATGCGACAATATTGCGGCAGTGGTGTTTGCCAACAGTGCCGGATGCGGCGAGAGCGCCTTATCCGGACTACGGTTGGGCATCGTGTGAGTc                                                                                                                                   >  6:366294/1‑149 (MQ=17)
                                 aTCGTTTGAGTCCCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGTCGTCGTATGAGTCACTGTAGGTCGGa                                                                                                                      >  2:321155/1‑148 (MQ=255)
                                    gTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                        >  4:244369/1‑127 (MQ=255)
                                    gTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                        <  3:244369/127‑1 (MQ=255)
                                       tGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGAt                                                                                                                <  3:176597/148‑1 (MQ=255)
                                         aGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                        <  1:396422/122‑1 (MQ=40)
                                         aGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                                                        >  2:396422/1‑122 (MQ=40)
                                          gTCACTGTCGTTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCCCTGGAGTTGGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTCCGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgcgc                                                                                                            <  7:145988/149‑1 (MQ=255)
                                          gTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgcg                                                                                                             >  7:117572/1‑148 (MQ=40)
                                              cTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGGTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgc                                                                                                        <  5:366294/149‑1 (MQ=38)
                                                 tcggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                     <  4:185396/149‑1 (MQ=38)
                                                        gATAAGATGCTCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGAAACAGAGCAGGATGCGGAGCGAGCGCCATAACCGGACTACCGTTGCGCATCGTTTCAGTCACTGTAGGTAAGATAAGAAGAGCAAGAATCGCATccc                                                                                              >  6:196040/1‑148 (MQ=11)
                                                        gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCg                                                                                              >  7:253985/1‑149 (MQ=37)
                                                        gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGCCATCGTAACAGTCACTGAATGTCGGATAAGATCCGCAAGCATCGCATcc                                                                                               >  2:68630/1‑148 (MQ=11)
                                                        gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGAGCATCTTTTGAGTCACTGTAGTTCGGATAAGATGCGAAAGCATCTCATcc                                                                                               >  6:406263/1‑148 (MQ=11)
                                                             gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGt                                                                                                                              >  1:131284/1‑112 (MQ=32)
                                                             gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGt                                                                                                                              <  2:131284/112‑1 (MQ=32)
                                                             gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgc                                                                                                              >  6:287950/1‑128 (MQ=35)
                                                             gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgc                                                                                                              <  5:287950/128‑1 (MQ=35)
                                                             gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACat                                                                                          >  8:88167/1‑148 (MQ=37)
                                                              aTGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACatta                                                                                        <  7:88167/149‑1 (MQ=37)
                                                                gcgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGTGCAACGTTAGAGTAACTGTAGGACGGATAAGATGCGAAAGCATCGCATCCGACattaat                                                                                      >  2:268666/1‑147 (MQ=11)
                                                                 cgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTg                                                                                     <  8:253985/149‑1 (MQ=37)
                                                                  gcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                                                                                                                 <  2:213873/88‑1 (MQ=33)
                                                                  gcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                                                                                                                 >  7:340968/1‑88 (MQ=33)
                                                                  gcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                                                                                                                 <  8:340968/88‑1 (MQ=33)
                                                                  gcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                                                                                                                 >  1:213873/1‑88 (MQ=33)
                                                                    aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCTTCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCGGACATTATTGCgg                                                                                  >  8:365916/1‑149 (MQ=255)
                                                                     aGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGa                                                                                                                 >  6:351781/1‑117 (MQ=35)
                                                                     aGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGa                                                                                                                 <  5:351781/117‑1 (MQ=35)
                                                                       tATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCAc                                                                               >  7:382464/1‑149 (MQ=37)
                                                                       tATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCAc                                                                               >  3:195446/1‑149 (MQ=37)
                                                                         tCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                                                                                                                 <  8:90501/81‑1 (MQ=33)
                                                                         tCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                                                                                                                 >  7:90501/1‑81 (MQ=33)
                                                                         tCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTg                                                                             >  8:8283/1‑149 (MQ=37)
                                                                         tCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATACGGCCTACGTTTGTGCATCGTTTAAGTCACTGTATTTCGTATAAGATACGCAAGCATAGCATCCGACAATAATGCGCCACTg                                                                             >  8:311946/1‑149 (MQ=11)
                                                                          cGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                     >  3:330530/1‑124 (MQ=38)
                                                                          cGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc                                                                                                     <  4:330530/124‑1 (MQ=38)
                                                                             aTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgcg                                                                                                             <  1:394759/113‑1 (MQ=35)
                                                                             aTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgcg                                                                                                             >  2:394759/1‑113 (MQ=35)
                                                                              tCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTGGCATCGGTAGAGACACTGTAGGTCGGATAAGATGCGAAAGCATCGCCTCAGACATTATTGAGGCAcatgattt                                                                        >  1:422329/1‑142 (MQ=11)
                                                                                cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTg                                                                      <  7:60633/149‑1 (MQ=37)
                                                                                    attattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGTGAATCGTTTGAGTCACTGTAGGTCGGAAAAGAAGCGCAAGCATCGCATCCGAAATTATTGCGGCACGGGAGTTTGGaaa                                                                  >  2:152581/1‑149 (MQ=11)
                                                                                       attGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTTTCCGGACTACGGTTGGGGAGCGTTTGAGTCACTGTAGGTCGGATATGATGCGAAAGCATCGCATCCGGCATTATTGCGGCACTGGAGTTTGGGAACAg                                                               >  2:163783/1‑149 (MQ=11)
                                                                                           cGGCACTGCAGTTTGGCAACAGGGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCCCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGcc                                                           <  1:163783/149‑1 (MQ=255)
                                                                                               aCTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGAt                                                       >  3:215853/1‑149 (MQ=39)
                                                                                                 tGGAGTTTTGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTa                                                                                                                             <  7:62382/77‑1 (MQ=25)
                                                                                                   gAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCgg                                                                                  <  1:203936/118‑1 (MQ=17)
                                                                                                   gAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCgg                                                                                  >  2:203936/1‑118 (MQ=17)
                                                                                                      tttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGgcgc                                                >  1:124502/1‑149 (MQ=39)
                                                                                                            aaCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAgcgc                                          <  2:124502/149‑1 (MQ=39)
                                                                                                             aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGcc                                         >  7:9135/1‑149 (MQ=39)
                                                                                                                aggCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGTTGGGCATCGTGTGAGTCACTGTAGGTAGGATAAGATGCGCAAGCATCGGATGCGACATTATTGCGGCACTGGAGTTTTGGAACAGCGCCGGATGCGGCGCGAGCGccgtt                                      >  2:152883/3‑146 (MQ=14)
                                                                                                                 tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                     >  6:164866/1‑149 (MQ=39)
                                                                                                                 ggCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                     >  8:359136/2‑149 (MQ=14)
                                                                                                                  gCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATc                                    >  4:181849/1‑149 (MQ=39)
                                                                                                                  gCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                  >  4:137365/1‑119 (MQ=11)
                                                                                                                         gcggcgCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTATGAGACACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                             >  8:205995/1‑149 (MQ=255)
                                                                                                                              gcgAGCGCATTATGGTGTCTACGGTTGGGGTTCGGTTGTGGCACTTTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGttt                        <  5:196040/149‑1 (MQ=14)
                                                                                                                                 aGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTaa                     <  8:9135/149‑1 (MQ=255)
                                                                                                                                  gggcCTTATCCGGCCTACGGGGGGGGATCGTGTGTGTCACTGGCGGCCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTaa                     <  5:406263/146‑1 (MQ=14)
                                                                                                                                      cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCa                <  1:152581/149‑1 (MQ=255)
                                                                                                                                      cTTATCCGGCCTACGGTTGGGAATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCa                <  8:117572/149‑1 (MQ=255)
                                                                                                                                          tCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCa             >  4:62358/1‑148 (MQ=255)
                                                                                                                                             ggCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCaaaaa         <  8:382464/149‑1 (MQ=255)
                                                                                                                                             ggCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCaaaaa         <  7:8283/149‑1 (MQ=255)
                                                                                                                                              gCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATGCGGCCTACGTTTTAATGCCTGCAAatat        >  2:261425/1‑149 (MQ=255)
                                                                                                                                              gCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACattatt                                                                                      <  8:392017/71‑1 (MQ=255)
                                                                                                                                              gCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACattatt                                                                                      >  7:392017/1‑71 (MQ=21)
                                                                                                                                                   cGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGt     <  4:195446/147‑1 (MQ=255)
                                                                                                                                                    ggTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGaa  >  5:361218/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAA  >  NC_000913/4153422‑4153718

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: