Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,434 T→G 7.9% intergenic (‑336/+262) 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,4340TG 7.9% 208.2 / 8.8 76intergenic (‑336/+262)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base T (26/44);  new base G (3/3);  total (29/47)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

CTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGT  >  NC_000913/4153288‑4153582
                                                                                                                                                  |                                                                                                                                                    
cTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCCCATACGACAATATTGCGGCACACGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                                                    >  1:179597/1‑149 (MQ=25)
 tttAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                                                                   >  2:175261/1‑149 (MQ=255)
 tttAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                                                                   >  1:59038/1‑149 (MQ=255)
    aaaTAAGCGTAACTTATGGAAATTTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                                                 <  4:24251/148‑1 (MQ=255)
      aTAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                                              >  7:66584/1‑149 (MQ=255)
         aGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                                                                                                                                           <  3:126750/149‑1 (MQ=255)
         aGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                                                                                                                                           <  7:247875/149‑1 (MQ=255)
            gTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAg                                                                                                                                        <  2:76957/149‑1 (MQ=255)
             tAACTTATGGAAATGTTAAAAAATCTCCCCAAGTAACACCACAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGt                                                                                                                                       <  6:346694/149‑1 (MQ=255)
             tAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGt                                                                                                                                       <  3:104769/149‑1 (MQ=255)
                cTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                                                                    <  4:332041/149‑1 (MQ=255)
                cTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                                                                    <  3:403/149‑1 (MQ=255)
                cTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                                                                    <  1:19029/149‑1 (MQ=255)
                   aTGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGc                                                                                                                                  <  1:175261/148‑1 (MQ=255)
                    tGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCa                                                                                                                                 <  8:44451/148‑1 (MQ=255)
                      gAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATc                                                                                                                               <  3:344640/148‑1 (MQ=255)
                         aTGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGttt                                                                                                                           <  4:256856/149‑1 (MQ=255)
                         aTGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGttt                                                                                                                           <  5:99886/149‑1 (MQ=255)
                            ttAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGGTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGa                                                                                                                         <  3:264107/148‑1 (MQ=255)
                                aaaaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCAc                                                                                                                    <  1:279142/149‑1 (MQ=255)
                                 aaaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACt                                                                                                                   <  2:186483/149‑1 (MQ=255)
                                   aTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgt                                                                                                                 >  6:351674/1‑149 (MQ=255)
                                    tCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtc                                                                                                                >  1:270745/1‑149 (MQ=255)
                                       ccccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtaggt                                                                                                             >  1:273486/1‑149 (MQ=255)
                                        cccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggtc                                                                                                            <  7:28712/149‑1 (MQ=255)
                                        cccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTAGGCAACAGTGACGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggtc                                                                                                            >  1:275925/1‑149 (MQ=255)
                                           aaGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGa                                                                                                         <  5:351674/149‑1 (MQ=255)
                                            aGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGAt                                                                                                        <  3:214478/149‑1 (MQ=255)
                                             gTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATa                                                                                                       <  8:119660/149‑1 (MQ=255)
                                              tAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGAt                                                                                                        >  3:180795/1‑147 (MQ=255)
                                                acacCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCAGCGTTTGAGTCACTGTCGGTCGGATAAg                                                                                                     <  1:165506/148‑1 (MQ=255)
                                                  acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTAAAGTTGGGCATCCTTTGCGTCACTCTCCGTCGGATAAGaag                                                                                                  >  7:337150/1‑147 (MQ=14)
                                                  acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTAGTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATg                                                                                                  >  7:350206/1‑149 (MQ=255)
                                                    cAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCGTGGTTTGAGTCACTGTCGGTCGGATAAGATgcg                                                                                                >  7:177902/1‑149 (MQ=255)
                                                       aGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCaa                                                                                             <  3:312738/149‑1 (MQ=39)
                                                          tgtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGt                                                                                           >  3:134124/1‑148 (MQ=38)
                                                            tAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCAc                                                                                                                    <  3:310353/121‑1 (MQ=34)
                                                            tAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCAc                                                                                                                    >  4:310353/1‑121 (MQ=34)
                                                             aGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATc                                                                                        <  7:321805/148‑1 (MQ=38)
                                                                 cGGATAAGATGCTCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATc                                                                                   <  1:102771/149‑1 (MQ=37)
                                                                     tAAGATGCGCAAGTATCGAATCCGACTTTATTGCGGCACTGGAGTTGTCCAACATTGCCGGATGCTTCGCGAGCGCCTTATCCGGACTACAGTTGGGCATCTTTTGAGTCTCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGAc                                                                               <  2:275925/149‑1 (MQ=255)
                                                                      aaGCTGCGCAAGTCTCGCATCCGACATTATTGCGGCACGGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTGTGAGTCACTGTCGGGCGGATAAGATGCGCAAGTATCGCAt                                                                                    <  7:366533/143‑1 (MQ=11)
                                                                      aaGATGCGCAAGTATCGCATCCGACATTAATGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTAGCCGGCATACAGTTGGGCATCGTGGGAGTCGCTGTCGGTCGGATAAGATGCGCAAGTATCGCag                                                                                    >  8:366533/1‑142 (MQ=11)
                                                                        gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACa                                                                              <  8:33425/147‑1 (MQ=37)
                                                                        gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACa                                                                              >  7:33425/1‑147 (MQ=37)
                                                                            cgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACattatt                                                                         <  6:153593/148‑1 (MQ=37)
                                                                               aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAg                                                                                                     >  8:323630/1‑117 (MQ=35)
                                                                               aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAg                                                                                                     <  7:323630/117‑1 (MQ=35)
                                                                                    tagaaTCCGACATTATTGCGGCACTGGAGTTTGTCAACAGTGCCGTATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTg                                                                <  1:235772/145‑1 (MQ=255)
                                                                                    tCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                                                                 >  7:335274/1‑116 (MQ=35)
                                                                                    ccgatctcGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTCGCCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTg                                                                <  5:288876/142‑1 (MQ=255)
                                                                                      gCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGGTGCGGCGCGAGCGGCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                                              <  7:14754/149‑1 (MQ=255)
                                                                                            gACATTATTGCGGCACTGGAGATTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATAAGGCGTACAGTTGGGCATCGGATGAGAAAATGTCGTTAGGATAGGATGCGCGGGGGTCGCCTGCGACAGTGTTGGGGCGCTGGAgatcgg                                                        >  6:288876/1‑144 (MQ=2)
                                                                                                  attGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGCTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACa                                                   <  4:180795/148‑1 (MQ=37)
                                                                                                      cGGCACTGGAGTTTGGCAACTGTGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGGCGGATAAGATGCGCAAGTATCGCATCCGACTTTATTGCGGCACTGGAGTTTGGCAACAGTGcc                                              <  2:385157/149‑1 (MQ=255)
                                                                                                            tGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCAt                                                                                    >  7:241528/1‑105 (MQ=17)
                                                                                                              gAGTTTGGCAACAGTTCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCTCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcg                                       <  8:157176/148‑1 (MQ=255)
                                                                                                               aGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcggc                                     <  4:134124/149‑1 (MQ=39)
                                                                                                                    ggCTACAGTGCAGGATGCGGCTCGAGCTCCTTATCATTCATACATTTGGGAATCTTTTGATTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGa                                 <  4:109672/148‑1 (MQ=14)
                                                                                                                     gCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgc                               <  8:110623/149‑1 (MQ=39)
                                                                                                                       ctgtaggCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcgc                             <  4:327867/143‑1 (MQ=14)
                                                                                                                       aaCAGTGGAGTATGCGGCTCGATCGCCTTATCAGGCCTACAGTTGGGCATAGTTTGAGTAACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcgc                             <  8:238498/149‑1 (MQ=14)
                                                                                                                       aaCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGt                                                            >  1:351282/1‑118 (MQ=39)
                                                                                                                        aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                                              >  1:326107/1‑115 (MQ=39)
                                                                                                                        aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc                            >  7:301419/1‑149 (MQ=39)
                                                                                                                        aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc                            <  8:66584/149‑1 (MQ=39)
                                                                                                                        aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGAACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGTGcc                            >  8:83259/1‑149 (MQ=255)
                                                                                                                             gCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                             >  2:372727/1‑127 (MQ=14)
                                                                                                                             gCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                             <  1:372727/127‑1 (MQ=14)
                                                                                                                             gCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACAGTTGGGCATCGTGTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATc                       >  3:262465/1‑149 (MQ=14)
                                                                                                                             gCCGGATAAGGCGCTCGCGCCGCATCCGGCCGACAGTTGAGCATCGTGTGAGTCCCTGTCGGGCGGATAAGATGCGCAAGTATCGCATCCGACCTTATTGCGGCCCTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATc                       >  5:165547/1‑149 (MQ=14)
                                                                                                                               atgaagcggccCGAGCTCCTTATCCGGCCTACAGTTGCGCATCCACTGATTCACTGTAGGTCGGATAAGATTCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATTCGGCGCGAGCGCCTTATcc                      <  7:182885/143‑1 (MQ=11)
                                                                                                                                          ctcgcgcCGCATCCGGCCTACCGTTGGGCCTCGTTTGAGTCACTGTCGGTCGTATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGCCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTg          <  4:262465/146‑1 (MQ=14)
                                                                                                                                           gAGCGCCTTATCCGGCCTACAGTTGCGCAGCGTGTGAGTCCCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg         <  8:177902/149‑1 (MQ=255)
                                                                                                                                            aGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTggg        >  1:10416/1‑149 (MQ=255)
                                                                                                                                                 cTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTgg                                                        >  6:272193/1‑96 (MQ=18)
                                                                                                                                                  ttATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGt  >  7:62505/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGT  >  NC_000913/4153288‑4153582

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: