Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,366,010:1 +T 5.5% coding (83/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,0101.T 5.5% 362.3 / 15.7 91coding (83/1140 nt)arnBuridine 5'‑(beta‑1‑threo‑pentapyranosyl‑4‑ulose diphosphate) aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base . (55/31);  new base T (2/3);  total (57/34)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

AAAATATACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTA  >  NC_000913/2365861‑2366159
                                                                                                                                                      |                                                                                                                                                     
aaaaTATACACAACACTCTTTATAGCAAATATTACTGGACAGGTATTCAATGTCGGAAGGAAAAGCTATGTCAGAATTTTTGCCTTTTTCTCGACCAGCAATGGGCGTGGAGTAACTCGCTTCAGTTAAAGAGGTTCTCGAATCCGGTTg                                                                                                                                                        <  1:394430/150‑1 (MQ=255)
aaaaTATACACAACACTCTTTATAGCAAATATAAGGGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTg                                                                                                                                                        <  1:1264462/150‑1 (MQ=255)
 aaaTATACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTgg                                                                                                                                                      <  1:1169606/150‑1 (MQ=255)
  aaTATACACAACACTCTATATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTATCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGa                                                                                                                                                     <  1:1603107/150‑1 (MQ=255)
   atatACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGAt                                                                                                                                                    <  1:1373559/150‑1 (MQ=255)
   atatACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGAt                                                                                                                                                    <  2:1269826/150‑1 (MQ=255)
   atatACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGAt                                                                                                                                                    >  1:450514/1‑150 (MQ=255)
        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:344603/150‑4 (MQ=255)
               ctctTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggatt                                                                                                                                        <  2:56445/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:441162/1‑150 (MQ=255)
                      taGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAAt                                                                                                                                 <  2:800903/150‑1 (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                                                                                                           >  2:997408/1‑150 (MQ=255)
                                            aTTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCttttt                                                                                                           >  1:521663/1‑150 (MQ=255)
                                             ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                          <  2:1126818/150‑1 (MQ=255)
                                             ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                          <  2:1190380/150‑1 (MQ=255)
                                             ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                          <  2:1065120/150‑1 (MQ=255)
                                             ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                          <  1:947686/150‑1 (MQ=255)
                                               cAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                                                                        >  2:1715546/1‑150 (MQ=255)
                                               cAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                                                                        >  1:1392708/1‑150 (MQ=255)
                                                    gCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTg                                                                                                   >  1:1691837/1‑150 (MQ=255)
                                                         aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                                                              >  2:873725/1‑150 (MQ=255)
                                                         aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                                                              >  1:1480769/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:1502977/1‑150 (MQ=255)
                                                                                                            ggaggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCAc                                           <  1:1036466/150‑1 (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            >  1:1202631/1‑150 (MQ=255)
                                                                                                                                           gAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGAt            >  2:1471126/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      <  2:1392708/150‑1 (MQ=255)
                                                                                                                                                 ggTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTg      <  1:629996/150‑1 (MQ=255)
                                                                                                                                                  gTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGATGGGAAATCAGCATGCCATCGCGGTCAGATCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGa     >  1:895069/1‑150 (MQ=255)
                                                                                                                                                     ggATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTa  <  2:521663/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAAATATACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTA  >  NC_000913/2365861‑2366159

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: