Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 291 A→T 5.5% intergenic (+36/‑46) 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_0009132910AT 5.5% 1111.3 / 39.1 312intergenic (+36/‑46)thrL/thrAthr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (148/145);  new base T (11/6);  total (159/153)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

Only 100 of 314 total aligned reads displayed.

ATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCA  >  NC_000913/142‑439
                                                                                                                                                     |                                                                                                                                                       
ataGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGa                                                                                                                                                         >  2:2479478/1‑150 (MQ=255)
ataGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGa                                                                                                                                                         >  1:3294512/1‑150 (MQ=255)
 taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGAc                                                                                                                                                        >  1:579041/1‑150 (MQ=255)
  aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACa                                                                                                                                                       >  1:1097979/1‑150 (MQ=255)
     cATAGCGCACAGAAATATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACTGTTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTg                                                                                                                                                   <  1:3756054/150‑1 (MQ=255)
      aTAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGc                                                                                                                                                  <  2:2873703/150‑1 (MQ=255)
          cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGc                                                                                                                                              <  2:5989434/150‑1 (MQ=255)
          cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGc                                                                                                                                              <  2:587198/150‑1 (MQ=255)
           gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑t                                                                                                                                            >  1:3496774/1‑150 (MQ=255)
                 acagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttt                                                                                                                                      >  1:4984953/1‑150 (MQ=255)
                  cataTAAAAATTACAGATTACACAACATCCATGAAACGACTTAGCACCACCATGACCCCCACCACCACCTTTACCACCGGTACGGGTGGGGGCTTACGCTAACAGTAAACACAGAAAAAAGCCCGCAGTTTCTG‑GTGCGGGC‑ttttttct                                                                                                                                     <  1:2959746/147‑3 (MQ=255)
                     aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCg                                                                                                                                 >  2:1051503/1‑150 (MQ=255)
                      tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCCCCACCATCCCCATTACCACCGGTAACGGTGCGGGCTGACGCGTAAAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGa                                                                                                                                <  1:1031582/150‑1 (MQ=255)
                      tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGa                                                                                                                                >  2:3629324/1‑150 (MQ=255)
                      tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGa                                                                                                                                >  2:334671/1‑150 (MQ=255)
                         aaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCa                                                                                                                             >  2:3172149/1‑150 (MQ=255)
                                 gagTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCCAAGGTAAc                                                                                                                     >  2:2006148/1‑150 (MQ=255)
                                      cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACTGGAAACACAGAAATTTGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGTCCAATGGTAACGAGGt                                                                                                                >  1:3182592/1‑150 (MQ=255)
                                      cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGt                                                                                                                <  2:3755686/150‑1 (MQ=255)
                                        caacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTaa                                                                                                              <  1:3893695/150‑1 (MQ=255)
                                         aacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTaac                                                                                                             >  2:1904828/1‑150 (MQ=255)
                                         aacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTaac                                                                                                             <  1:5952975/150‑1 (MQ=255)
                                         aacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTaac                                                                                                             <  1:5571987/150‑1 (MQ=255)
                                           caTCCATGAATCGCATTAGCACCACCATTACCACCACAATCACCATTACCACAGTTATCGGTGCGGGCTGACGCGTACATTAAAAACAGAAAAAAGCCCGCACCTGACA‑TTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTaacaa                                                                                                           <  1:4748322/150‑1 (MQ=255)
                                           caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTaacaa                                                                                                           <  1:1426167/150‑1 (MQ=255)
                                            aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTc                                                                                                                                  >  2:497261/1‑126 (MQ=255)
                                               cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                       >  1:1478216/1‑150 (MQ=255)
                                                aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                      >  1:268318/1‑150 (MQ=255)
                                                aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                      >  1:1864407/1‑150 (MQ=255)
                                                    aaCGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                                                                                  <  1:1935915/150‑1 (MQ=255)
                                                    aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                                                                                  >  2:1700363/1‑150 (MQ=255)
                                                    aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                                                                                  >  2:4371129/1‑150 (MQ=255)
                                                     aCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                                                                                 <  2:6433495/150‑1 (MQ=255)
                                                         aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                             >  2:5578254/1‑150 (MQ=255)
                                                         aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                             >  2:4190690/1‑150 (MQ=255)
                                                          ttAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                                                                            <  1:5321388/150‑1 (MQ=255)
                                                           tcgCACCACCAGTCCCACCACCATCACCATTACCACCGGTACCGGTGCGGGCGGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACT‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                           <  1:4416654/148‑1 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                           <  2:2085677/150‑1 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                           <  2:5346533/150‑1 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                           <  1:675644/150‑1 (MQ=255)
                                                                ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                      >  1:183115/1‑150 (MQ=255)
                                                                 caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTCACAACCATGCGAGTGTTGAAGGTcgg                                                                                     >  2:4676801/1‑150 (MQ=255)
                                                                 caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                     <  1:6514977/150‑1 (MQ=255)
                                                                      ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                <  1:3613210/150‑1 (MQ=255)
                                                                          caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGGTCGGCGGTACATc                                                                            >  2:4882522/1‑150 (MQ=255)
                                                                            ccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGACGTTCGGCGGTACATCAg                                                                          >  2:3793177/1‑150 (MQ=255)
                                                                               ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                                                                       <  1:6035251/150‑1 (MQ=255)
                                                                                catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                                                                      <  2:295046/150‑1 (MQ=255)
                                                                                  tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTT‑TTCGCCCCAAGGTAACGAGGTAACAACCCTGCGGGGGTTGAAGTTCGGCGGTACAACAGTGGCa                                                                     >  1:3049377/1‑150 (MQ=255)
                                                                                     ccatTCCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACT‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg                                                                 <  2:3258965/150‑1 (MQ=255)
                                                                                     ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTACCGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg                                                                 >  2:1662194/1‑150 (MQ=255)
                                                                                       ataacaacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCAGACA‑GTGCGGGC‑TTTTTTTCTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                >  2:306416/7‑150 (MQ=255)
                                                                                       atTACCACAGGTACCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGGACAGCCGTGGCAAATGCa                                                               >  1:4529842/1‑150 (MQ=255)
                                                                                       atTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa                                                               <  1:549435/150‑1 (MQ=255)
                                                                                           ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACC‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc                                                           >  1:4309925/1‑150 (MQ=255)
                                                                                           ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTT‑TTCCACCAACGGTAACGAGGTACCAACCATCCGAGGGTTGAAGTTCGGCGGTACATCATTGGCAATTGCAGaa                                                            >  2:3988615/1‑150 (MQ=255)
                                                                                             acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGt                                                         >  1:1242870/1‑150 (MQ=255)
                                                                                                 gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCaa                                                                                                                            <  2:714032/79‑1 (MQ=255)
                                                                                                     cGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCg                                                 >  1:2819565/1‑150 (MQ=255)
                                                                                                      ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                <  1:2479478/150‑1 (MQ=255)
                                                                                                        tGCGGGCTGACGCGGCCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgtgt                                              <  2:3413163/150‑1 (MQ=255)
                                                                                                        ggCGGGCTGTCGCGTACTGGAAACACACAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGTGGTAACAACCATGCGAGTGTTTAAGTTCGGCGGTACATCAGTGGCAAATGCAGATCGTTTTCTGCgtgt                                              <  2:5134568/149‑1 (MQ=255)
                                                                                                           gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                           <  1:5160014/150‑1 (MQ=255)
                                                                                                             gCTGACGCGTACAGGAAACACAGAATACAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCACAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCg                                         <  2:4283055/150‑1 (MQ=255)
                                                                                                              cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGa                                        <  1:952410/150‑1 (MQ=255)
                                                                                                                 aCGCGTACAGGAAACACAGAAAAATGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGTCCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGatat                                     >  1:1542381/1‑150 (MQ=255)
                                                                                                                   gcgTACAGGAAACACAGAAAAAAGCCCGCCCCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCCAAGGTACCGCGGGAACCACCATGCGAGGGTTGAAGTTCGGCGGGACCACAGGGGCAAATGCCGGCCGGTTTCGGCGTGTTGCCGaacgtc                                   >  1:5256998/1‑145 (MQ=255)
                                                                                                                   gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTc                                   >  2:3105255/1‑150 (MQ=255)
                                                                                                                   gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTc                                   >  2:153821/1‑150 (MQ=255)
                                                                                                                    cgTACAGGAAACACAGAAAAAAGCCCGCACCTTACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAAAAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCTAATTCAGAACGTTTTGTGCGTGTTTGCGATATTCt                                  >  1:2654417/1‑150 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGttt                                                       <  1:2198288/128‑1 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGttt                                                       >  2:2198839/1‑128 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 <  2:2516006/150‑1 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 <  1:1982566/150‑1 (MQ=255)
                                                                                                                       aCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCACAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                               <  1:6276865/150‑1 (MQ=255)
                                                                                                                        cAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaa                              <  2:1591453/150‑1 (MQ=255)
                                                                                                                         aGGAAACACAGAAAAATGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGTAGTTCGGCGGTACAGCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                             >  2:6203349/1‑150 (MQ=255)
                                                                                                                              acacAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAt                        >  2:530480/1‑150 (MQ=255)
                                                                                                                                gaaaaaaaaaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaa                              <  1:3157145/138‑1 (MQ=38)
                                                                                                                                acaGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                      >  2:55395/1‑150 (MQ=255)
                                                                                                                                acaGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                      <  1:2063963/150‑1 (MQ=255)
                                                                                                                                acaGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                      >  2:6556197/1‑150 (MQ=255)
                                                                                                                                 aaaaaaaaaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 >  2:3565358/4‑138 (MQ=38)
                                                                                                                                  aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                    <  2:1242870/150‑1 (MQ=255)
                                                                                                                                  aGAAAAAAGCCCGAACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                    <  1:495498/150‑1 (MQ=255)
                                                                                                                                   gAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcag                   <  1:6663898/150‑1 (MQ=255)
                                                                                                                                    aaaaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                  >  2:1730430/1‑150 (MQ=255)
                                                                                                                                    aaaaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                  <  1:4190690/150‑1 (MQ=255)
                                                                                                                                      aaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggca                <  2:54061/150‑1 (MQ=255)
                                                                                                                                       aaaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcag               >  2:6040029/1‑150 (MQ=255)
                                                                                                                                        aaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg              <  2:1480886/150‑1 (MQ=255)
                                                                                                                                        aaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg              >  2:3302720/1‑150 (MQ=255)
                                                                                                                                         aGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCTGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggg             >  2:86178/1‑150 (MQ=255)
                                                                                                                                          gCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg            >  2:2385618/1‑150 (MQ=255)
                                                                                                                                            ccGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca          >  2:1865063/1‑150 (MQ=255)
                                                                                                                                            ccGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAAATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTTCCGATATTCTGGAAAGCAATGCCAggcaggggca          >  1:5102858/1‑150 (MQ=255)
                                                                                                                                             cGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag         >  1:882914/1‑150 (MQ=255)
                                                                                                                                               cACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGt       <  1:3090262/150‑1 (MQ=255)
                                                                                                                                                aCCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg      <  2:2931537/150‑1 (MQ=255)
                                                                                                                                                 ccTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg     <  2:268318/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
ATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTT‑TTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCA  >  NC_000913/142‑439

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: ‑