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

AAAATATACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTA  >  NC_000913/2365861‑2366159
                                                                                                                                                      |                                                                                                                                                     
aaaaTATACACAACACTCTTTATAGCAAATATTACTGGACAGGTATTCAATGTCGGAAGGAAAAGCTATGTCAGAATTTTTGCCTTTTTCTCGACCAGCAATGGGCGTGGAGTAACTCGCTTCAGTTAAAGAGGTTCTCGAATCCGGTTg                                                                                                                                                        <  1:394430/150‑1 (MQ=255)
aaaaTATACACAACACTCTTTATAGCAAATATAAGGGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTg                                                                                                                                                        <  1:1264462/150‑1 (MQ=255)
 aaaTATACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTg‑g                                                                                                                                                      <  1:1169606/150‑1 (MQ=255)
  aaTATACACAACACTCTATATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTATCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑Ga                                                                                                                                                     <  1:1603107/150‑1 (MQ=255)
   atatACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GAt                                                                                                                                                    <  1:1373559/150‑1 (MQ=255)
   atatACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GAt                                                                                                                                                    <  2:1269826/150‑1 (MQ=255)
   atatACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GAt                                                                                                                                                    >  1:450514/1‑150 (MQ=255)
        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:344603/150‑4 (MQ=255)
               ctctTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggatt                                                                                                                                        <  2:56445/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:441162/1‑150 (MQ=255)
                      taGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAAt                                                                                                                                 <  2:800903/150‑1 (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                                                                                                           >  2:997408/1‑150 (MQ=255)
                                            aTTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCttttt                                                                                                           >  1:521663/1‑150 (MQ=255)
                                             ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                          <  2:1126818/150‑1 (MQ=255)
                                             ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                          <  2:1190380/150‑1 (MQ=255)
                                             ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                          <  2:1065120/150‑1 (MQ=255)
                                             ttCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                                                                          <  1:947686/150‑1 (MQ=255)
                                               cAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                                                                        >  2:1715546/1‑150 (MQ=255)
                                               cAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                                                                        >  1:1392708/1‑150 (MQ=255)
                                                    gCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTg                                                                                                   >  1:1691837/1‑150 (MQ=255)
                                                         aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                                                              >  2:873725/1‑150 (MQ=255)
                                                         aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                                                              >  1:1480769/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:1502977/1‑150 (MQ=255)
                                                                                                            ggaggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCAc                                           <  1:1036466/150‑1 (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            >  1:1202631/1‑150 (MQ=255)
                                                                                                                                           gAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGAt            >  2:1471126/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      <  2:1392708/150‑1 (MQ=255)
                                                                                                                                                 ggTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTg      <  1:629996/150‑1 (MQ=255)
                                                                                                                                                  gTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGATGGGAAATCAGCATGCCATCGCGGTCAGATCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGa     >  1:895069/1‑150 (MQ=255)
                                                                                                                                                     g‑gATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTa  <  2:521663/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAAATATACACAACACTCTTTATAGCAAATATAAGTGGACAGGTATTCAATGGCGGAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTA  >  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: ‑