Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,366,015 Δ1 bp 5.4% coding (88/1140 nt) arnB → uridine 5'‑(beta‑1‑threo‑pentapyranosyl‑4‑ulose diphosphate) aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,366,0150A. 5.4% 258.9 / 9.4 92coding (88/1140 nt)arnBuridine 5'‑(beta‑1‑threo‑pentapyranosyl‑4‑ulose diphosphate) aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (54/30);  new base . (2/3);  total (57/35)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

CACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCT  >  NC_000913/2365869‑2366164
                                                                                                                                                   |                                                                                                                                                     
cacaACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCTCTGAAGTTAAAGAGGTTATCGAATCCGGTTGGATAACaa                                                                                                                                                    >  1:1043225/1‑150 (MQ=255)
 acaacaCTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAAc                                                                                                                                                   <  1:1700511/150‑1 (MQ=255)
 acaacaCTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAAc                                                                                                                                                   >  2:429199/1‑150 (MQ=255)
  caacaCTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAAcc                                                                                                                                                  <  1:1819583/150‑1 (MQ=255)
  caacaCTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGAGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAAcc                                                                                                                                                  >  1:948155/1‑150 (MQ=255)
  caacaCTCTATATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAAcc                                                                                                                                                  <  2:39744/150‑1 (MQ=255)
    aaacTCTTTATAGCATATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGGCTTTTTCGCGAGCAGCTATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCgg                                                                                                                                                <  2:1794014/148‑1 (MQ=255)
       ctctTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggatt                                                                                                                                             <  2:56445/150‑4 (MQ=255)
       ctctTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggatt                                                                                                                                             <  2:344603/150‑4 (MQ=255)
       ctctTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTcc                                                                                                                                             >  2:1577889/1‑150 (MQ=255)
        tctTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCg                                                                                                                                            <  2:139447/150‑1 (MQ=255)
           ttATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCgaag                                                                                                                                         >  2:1352368/1‑150 (MQ=255)
              taGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAAt                                                                                                                                      <  2:800903/150‑1 (MQ=255)
              taGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAAt                                                                                                                                      >  2:441162/1‑150 (MQ=255)
               aGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATc                                                                                                                                     <  2:943705/150‑1 (MQ=255)
                gCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCa                                                                                                                                    >  2:194681/1‑150 (MQ=255)
                       taAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCtt                                                                                                                             >  2:475012/1‑150 (MQ=255)
                        aaGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTg                                                                                                                            <  2:989278/150‑1 (MQ=255)
                        aaGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTGGGATCACAACCGGTACGAAGAATCAGGCGCTTg                                                                                                                            >  1:110359/1‑150 (MQ=255)
                         aGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTAGCTGCAGTTAACGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGa                                                                                                                           >  2:1118251/1‑150 (MQ=255)
                                 ggTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCtt                                                                                                                   >  2:1154606/1‑150 (MQ=255)
                                 ggTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCtt                                                                                                                   >  2:1055765/1‑150 (MQ=255)
                                  gTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCtgt                                                                                                                  >  2:1067204/1‑148 (MQ=255)
                                   tATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCtttt                                                                                                                 <  2:777066/150‑1 (MQ=255)
                                    aTTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCttttt                                                                                                                >  1:521663/1‑150 (MQ=255)
                                    aTTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCttttt                                                                                                                >  2:997408/1‑150 (MQ=255)
                                     ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                               <  2:1065120/150‑1 (MQ=255)
                                     ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                               <  2:1126818/150‑1 (MQ=255)
                                     ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                               <  1:947686/150‑1 (MQ=255)
                                     ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                               <  2:1190380/150‑1 (MQ=255)
                                       cAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                                                                             >  1:1392708/1‑150 (MQ=255)
                                       cAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                                                                             >  2:1715546/1‑150 (MQ=255)
                                            gCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTg                                                                                                        >  1:1691837/1‑150 (MQ=255)
                                                 aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                                                                   >  1:1480769/1‑150 (MQ=255)
                                                 aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                                                                   >  2:873725/1‑150 (MQ=255)
                                                 aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                                                                   >  1:1219773/1‑150 (MQ=255)
                                                 aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTACGAAGAATCAGGCGCTTGAGCTAGCTTTTTGCCAGTTGACggg                                                                                                   >  1:1131424/1‑150 (MQ=255)
                                                   gtaaaGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGGGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGGTACCCACCGGTCCGAAGTATGAGGCGGTTGAGCAAGCTTTTTGCCAGTTGACGGGaa                                                                                                 >  2:1447277/3‑150 (MQ=255)
                                                   gaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGTCGCTTGAGCAAGCTTTTTGCCAGTTGACGTGaa                                                                                                 >  2:1126425/1‑150 (MQ=255)
                                                   gaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGTGaa                                                                                                 >  2:49436/1‑150 (MQ=255)
                                                   gaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGaa                                                                                                 >  2:1132112/1‑150 (MQ=255)
                                                   gaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGATAAAGAGGTTCTCGAAACCGGTTGGATCACAACCGGTCCGATGAATCAGGCGCTTGAGCAAGATTTTTGCCAGTTGACGGGaa                                                                                                 >  2:1666793/1‑150 (MQ=255)
                                                     aaaGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAAt                                                                                               >  1:613088/1‑150 (MQ=255)
                                                      aaGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATc                                                                                              <  1:510908/150‑1 (MQ=255)
                                                        gCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAg                                                                                            >  2:517906/1‑150 (MQ=255)
                                                           aTGTGATAATTTTTGCCTTTTTCGCGACCATCAATGGGCGTGGAGGAACTCGCTGCATTTAAAGAGGTTCTCGAATCCGGTTTTATCTCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCAt                                                                                         <  1:1608865/150‑1 (MQ=255)
                                                             gTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGc                                                                                       >  2:1389806/1‑150 (MQ=255)
                                                             gTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGc                                                                                       >  1:421998/1‑150 (MQ=255)
                                                               cAGAATTTTTGCCTTTTTCGCGACCAGAAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGATATCAGCATGCCa                                                                                     >  1:142054/1‑150 (MQ=255)
                                                                aGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCAt                                                                                    >  2:1430099/1‑150 (MQ=255)
                                                                     ttttGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCgg                                                                               <  2:1047033/150‑1 (MQ=255)
                                                                        tGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCa                                                                            >  2:15197/1‑150 (MQ=255)
                                                                         gCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCATGAGCAAGCTTTTTGCCAGTTGACGGGAAATAAGCATGCCATCGCGGTCAg                                                                           >  2:1356321/1‑150 (MQ=255)
                                                                          ccTTTTTCGCGCCCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGt                                                                          >  1:1173889/1‑150 (MQ=255)
                                                                           cTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGtt                                                                         >  2:467643/1‑150 (MQ=255)
                                                                               ttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGc                                                                     <  1:1568138/150‑1 (MQ=255)
                                                                                      ccAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAccgcc                                                              >  2:345813/1‑150 (MQ=255)
                                                                                             tGGGCGTGGAGGAACTAGCTGCAGTTAAAGAGGTTCTCGAATCCGGTCGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTGTTGCCAGTTGACGGGAAATCAGCATGCCTTCGCGGTCAGTTCAGCCACCGCCGGAATGc                                                       <  1:302421/150‑1 (MQ=255)
                                                                                                gCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGCATGCata                                                    >  2:358506/1‑150 (MQ=255)
                                                                                                 cGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATTGCGGTCAGTTCAGCCACCGCCGGAATGCatat                                                   >  2:568651/1‑150 (MQ=255)
                                                                                                   tGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACAGCCGGAATGCATATCa                                                 <  1:1352737/150‑1 (MQ=255)
                                                                                                   tGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCACTTGACGCGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCa                                                 >  1:252472/1‑150 (MQ=255)
                                                                                                    ggaggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCAc                                                <  1:1036466/150‑1 (MQ=255)
                                                                                                    ggaggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCAc                                                >  1:1502977/1‑150 (MQ=255)
                                                                                                        gaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTa                                            <  2:450514/150‑1 (MQ=255)
                                                                                                             cGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGc                                       >  1:1713745/1‑150 (MQ=255)
                                                                                                                tGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCAAGTGGAGGGGAATTCCGCATGCCTGCGGGGTCTGTTCGGCCACCGCCGGATTCCATATCAGGCAAATGGCGtg                                    >  1:165882/1‑149 (MQ=255)
                                                                                                                     ttAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGaaaa                               >  2:1323204/1‑150 (MQ=255)
                                                                                                                      tAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAAt                              <  1:1792851/150‑1 (MQ=255)
                                                                                                                        aaGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTg                            <  1:873725/150‑1 (MQ=255)
                                                                                                                         agagGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTgg                           <  2:1034787/150‑1 (MQ=255)
                                                                                                                             gTTCTCGAATCCGGTTGGATCACAACCGGTCTGAACAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGATTGCATATCACGCTAATGGCGTTGAAAATTGGCAAg                       >  2:768649/1‑150 (MQ=255)
                                                                                                                             gTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAg                       >  1:834715/1‑150 (MQ=255)
                                                                                                                              ttCTCGAATCCGGTTGGATCATAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAgg                      >  1:860834/1‑150 (MQ=255)
                                                                                                                               tctcGAATCCGGTTGGATCACAACCGGGCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGCCGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAggg                     <  1:980336/150‑1 (MQ=255)
                                                                                                                                ctcGAATCCGGTTGGATCTCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGc                    <  2:964223/150‑1 (MQ=255)
                                                                                                                                ctcGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGc                    >  1:1622680/1‑150 (MQ=255)
                                                                                                                                   gAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGAt                 >  2:1471126/1‑150 (MQ=255)
                                                                                                                                   gAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGAt                 >  1:1202631/1‑150 (MQ=255)
                                                                                                                                   cggaccGGTTGTGATCCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTAGCAAGGGCGAt                 >  1:611437/5‑150 (MQ=255)
                                                                                                                                    ggaccGGTTGTGATCCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATg                >  2:1597651/4‑150 (MQ=255)
                                                                                                                                     aTCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGa               <  1:1568098/150‑1 (MQ=255)
                                                                                                                                       ccGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAg             >  2:470703/1‑150 (MQ=255)
                                                                                                                                        cGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGt            >  1:1252754/1‑150 (MQ=255)
                                                                                                                                         ggTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTg           <  1:629996/150‑1 (MQ=255)
                                                                                                                                         ggTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTg           <  2:1392708/150‑1 (MQ=255)
                                                                                                                                          gTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGATGGGAAATCAGCATGCCATCGCGGTCAGATCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGa          >  1:895069/1‑150 (MQ=255)
                                                                                                                                             ggATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTa       <  2:521663/150‑1 (MQ=255)
                                                                                                                                               gATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTAc      <  1:1252799/150‑1 (MQ=255)
                                                                                                                                               gATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTAc      <  1:917860/150‑1 (MQ=255)
                                                                                                                                               gATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTAc      <  1:1517141/150‑1 (MQ=255)
                                                                                                                                                 tCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGc    >  2:1061028/1‑150 (MQ=255)
                                                                                                                                                 ccacaACCGGGCCGACCAACCCGCCGCTTGACCATCCTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGcc   <  1:1252183/149‑1 (MQ=255)
                                                                                                                                                  cacaACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGcc   <  1:891537/150‑1 (MQ=255)
                                                                                                                                                   acaACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATAACACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCt  >  1:1130866/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCT  >  NC_000913/2365869‑2366164

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