Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,997,433 T→A 6.2% M84K (ATG→AAG)  yigB → 5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase; pyrimidine phosphatase; riboflavin synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,997,4330TA 6.2% 312.1 / 7.6 96M84K (ATG→AAG) yigB5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase; pyrimidine phosphatase; riboflavin synthesis
Reads supporting (aligned to +/- strand):  ref base T (40/50);  new base A (3/3);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACAGCTGGCGAAGAAATGGCCG  >  NC_000913/3997287‑3997572
                                                                                                                                                  |                                                                                                                                           
gcTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGc                                                                                                                                           >  2:370997/1‑149 (MQ=255)
     ccTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGAc                                                                                                                                      >  6:30755/1‑149 (MQ=255)
      cTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACg                                                                                                                                     >  3:95463/1‑149 (MQ=255)
            gCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGAAGCTCGACGCCGGGc                                                                                                                               >  7:305436/1‑149 (MQ=255)
            gCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCCGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGAAGCTCGACGCCGGGc                                                                                                                               >  1:292913/1‑149 (MQ=255)
             cAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCt                                                                                                                              >  1:108414/1‑149 (MQ=255)
              aaaaTTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCt                                                                                                                              >  3:9382/1‑148 (MQ=255)
              aaaaTTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCt                                                                                                                              >  4:285824/1‑148 (MQ=255)
              aaaaTTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCt                                                                                                                              >  4:296105/1‑148 (MQ=255)
                aaTTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAg                                                                                                                           <  8:227850/149‑1 (MQ=255)
                 aTTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGt                                                                                                                          <  7:212456/149‑1 (MQ=255)
                     tcatcCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGcc                                                                                                                       <  1:80020/148‑1 (MQ=255)
                     tcatcCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGAAGCTCGACGCCGGGCTGAGTGCCg                                                                                                                      <  8:305436/149‑1 (MQ=255)
                     tcatcCGGCGCTGCGCAGCTTCCAGAATGAAGATCCGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGAAGCTCGACGCCGGGCTGAGTGCCg                                                                                                                      <  2:292913/149‑1 (MQ=255)
                      catcCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACTCCGGGCTGAGTGCCg                                                                                                                      >  7:91303/1‑148 (MQ=255)
                       atcCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCgaa                                                                                                                    <  1:73948/149‑1 (MQ=255)
                       atcCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCgaa                                                                                                                    <  3:323979/149‑1 (MQ=255)
                          cggcggTGCGCAGCTTCCAGAATGAATATATGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCgaagaa                                                                                                                 <  4:9382/149‑1 (MQ=255)
                          cGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCgaagaa                                                                                                                 <  7:38861/149‑1 (MQ=255)
                              gcTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCa                                                                                                             <  1:65024/149‑1 (MQ=255)
                              gcTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCa                                                                                                             <  4:95463/149‑1 (MQ=255)
                              gcTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCa                                                                                                             <  5:30755/149‑1 (MQ=255)
                                ggcgcAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGt                                                                                                           <  5:44560/148‑1 (MQ=255)
                                  cgcAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGt                                                                                                           <  8:333595/147‑1 (MQ=255)
                                    cAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAg                                                                                                       <  5:148985/149‑1 (MQ=255)
                                    cAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAg                                                                                                       <  2:234311/149‑1 (MQ=255)
                                        ttCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGa                                                                                                                                       <  2:340931/113‑1 (MQ=255)
                                        ttCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGa                                                                                                                                       >  1:340931/1‑113 (MQ=255)
                                        ttCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGgcg                                                                                                    >  1:118443/1‑148 (MQ=255)
                                          ccAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGca                                                                                                  >  3:139929/1‑148 (MQ=255)
                                             gAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACg                                                                                              >  4:115425/1‑149 (MQ=255)
                                              aaTGAACATCTGCAACCCCGGGGCCTGGGGGTACGGGACGCGTAACCCGAGATTTTTCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACTCCGGGCTTAGTGCCGAAGAAGCCAGTGCAGGCGCACACgc                                                                                             <  7:264488/149‑1 (MQ=255)
                                                tGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACgcag                                                                                           >  7:350029/1‑149 (MQ=255)
                                                   aGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCaa                                                                                        <  3:147589/149‑1 (MQ=255)
                                                    gATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAAt                                                                                       <  2:118443/149‑1 (MQ=255)
                                                    gATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAAt                                                                                       >  4:247600/1‑149 (MQ=255)
                                                         gCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCgaagaa                                                                                                                 >  5:187697/1‑118 (MQ=255)
                                                         gCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCgaagaa                                                                                                                 <  6:187697/118‑1 (MQ=255)
                                                          cAACGCCTGCGCCAGGCGGTACTGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTAGACGCCGGGCTGAGTTCCTAAGAAGCCAGTGCAGGCTCACACGCAGCAATGATCaa                                                                                 >  8:196953/1‑149 (MQ=255)
                                                          cAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATg                                                                                      <  6:298424/144‑1 (MQ=255)
                                                          cAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGa                                                                                     <  7:258590/145‑1 (MQ=255)
                                                          cAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCaa                                                                                 >  4:340865/1‑149 (MQ=255)
                                                          cAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCGGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGa                                                                                     <  5:99183/145‑1 (MQ=255)
                                                           aaCGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCgaa                                                                                                                    >  7:287545/1‑113 (MQ=255)
                                                           aaCGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCgaa                                                                                                                    <  8:287545/113‑1 (MQ=255)
                                                            aCGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACt                                                                               >  8:67824/1‑149 (MQ=255)
                                                              gCCTGCGCCAGGCGGTACGGGATGAGGAACCAGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACTAGCGATGCTCGACTCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACttt                                                                             <  8:91303/149‑1 (MQ=255)
                                                              gCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACttt                                                                             <  1:204789/149‑1 (MQ=255)
                                                                cTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGc                                                                           <  3:285824/149‑1 (MQ=255)
                                                                     ccAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAAt                                                                      <  3:115425/149‑1 (MQ=255)
                                                                     ccAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAAt                                                                      >  5:301921/1‑149 (MQ=255)
                                                                      cAGGCGGTACGGGAAGGGGAAGCCGAGATTTATCACGACGTGACGCGGTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGTCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATg                                                                     <  4:215545/149‑1 (MQ=255)
                                                                      cAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATg                                                                     >  5:102644/1‑149 (MQ=255)
                                                                      cAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATg                                                                     >  3:231371/1‑149 (MQ=255)
                                                                      cAGGCGGTACGGGAAGCGGAACCCGAGATTTATAACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATg                                                                     >  5:328319/1‑149 (MQ=255)
                                                                       aGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCaa                                                                                        <  3:140048/129‑1 (MQ=255)
                                                                       aGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCaa                                                                                        >  4:140048/1‑129 (MQ=255)
                                                                        ggCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGgc                                                                   >  2:221709/1‑149 (MQ=255)
                                                                                  gAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATg                                                                     >  6:181072/1‑137 (MQ=255)
                                                                                  gAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATg                                                                     <  5:181072/137‑1 (MQ=255)
                                                                                  gAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCa                                                                >  8:228893/1‑142 (MQ=255)
                                                                                  gAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCa                                                                <  7:228893/142‑1 (MQ=255)
                                                                                  gAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAAt                                                         >  7:23142/1‑149 (MQ=255)
                                                                                      cGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGAc                                                     <  8:23142/149‑1 (MQ=255)
                                                                                      cGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGAc                                                     >  8:23166/1‑149 (MQ=255)
                                                                                      cGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGAc                                                     <  8:350029/149‑1 (MQ=255)
                                                                                      cGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGAc                                                     <  3:247600/149‑1 (MQ=255)
                                                                                           cccGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTccc                                                >  1:375510/1‑149 (MQ=255)
                                                                                            ccGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCg                                               >  8:91449/1‑149 (MQ=255)
                                                                                              gagaTTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCgca                                             <  4:139929/149‑1 (MQ=255)
                                                                                              gagaTTTATCACGACGTGACGCGCTGGCGGTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCg                                               >  3:229630/1‑147 (MQ=255)
                                                                                                     aTCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACt                                      >  1:74818/1‑149 (MQ=255)
                                                                                                        acgacgTGACGCGCTGGCGTTTTCGTTCGATTGAAGAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCAc                                   <  6:328319/149‑1 (MQ=255)
                                                                                                        acgacgTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGTTGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGgc                                                                   <  7:130924/117‑1 (MQ=255)
                                                                                                        acgacgTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGCTGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGgc                                                                   >  8:130924/1‑117 (MQ=255)
                                                                                                        acgacgTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCAc                                   <  3:340865/149‑1 (MQ=255)
                                                                                                          gacgTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGa                                 >  8:321595/1‑149 (MQ=255)
                                                                                                               gACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCt                            >  5:135726/1‑149 (MQ=255)
                                                                                                               gACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCt                            >  5:32013/1‑149 (MQ=255)
                                                                                                                  gcgcTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCAc                                   >  3:84450/1‑139 (MQ=255)
                                                                                                                  gcgcTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCAc                                   <  4:84450/139‑1 (MQ=255)
                                                                                                                    gcTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTaaaa                       <  2:108414/149‑1 (MQ=255)
                                                                                                                    gcTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTaaaa                       <  7:196953/149‑1 (MQ=255)
                                                                                                                        gCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACa                     <  6:102644/147‑1 (MQ=255)
                                                                                                                        gCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACa                     <  1:221709/147‑1 (MQ=255)
                                                                                                                        gCGTTTTCGTTCGATTGAACAAGCGATGCGCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACAGc                   <  8:339166/149‑1 (MQ=255)
                                                                                                                          gTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACAGCTg                 <  6:133472/149‑1 (MQ=255)
                                                                                                                                gttcgATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACAGCTGGCgaa            <  3:296105/148‑1 (MQ=255)
                                                                                                                                    gATTGAAGAAGCGATGCTGGACGCCGGTCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACAGCTGGCGAAGAAAt       <  2:74818/149‑1 (MQ=255)
                                                                                                                                    gATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACAGCTGGCGAAGAAAt       <  2:325437/149‑1 (MQ=255)
                                                                                                                                    gATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACAGCTGGCGAAGAAAt       >  7:295650/1‑149 (MQ=255)
                                                                                                                                    gATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACAGCTGGCGAAGAAAt       >  2:146317/1‑149 (MQ=255)
                                                                                                                                       tGAACAAGCGAAGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGc                                                              <  5:185068/91‑1 (MQ=255)
                                                                                                                                       tGAACAAGCGAAGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGc                                                              >  6:185068/1‑91 (MQ=255)
                                                                                                                                         aaCAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACAGCTGGCGAAGAAATGGCCg  <  6:135726/149‑1 (MQ=255)
                                                                                                                                         aaCAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACAGCTGGCGAAGAAATGGCCg  <  6:32013/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGCGGTACGGGAAGCGGAACCCGAGATTTATCACGACGTGACGCGCTGGCGTTTTCGTTCGATTGAACAAGCGATGCTCGACGCCGGGCTGAGTGCCGAAGAAGCCAGTGCAGGCGCACACGCAGCAATGATCAACTTTGCCAAATGGCGCAGCCGAATCGACGTCCCGCAGCAAACTCACGACACCTTAAAACAGCTGGCGAAGAAATGGCCG  >  NC_000913/3997287‑3997572

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