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

CACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCT  >  NC_000913/2365869‑2366164
                                                                                                                                                   |                                                                                                                                                     
cacaACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCTCTGAAGTTAAAGAGGTTATCGAATCCGGTTG‑GATAACaa                                                                                                                                                    >  1:1043225/1‑150 (MQ=255)
 acaacaCTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAAc                                                                                                                                                   <  1:1700511/150‑1 (MQ=255)
 acaacaCTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAAc                                                                                                                                                   >  2:429199/1‑150 (MQ=255)
  caacaCTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAAcc                                                                                                                                                  <  1:1819583/150‑1 (MQ=255)
  caacaCTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGAGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAAcc                                                                                                                                                  >  1:948155/1‑150 (MQ=255)
  caacaCTCTATATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAAcc                                                                                                                                                  <  2:39744/150‑1 (MQ=255)
    aaacTCTTTATAGCATATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGGCTTTTTCGCGAGCAGCTATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCgg                                                                                                                                                <  2:1794014/148‑1 (MQ=255)
       ctctTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggatt                                                                                                                                             <  2:56445/150‑4 (MQ=255)
       ctctTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggatt                                                                                                                                             <  2:344603/150‑4 (MQ=255)
       ctctTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTcc                                                                                                                                             >  2:1577889/1‑150 (MQ=255)
        tctTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCg                                                                                                                                            <  2:139447/150‑1 (MQ=255)
           ttATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCgaag                                                                                                                                         >  2:1352368/1‑150 (MQ=255)
              taGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAAt                                                                                                                                      <  2:800903/150‑1 (MQ=255)
              taGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAAt                                                                                                                                      >  2:441162/1‑150 (MQ=255)
               aGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATc                                                                                                                                     <  2:943705/150‑1 (MQ=255)
                gCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCa                                                                                                                                    >  2:194681/1‑150 (MQ=255)
                       taAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCtt                                                                                                                             >  2:475012/1‑150 (MQ=255)
                        aaGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTg                                                                                                                            <  2:989278/150‑1 (MQ=255)
                        aaGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTGG‑GATCACAACCGGTACGAAGAATCAGGCGCTTg                                                                                                                            >  1:110359/1‑150 (MQ=255)
                         aGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTAGCTGCAGTTAACGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGa                                                                                                                           >  2:1118251/1‑150 (MQ=255)
                                 ggTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCtt                                                                                                                   >  2:1154606/1‑150 (MQ=255)
                                 ggTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCtt                                                                                                                   >  2:1055765/1‑150 (MQ=255)
                                  gTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCtgt                                                                                                                  >  2:1067204/1‑148 (MQ=255)
                                   tATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCtttt                                                                                                                 <  2:777066/150‑1 (MQ=255)
                                    aTTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCttttt                                                                                                                >  1:521663/1‑150 (MQ=255)
                                    aTTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCttttt                                                                                                                >  2:997408/1‑150 (MQ=255)
                                     ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                               <  2:1065120/150‑1 (MQ=255)
                                     ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                               <  2:1126818/150‑1 (MQ=255)
                                     ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                               <  1:947686/150‑1 (MQ=255)
                                     ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                               <  2:1190380/150‑1 (MQ=255)
                                       cAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                                                                             >  1:1392708/1‑150 (MQ=255)
                                       cAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                                                                             >  2:1715546/1‑150 (MQ=255)
                                            gCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTg                                                                                                        >  1:1691837/1‑150 (MQ=255)
                                                 aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                                                                   >  1:1480769/1‑150 (MQ=255)
                                                 aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                                                                   >  2:873725/1‑150 (MQ=255)
                                                 aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                                                                   >  1:1219773/1‑150 (MQ=255)
                                                 aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTACGAAGAATCAGGCGCTTGAGCTAGCTTTTTGCCAGTTGACggg                                                                                                   >  1:1131424/1‑150 (MQ=255)
                                                   gtaaaGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGGGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GGTACCCACCGGTCCGAAGTATGAGGCGGTTGAGCAAGCTTTTTGCCAGTTGACGGGaa                                                                                                 >  2:1447277/3‑150 (MQ=255)
                                                   gaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGTCGCTTGAGCAAGCTTTTTGCCAGTTGACGTGaa                                                                                                 >  2:1126425/1‑150 (MQ=255)
                                                   gaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGTGaa                                                                                                 >  2:49436/1‑150 (MQ=255)
                                                   gaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGaa                                                                                                 >  2:1132112/1‑150 (MQ=255)
                                                   gaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGATAAAGAGGTTCTCGAAACCGGTTG‑GATCACAACCGGTCCGATGAATCAGGCGCTTGAGCAAGATTTTTGCCAGTTGACGGGaa                                                                                                 >  2:1666793/1‑150 (MQ=255)
                                                     aaaGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAAt                                                                                               >  1:613088/1‑150 (MQ=255)
                                                      aaGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATc                                                                                              <  1:510908/150‑1 (MQ=255)
                                                        gCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAg                                                                                            >  2:517906/1‑150 (MQ=255)
                                                           aTGTGATAATTTTTGCCTTTTTCGCGACCATCAATGGGCGTGGAGGAACTCGCTGCATTTAAAGAGGTTCTCGAATCCGGTTT‑TATCTCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCAt                                                                                         <  1:1608865/150‑1 (MQ=255)
                                                             gTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGc                                                                                       >  2:1389806/1‑150 (MQ=255)
                                                             gTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGc                                                                                       >  1:421998/1‑150 (MQ=255)
                                                               cAGAATTTTTGCCTTTTTCGCGACCAGAAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGATATCAGCATGCCa                                                                                     >  1:142054/1‑150 (MQ=255)
                                                                aGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCAt                                                                                    >  2:1430099/1‑150 (MQ=255)
                                                                     ttttGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCgg                                                                               <  2:1047033/150‑1 (MQ=255)
                                                                        tGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCa                                                                            >  2:15197/1‑150 (MQ=255)
                                                                         gCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCATGAGCAAGCTTTTTGCCAGTTGACGGGAAATAAGCATGCCATCGCGGTCAg                                                                           >  2:1356321/1‑150 (MQ=255)
                                                                          ccTTTTTCGCGCCCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGt                                                                          >  1:1173889/1‑150 (MQ=255)
                                                                           cTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGtt                                                                         >  2:467643/1‑150 (MQ=255)
                                                                               ttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGc                                                                     <  1:1568138/150‑1 (MQ=255)
                                                                                      ccAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAccgcc                                                              >  2:345813/1‑150 (MQ=255)
                                                                                             tGGGCGTGGAGGAACTAGCTGCAGTTAAAGAGGTTCTCGAATCCGGTCG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTGTTGCCAGTTGACGGGAAATCAGCATGCCTTCGCGGTCAGTTCAGCCACCGCCGGAATGc                                                       <  1:302421/150‑1 (MQ=255)
                                                                                                gCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGCATGCata                                                    >  2:358506/1‑150 (MQ=255)
                                                                                                 cGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATTGCGGTCAGTTCAGCCACCGCCGGAATGCatat                                                   >  2:568651/1‑150 (MQ=255)
                                                                                                   tGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACAGCCGGAATGCATATCa                                                 <  1:1352737/150‑1 (MQ=255)
                                                                                                   tGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCACTTGACGCGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCa                                                 >  1:252472/1‑150 (MQ=255)
                                                                                                    ggaggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCAc                                                <  1:1036466/150‑1 (MQ=255)
                                                                                                    ggaggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCAc                                                >  1:1502977/1‑150 (MQ=255)
                                                                                                        gaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTa                                            <  2:450514/150‑1 (MQ=255)
                                                                                                             cGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGc                                       >  1:1713745/1‑150 (MQ=255)
                                                                                                                tGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCAAGTGGAGGGGAATTCCGCATGCCTGCGGGGTCTGTTCGGCCACCGCCGGATTCCATATCAGGCAAATGGCGtg                                    >  1:165882/1‑149 (MQ=255)
                                                                                                                     ttAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGaaaa                               >  2:1323204/1‑150 (MQ=255)
                                                                                                                      tAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAAt                              <  1:1792851/150‑1 (MQ=255)
                                                                                                                        aaGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTg                            <  1:873725/150‑1 (MQ=255)
                                                                                                                         agagGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTgg                           <  2:1034787/150‑1 (MQ=255)
                                                                                                                             gTTCTCGAATCCGGTTG‑GATCACAACCGGTCTGAACAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGATTGCATATCACGCTAATGGCGTTGAAAATTGGCAAg                       >  2:768649/1‑150 (MQ=255)
                                                                                                                             gTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAg                       >  1:834715/1‑150 (MQ=255)
                                                                                                                              ttCTCGAATCCGGTTG‑GATCATAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAgg                      >  1:860834/1‑150 (MQ=255)
                                                                                                                               tctcGAATCCGGTTG‑GATCACAACCGGGCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGCCGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAggg                     <  1:980336/150‑1 (MQ=255)
                                                                                                                                ctcGAATCCGGTTG‑GATCTCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGc                    <  2:964223/150‑1 (MQ=255)
                                                                                                                                ctcGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGc                    >  1:1622680/1‑150 (MQ=255)
                                                                                                                                   gAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGAt                 >  2:1471126/1‑150 (MQ=255)
                                                                                                                                   gAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGAt                 >  1:1202631/1‑150 (MQ=255)
                                                                                                                                   cggaccGGTTGTGATC‑CAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTAGCAAGGGCGAt                 >  1:611437/5‑150 (MQ=255)
                                                                                                                                    ggaccGGTTGTGATC‑CAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATg                >  2:1597651/4‑150 (MQ=255)
                                                                                                                                     aTCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGa               <  1:1568098/150‑1 (MQ=255)
                                                                                                                                       ccGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAg             >  2:470703/1‑150 (MQ=255)
                                                                                                                                        cGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGt            >  1:1252754/1‑150 (MQ=255)
                                                                                                                                         ggTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTg           <  1:629996/150‑1 (MQ=255)
                                                                                                                                         ggTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTg           <  2:1392708/150‑1 (MQ=255)
                                                                                                                                          gTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGATGGGAAATCAGCATGCCATCGCGGTCAGATCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGa          >  1:895069/1‑150 (MQ=255)
                                                                                                                                             g‑gATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTa       <  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)
                                                                                                                                                 ccacaACCGGGCCGACCAACCCGCCGCTTGA‑CCATCCTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGcc   <  1:1252183/149‑1 (MQ=255)
                                                                                                                                                  cacaACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGcc   <  1:891537/150‑1 (MQ=255)
                                                                                                                                                   acaACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATAACACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCt  >  1:1130866/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCT  >  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: ‑