Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 293:1 +G 6.4% intergenic (+38/‑44) thrL → / → thrA thr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132931.G 6.4% 1296.8 / 80.4 283intergenic (+38/‑44)thrL/thrAthr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base . (137/128);  new base G (12/6);  total (149/134)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-01

Only 100 of 286 total aligned reads displayed.

AGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACC  >  NC_000913/144‑441
                                                                                                                                                      |                                                                                                                                                     
  gCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGt                                                                                                                                                     >  2:1045175/1‑150 (MQ=255)
     tAGCGCACAGACAGATAAAAATTACAGATTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACCGGTAAGGGTGCGGGCTGACGCGTACAGTAAACAAATAAAAAAGCCCGCACCTGACAGTGCg                                                                                                                                                  <  1:1672712/150‑1 (MQ=255)
      aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCgg                                                                                                                                                 <  1:5235622/150‑1 (MQ=255)
      aGCGCACAGAAATATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCTTACAGTAAACACAGAAAAAAGCCCGCACCTGACAGTGCgg                                                                                                                                                 <  1:5878531/150‑1 (MQ=255)
       gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                                                                                <  2:5932608/150‑1 (MQ=255)
                cagaTAAACATTACAGAGTACACCACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtttttttt                                                                                                                                      <  2:5086576/150‑1 (MQ=255)
                 agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttt                                                                                                                                     >  1:4289894/1‑150 (MQ=255)
                 agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACACAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttt                                                                                                                                     <  1:728572/150‑1 (MQ=255)
                  gaTAAAAATTACAGCGTCCACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCCCAGGTAACGGTGCGGGCTGACGCGTACCGGAAACACACAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTc                                                                                                                                    <  1:193706/150‑1 (MQ=255)
                  gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTc                                                                                                                                    <  1:1011453/150‑1 (MQ=255)
                   aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCg                                                                                                                                   >  2:2990039/1‑150 (MQ=255)
                    tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGa                                                                                                                                  >  2:222271/1‑150 (MQ=255)
                      aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGAcc                                                                                                                                >  1:1243933/1‑150 (MQ=255)
                           tACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAATAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAgg                                                                                                                           >  1:2309024/1‑150 (MQ=255)
                                    cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGt                                                                                                                  <  1:4379113/150‑1 (MQ=255)
                                     acaacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTa                                                                                                                 >  2:6626959/1‑150 (MQ=255)
                                       aacaTCAATTAAACTCATTAGCAACACCATTACCACCACAATCAACATTACCACAGGTAAGGGTGCGGGCTGACGCGTACAGGAAACACATAAAAAAGCCCGAACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                                                                               <  1:6248707/150‑1 (MQ=255)
                                          aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAATGCCCGCACCTGACAGTGCGGGCTTTTTTTTTAGTACAAAGGTAACGAGGTaacaac                                                                                                            >  2:2848004/1‑150 (MQ=255)
                                           tCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAAAGGAGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAAAGAGGTAACAAcc                                                                                                           >  2:2236407/1‑150 (MQ=255)
                                            ccATGAAACGCATTAGCACCACCATTACCACCACCCTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                          <  2:2498114/150‑1 (MQ=255)
                                            ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                          >  1:6449921/1‑150 (MQ=255)
                                            ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGGAACAACCa                                                                                                          >  2:15750/1‑150 (MQ=255)
                                             cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                         >  1:3669718/1‑150 (MQ=255)
                                             cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                         <  1:3841669/150‑1 (MQ=255)
                                              aTGAAACGCATTAGCACCACCCTGACCACCACCATCACCATGACCACCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACCGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                        <  1:5655495/150‑1 (MQ=255)
                                               tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                       >  1:3282694/1‑150 (MQ=255)
                                                 aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCgc                                                                                                     >  2:5331141/1‑149 (MQ=255)
                                                  aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                                                                                     <  2:3368225/150‑1 (MQ=255)
                                                   aCGCATTAGCACCACCATTACCCCCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                                                                                   <  2:3157268/150‑1 (MQ=255)
                                                      cATTAGCACCCCCATTACCACCACCCTCCCCATTACCCCAGGTAACGGTGCGGGCTGACGCGGACAGGAAACACCGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                                <  2:5352481/150‑1 (MQ=255)
                                                       aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAATTGCCCGCAGCTGACAGTGCGGGCTTTTTTTTTCGACCAAATGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                               >  1:1487937/1‑150 (MQ=255)
                                                       aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                               <  1:3402985/150‑1 (MQ=255)
                                                        ttAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                                                                              <  1:4480164/150‑1 (MQ=255)
                                                         tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGTAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                             <  1:4951542/150‑1 (MQ=255)
                                                          aGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                                                                            >  1:3018042/1‑150 (MQ=255)
                                                             accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                         >  1:913033/1‑150 (MQ=255)
                                                               caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                       >  1:2482163/1‑150 (MQ=255)
                                                               caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGc                                                                                                                                               >  1:6674037/1‑95 (MQ=38)
                                                                accaTTCCCACCACCATCACCATTACCACAGTTAACGGTGCGGGATGACCCTTACAGTAAACACACAAAAAACCCCCAACCTTACATTGCGGGATTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                      <  1:2848004/150‑1 (MQ=255)
                                                                accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                      >  2:744603/1‑150 (MQ=255)
                                                                        caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                              <  1:1864563/150‑1 (MQ=255)
                                                                        caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                              >  1:5916094/1‑150 (MQ=255)
                                                                         accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                                                                             >  2:2170256/1‑150 (MQ=255)
                                                                           caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                           <  2:913033/150‑1 (MQ=255)
                                                                           caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGGACATCAGt                                                                           >  1:1562459/1‑150 (MQ=255)
                                                                             ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                                                                         >  2:4846576/1‑150 (MQ=255)
                                                                              catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttttc                                                                                                                                   >  1:2116983/1‑89 (MQ=255)
                                                                              catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                                                                        >  2:5825351/1‑150 (MQ=255)
                                                                                tcaccatTACCCCAGGTAACGGTGCGGGCTGACGCGGACAGGACCCACAGAAAAAAGCCCGCACCTGACCGTGCGGGCTTTTTTGTTCGACCAAAGGTACCGAGGGACCACCCATGCGAGGGTTGAAGTTCGGCGGTCCATCAGTGGCac                                                                      >  2:391265/1‑149 (MQ=255)
                                                                                tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa                                                                      >  1:6130442/1‑150 (MQ=255)
                                                                                  accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt                                                                    >  1:637050/1‑150 (MQ=255)
                                                                                      ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCCAAGGTAACGAGGTAACAACCATGCGAGTGGTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg                                                                >  2:6209750/1‑150 (MQ=255)
                                                                                          cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACCGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACg                                                            >  1:2629846/1‑150 (MQ=255)
                                                                                           acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGAc                                                                                                                                 >  1:571853/1‑80 (MQ=255)
                                                                                               gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCaaa                                                                                                                             <  2:1413701/80‑1 (MQ=255)
                                                                                               gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                       >  2:2958611/1‑150 (MQ=255)
                                                                                               gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCAGACAGTGCGGGCATTTTTTTTCGACCAAAGGTAACGAGGTAACCACCAAGCGCGAGTTGAAGTTCGGCGGTACATCAGTGGCAAATTCCGAACGtttta                                                       >  1:3470938/1‑149 (MQ=255)
                                                                                                    ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACCGTGCGGGCTTTTTTTTTCGACCCAAGGTAACGAGGTAACAACCATGCGCGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                  >  1:5767510/1‑150 (MQ=255)
                                                                                                     gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgtg                                                 >  1:5479770/1‑150 (MQ=255)
                                                                                                        cGGGCTGACGCGTACAGGAAACCCCGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCACAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTg                                              <  1:1779954/150‑1 (MQ=255)
                                                                                                         gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCCAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCCAATGCAGAACGTTTTCTGCGTGTTGc                                             >  1:2148338/1‑150 (MQ=255)
                                                                                                         gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCGAATGCAGAACGTTTTCTGCGTGTTGc                                             <  2:2724528/150‑1 (MQ=255)
                                                                                                         gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                             <  1:5825379/150‑1 (MQ=255)
                                                                                                             tGACGCGTACAGGAAACACAGAAAATAGCCCGCACCTGACAGTTCGGGCTTTTTTTTTGGACCAAATGTAACGAGGTAAAAACCATGGGAGTGTTGTAGTTGGGGGGTACATCAGTGGGAAATGCAGAAAGTTTTGTGGGTGTTGGCgtt                                         >  1:6097432/1‑148 (MQ=255)
                                                                                                               aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGCGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAATCAGAACGTTTTCTGCGTGTTGCCGatat                                       >  1:3738145/1‑150 (MQ=255)
                                                                                                                cgcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATAtt                                      <  1:2337356/150‑1 (MQ=255)
                                                                                                                  cgTACAGGAAACACAGAAAAAATCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAAAAACCAAGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAAATCAGAAAGTTTTGTGCGTGTTTCCGATAttat                                    >  1:6511984/1‑148 (MQ=255)
                                                                                                                   gTACAGGAAACCCAGACCACAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAACGGTAACGAGGTAACACCCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAACTGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                   <  2:5767510/150‑1 (MQ=255)
                                                                                                                   gTACAGGAAACACAGAAAATAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGGCCAAAGGTAAACAGGTAACAACCATGCGAGGGTTTAAGTTCGGCGGTACATCCGTTGCAAAAGCAGAAAGTTTTGTGCGGTTTTCCGAGATTCTg                                   >  1:4991310/1‑150 (MQ=255)
                                                                                                                   gTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                   <  1:1198229/150‑1 (MQ=255)
                                                                                                                    tACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                                  <  1:4404085/150‑1 (MQ=255)
                                                                                                                     aCAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                 <  2:1839639/150‑1 (MQ=255)
                                                                                                                     aCAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                 >  2:6090087/1‑150 (MQ=255)
                                                                                                                     aCAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                 <  1:4634701/150‑1 (MQ=255)
                                                                                                                     aCAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                 <  1:6083390/150‑1 (MQ=255)
                                                                                                                       aGGAAACACAGCAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                               <  2:4426181/150‑1 (MQ=255)
                                                                                                                         gAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTATGGAAAGc                             >  1:4132860/1‑150 (MQ=255)
                                                                                                                          aaaCACAGAAAAAAGCCCGCACCTGACAGTGCGGGATTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCa                            <  2:5686481/150‑1 (MQ=255)
                                                                                                                           aaCACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATAGTCTGGTAAGcga                           >  2:6039157/1‑148 (MQ=255)
                                                                                                                            acacAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAt                          <  1:1422153/150‑1 (MQ=255)
                                                                                                                             cacaGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATg                         <  2:40039/150‑1 (MQ=255)
                                                                                                                              acaGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCCTGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCacggc                        >  2:4757632/1‑146 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                      >  2:3608490/1‑150 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                      <  1:2207011/150‑1 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAGGCca                      >  2:6535678/1‑150 (MQ=255)
                                                                                                                                  aaaaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                    >  2:4875209/1‑150 (MQ=255)
                                                                                                                                   aaaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggc                   >  2:5319588/1‑150 (MQ=255)
                                                                                                                                      aaGCCCGCTCCTGACAGTGGGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg                <  2:4366796/150‑1 (MQ=255)
                                                                                                                                      aaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGGTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATGCTGGAAAGCAATGCcaggccg                 >  2:1227039/1‑148 (MQ=255)
                                                                                                                                       aGCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAAAGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                   >  2:4548903/1‑130 (MQ=255)
                                                                                                                                        gCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg              <  1:2004780/150‑1 (MQ=255)
                                                                                                                                        gCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATCTTCTGGAAAGCAATGCCAggcagggg              >  2:1596963/1‑150 (MQ=255)
                                                                                                                                         cccGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                  >  2:1758511/1‑113 (MQ=255)
                                                                                                                                          ccGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca            >  2:1121851/1‑150 (MQ=255)
                                                                                                                                          ccGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGAGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca            >  2:6391427/1‑150 (MQ=255)
                                                                                                                                             cACCTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGTAACAAACATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTGCCGATATTATGGAAAGCAATGCCAGGCAGGGGCAGGt         >  1:51843/1‑150 (MQ=255)
                                                                                                                                                cTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc      >  2:2101567/1‑150 (MQ=255)
                                                                                                                                                cTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGCTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc      <  2:1301947/150‑1 (MQ=255)
                                                                                                                                                cTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAAATCAGAACGTTTTATGCGTGTTGCCGATATTATGGAAAGCAATGCCAGGCAGGGGCAGGTGGc      >  1:4021297/1‑150 (MQ=255)
                                                                                                                                                    cAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAcc  <  2:3782536/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACC  >  NC_000913/144‑441

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