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.

AGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACC  >  NC_000913/144‑441
                                                                                                                                                      |                                                                                                                                                     
  gCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑Gt                                                                                                                                                     >  2:1045175/1‑150 (MQ=255)
     tAGCGCACAGACAGATAAAAATTACAGATTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACCGGTAAGGGTGCGGGCTGACGCGTACAGTAAACAAATAAAAAAGCCCGCACCTGACA‑GTGCg                                                                                                                                                  <  1:1672712/150‑1 (MQ=255)
      aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                 <  1:5235622/150‑1 (MQ=255)
      aGCGCACAGAAATATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCTTACAGTAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                 <  1:5878531/150‑1 (MQ=255)
       gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                                <  2:5932608/150‑1 (MQ=255)
                cagaTAAACATTACAGAGTACACCACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑tttttttt                                                                                                                                      <  2:5086576/150‑1 (MQ=255)
                 agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttttt                                                                                                                                     >  1:4289894/1‑150 (MQ=255)
                 agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACACAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttttt                                                                                                                                     <  1:728572/150‑1 (MQ=255)
                  gaTAAAAATTACAGCGTCCACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCCCAGGTAACGGTGCGGGCTGACGCGTACCGGAAACACACAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTc                                                                                                                                    <  1:193706/150‑1 (MQ=255)
                  gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTc                                                                                                                                    <  1:1011453/150‑1 (MQ=255)
                   aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCg                                                                                                                                   >  2:2990039/1‑150 (MQ=255)
                    tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGa                                                                                                                                  >  2:222271/1‑150 (MQ=255)
                      aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGAcc                                                                                                                                >  1:1243933/1‑150 (MQ=255)
                           tACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAATAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAgg                                                                                                                           >  1:2309024/1‑150 (MQ=255)
                                    cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGt                                                                                                                  <  1:4379113/150‑1 (MQ=255)
                                     acaacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTa                                                                                                                 >  2:6626959/1‑150 (MQ=255)
                                       aacaTCAATTAAACTCATTAGCAACACCATTACCACCACAATCAACATTACCACAGGTAAGGGTGCGGGCTGACGCGTACAGGAAACACATAAAAAAGCCCGAACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                                                                               <  1:6248707/150‑1 (MQ=255)
                                          aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAATGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGTACAAAGGTAACGAGGTaacaac                                                                                                            >  2:2848004/1‑150 (MQ=255)
                                           tCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAAAGGAGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAAAGAGGTAACAAcc                                                                                                           >  2:2236407/1‑150 (MQ=255)
                                            ccATGAAACGCATTAGCACCACCATTACCACCACCCTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                          <  2:2498114/150‑1 (MQ=255)
                                            ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                          >  1:6449921/1‑150 (MQ=255)
                                            ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGGAACAACCa                                                                                                          >  2:15750/1‑150 (MQ=255)
                                             cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                         >  1:3669718/1‑150 (MQ=255)
                                             cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                         <  1:3841669/150‑1 (MQ=255)
                                              aTGAAACGCATTAGCACCACCCTGACCACCACCATCACCATGACCACCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACCGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                        <  1:5655495/150‑1 (MQ=255)
                                               tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                       >  1:3282694/1‑150 (MQ=255)
                                                 aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCgc                                                                                                     >  2:5331141/1‑149 (MQ=255)
                                                  aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                                                                                     <  2:3368225/150‑1 (MQ=255)
                                                   aCGCATTAGCACCACCATTACCCCCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                                                                                   <  2:3157268/150‑1 (MQ=255)
                                                      cATTAGCACCCCCATTACCACCACCCTCCCCATTACCCCAGGTAACGGTGCGGGCTGACGCGGACAGGAAACACCGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                                <  2:5352481/150‑1 (MQ=255)
                                                       aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAATTGCCCGCAGCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAATGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                               >  1:1487937/1‑150 (MQ=255)
                                                       aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                               <  1:3402985/150‑1 (MQ=255)
                                                        ttAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                                                                              <  1:4480164/150‑1 (MQ=255)
                                                         tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGTAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                             <  1:4951542/150‑1 (MQ=255)
                                                          aGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                                                                            >  1:3018042/1‑150 (MQ=255)
                                                             accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                         >  1:913033/1‑150 (MQ=255)
                                                               caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                       >  1:2482163/1‑150 (MQ=255)
                                                               caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGc                                                                                                                                               >  1:6674037/1‑95 (MQ=38)
                                                                accaTTCCCACCACCATCACCATTACCACAGTTAACGGTGCGGGATGACCCTTACAGTAAACACACAAAAAACCCCCAACCTTACA‑TTGCGGGA‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                      <  1:2848004/150‑1 (MQ=255)
                                                                accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                      >  2:744603/1‑150 (MQ=255)
                                                                        caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                              <  1:1864563/150‑1 (MQ=255)
                                                                        caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                              >  1:5916094/1‑150 (MQ=255)
                                                                         accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                                                                             >  2:2170256/1‑150 (MQ=255)
                                                                           caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                           <  2:913033/150‑1 (MQ=255)
                                                                           caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGGACATCAGt                                                                           >  1:1562459/1‑150 (MQ=255)
                                                                             ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                                                                         >  2:4846576/1‑150 (MQ=255)
                                                                              catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttttttc                                                                                                                                   >  1:2116983/1‑89 (MQ=255)
                                                                              catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                                                                        >  2:5825351/1‑150 (MQ=255)
                                                                                tcaccatTACCCCAGGTAACGGTGCGGGCTGACGCGGACAGGACCCACAGAAAAAAGCCCGCACCTGACC‑GTGCGGGC‑TTTTTTGTTCGACCAAAGGTACCGAGGGACCACCCATGCGAGGGTTGAAGTTCGGCGGTCCATCAGTGGCac                                                                      >  2:391265/1‑149 (MQ=255)
                                                                                tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa                                                                      >  1:6130442/1‑150 (MQ=255)
                                                                                  accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt                                                                    >  1:637050/1‑150 (MQ=255)
                                                                                      ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCCAAGGTAACGAGGTAACAACCATGCGAGTGGTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg                                                                >  2:6209750/1‑150 (MQ=255)
                                                                                          cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACC‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACg                                                            >  1:2629846/1‑150 (MQ=255)
                                                                                           acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGAc                                                                                                                                 >  1:571853/1‑80 (MQ=255)
                                                                                               gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCaaa                                                                                                                             <  2:1413701/80‑1 (MQ=255)
                                                                                               gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                       >  2:2958611/1‑150 (MQ=255)
                                                                                               gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCAGACA‑GTGCGGGC‑ATTTTTTTTCGACCAAAGGTAACGAGGTAACCACCAAGCGCGAGTTGAAGTTCGGCGGTACATCAGTGGCAAATTCCGAACGtttta                                                       >  1:3470938/1‑149 (MQ=255)
                                                                                                    ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACC‑GTGCGGGC‑TTTTTTTTTCGACCCAAGGTAACGAGGTAACAACCATGCGCGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                  >  1:5767510/1‑150 (MQ=255)
                                                                                                     gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgtg                                                 >  1:5479770/1‑150 (MQ=255)
                                                                                                        cGGGCTGACGCGTACAGGAAACCCCGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCACAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTg                                              <  1:1779954/150‑1 (MQ=255)
                                                                                                         gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCCAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCCAATGCAGAACGTTTTCTGCGTGTTGc                                             >  1:2148338/1‑150 (MQ=255)
                                                                                                         gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCGAATGCAGAACGTTTTCTGCGTGTTGc                                             <  2:2724528/150‑1 (MQ=255)
                                                                                                         gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                             <  1:5825379/150‑1 (MQ=255)
                                                                                                             tGACGCGTACAGGAAACACAGAAAATAGCCCGCACCTGACA‑GTTCGGGC‑TTTTTTTTTGGACCAAATGTAACGAGGTAAAAACCATGGGAGTGTTGTAGTTGGGGGGTACATCAGTGGGAAATGCAGAAAGTTTTGTGGGTGTTGGCgtt                                         >  1:6097432/1‑148 (MQ=255)
                                                                                                               aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGCGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAATCAGAACGTTTTCTGCGTGTTGCCGatat                                       >  1:3738145/1‑150 (MQ=255)
                                                                                                                cgcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATAtt                                      <  1:2337356/150‑1 (MQ=255)
                                                                                                                  cgTACAGGAAACACAGAAAAAATCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAAAAACCAAGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAAATCAGAAAGTTTTGTGCGTGTTTCCGATAttat                                    >  1:6511984/1‑148 (MQ=255)
                                                                                                                   gTACAGGAAACCCAGACCACAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAACGGTAACGAGGTAACACCCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAACTGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                   <  2:5767510/150‑1 (MQ=255)
                                                                                                                   gTACAGGAAACACAGAAAATAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGGCCAAAGGTAAACAGGTAACAACCATGCGAGGGTTTAAGTTCGGCGGTACATCCGTTGCAAAAGCAGAAAGTTTTGTGCGGTTTTCCGAGATTCTg                                   >  1:4991310/1‑150 (MQ=255)
                                                                                                                   gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                   <  1:1198229/150‑1 (MQ=255)
                                                                                                                    tACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                                  <  1:4404085/150‑1 (MQ=255)
                                                                                                                     aCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                 <  2:1839639/150‑1 (MQ=255)
                                                                                                                     aCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                 >  2:6090087/1‑150 (MQ=255)
                                                                                                                     aCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                 <  1:4634701/150‑1 (MQ=255)
                                                                                                                     aCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                 <  1:6083390/150‑1 (MQ=255)
                                                                                                                       aGGAAACACAGCAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                               <  2:4426181/150‑1 (MQ=255)
                                                                                                                         gAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTATGGAAAGc                             >  1:4132860/1‑150 (MQ=255)
                                                                                                                          aaaCACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGA‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCa                            <  2:5686481/150‑1 (MQ=255)
                                                                                                                           aaCACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATAGTCTGGTAAGcga                           >  2:6039157/1‑148 (MQ=255)
                                                                                                                            acacAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAt                          <  1:1422153/150‑1 (MQ=255)
                                                                                                                             cacaGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATg                         <  2:40039/150‑1 (MQ=255)
                                                                                                                              acaGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCCTGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCacggc                        >  2:4757632/1‑146 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                      >  2:3608490/1‑150 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                      <  1:2207011/150‑1 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAGGCca                      >  2:6535678/1‑150 (MQ=255)
                                                                                                                                  aaaaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                    >  2:4875209/1‑150 (MQ=255)
                                                                                                                                   aaaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggc                   >  2:5319588/1‑150 (MQ=255)
                                                                                                                                      aaGCCCGCTCCTGACA‑GTGGGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg                <  2:4366796/150‑1 (MQ=255)
                                                                                                                                      aaGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGGTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATGCTGGAAAGCAATGCcaggccg                 >  2:1227039/1‑148 (MQ=255)
                                                                                                                                       aGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAAAGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                   >  2:4548903/1‑130 (MQ=255)
                                                                                                                                        gCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg              <  1:2004780/150‑1 (MQ=255)
                                                                                                                                        gCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATCTTCTGGAAAGCAATGCCAggcagggg              >  2:1596963/1‑150 (MQ=255)
                                                                                                                                         cccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                  >  2:1758511/1‑113 (MQ=255)
                                                                                                                                          ccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca            >  2:1121851/1‑150 (MQ=255)
                                                                                                                                          ccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGAGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca            >  2:6391427/1‑150 (MQ=255)
                                                                                                                                             cACCTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAAACATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTGCCGATATTATGGAAAGCAATGCCAGGCAGGGGCAGGt         >  1:51843/1‑150 (MQ=255)
                                                                                                                                                cTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc      >  2:2101567/1‑150 (MQ=255)
                                                                                                                                                cTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGCTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc      <  2:1301947/150‑1 (MQ=255)
                                                                                                                                                cTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAAATCAGAACGTTTTATGCGTGTTGCCGATATTATGGAAAGCAATGCCAGGCAGGGGCAGGTGGc      >  1:4021297/1‑150 (MQ=255)
                                                                                                                                                    cA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAcc  <  2:3782536/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACC  >  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: ‑