Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 293:1 +G 5.6% 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 5.6% 1406.8 / 74.9 305intergenic (+38/‑44)thrL/thrAthr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base . (142/146);  new base G (11/6);  total (153/152)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

Only 100 of 312 total aligned reads displayed.

AGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCG  >  NC_000913/144‑442
                                                                                                                                                       |                                                                                                                                                      
aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACa                                                                                                                                                         >  1:1097979/1‑150 (MQ=255)
 ggCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑g                                                                                                                                                       >  1:5270650/1‑150 (MQ=255)
   cATAGCGCACAGAAATATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACTGTTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTg                                                                                                                                                     <  1:3756054/150‑1 (MQ=255)
    aTAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGc                                                                                                                                                    <  2:2873703/150‑1 (MQ=255)
    aTAGCGCACAGAAATATAAAAATTCCAGATTACACAACATCCATGAAACGCATTCGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGc                                                                                                                                                    <  1:3301309/150‑1 (MQ=255)
     tAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCg                                                                                                                                                   >  2:1317941/1‑150 (MQ=255)
       gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                                 <  2:2488191/150‑1 (MQ=255)
        cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGc                                                                                                                                                <  2:5989434/150‑1 (MQ=255)
         gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑t                                                                                                                                              >  1:3496774/1‑150 (MQ=255)
         gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑t                                                                                                                                              >  2:3071901/1‑150 (MQ=255)
               acagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttt                                                                                                                                        >  1:4984953/1‑150 (MQ=255)
                 agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACCCAT‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttttt                                                                                                                                      <  2:3123552/150‑1 (MQ=255)
                    tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCCCCACCATCCCCATTACCACCGGTAACGGTGCGGGCTGACGCGTAAAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGa                                                                                                                                   <  1:1031582/150‑1 (MQ=255)
                    tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGa                                                                                                                                   >  2:334671/1‑150 (MQ=255)
                                    cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACTGGAAACACAG‑AAATTTGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGTCCAATGGTAACGAGGt                                                                                                                   >  1:3182592/1‑150 (MQ=255)
                                     acaacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTa                                                                                                                  >  2:1095487/1‑150 (MQ=255)
                                       aacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                                                                                <  1:5571987/150‑1 (MQ=255)
                                       aacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                                                                                >  2:1904828/1‑150 (MQ=255)
                                         caTCCATGAATCGCATTAGCACCACCATTACCACCACAATCACCATTACCACAGTTATCGGTGCGGGCTGACGCGTACATTAAAAACAG‑AAAAAAGCCCGCACCTGACA‑TTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                                                                              <  1:4748322/150‑1 (MQ=255)
                                         caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                                                                              <  1:1426167/150‑1 (MQ=255)
                                         caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                                                                              <  1:442134/150‑1 (MQ=255)
                                             cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                          >  1:1478216/1‑150 (MQ=255)
                                              aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                         >  1:268318/1‑150 (MQ=255)
                                              aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                         >  1:1864407/1‑150 (MQ=255)
                                               tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                        <  2:2687219/150‑1 (MQ=255)
                                                 aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGAACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                                                                                      <  1:711836/150‑1 (MQ=255)
                                                  aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                                                                                     >  2:1700363/1‑150 (MQ=255)
                                                  aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                                                                                     >  2:4371129/1‑150 (MQ=255)
                                                   aCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                                                                                    <  2:6433495/150‑1 (MQ=255)
                                                       aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAATGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                                >  2:1998285/1‑150 (MQ=255)
                                                       aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                                >  2:4190690/1‑150 (MQ=255)
                                                       aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                                >  2:5578254/1‑150 (MQ=255)
                                                        ttAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                                                                               <  1:5321388/150‑1 (MQ=255)
                                                         tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                              <  2:2085677/150‑1 (MQ=255)
                                                         tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                              <  2:5346533/150‑1 (MQ=255)
                                                              ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                         >  1:183115/1‑150 (MQ=255)
                                                               caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTCACAACCATGCGAGTGTTGAAGGTcgg                                                                                        >  2:4676801/1‑150 (MQ=255)
                                                                accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                       <  1:6687343/150‑1 (MQ=255)
                                                                    ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                   <  1:3613210/150‑1 (MQ=255)
                                                                               atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                                                                        >  2:3758342/1‑150 (MQ=255)
                                                                               atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGTACATCAGTGGcc                                                                        >  1:6547831/1‑149 (MQ=255)
                                                                                tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTCGCCCCAAGGTAACGAGGTAACAACCCTGCGGGGGTTGAAGTTCGGCGGTACAACAGTGGCa                                                                        >  1:3049377/1‑150 (MQ=255)
                                                                                   ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTACCGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg                                                                    >  2:1662194/1‑150 (MQ=255)
                                                                                    catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                   >  2:3614339/1‑150 (MQ=255)
                                                                                     atTACCACAGGTACCGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCCCCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGGACAGCCGTGGCAAATGCa                                                                  >  1:4529842/1‑150 (MQ=255)
                                                                                     atTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa                                                                  <  1:549435/150‑1 (MQ=255)
                                                                                     atTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa                                                                  <  2:3276340/150‑1 (MQ=255)
                                                                                      ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAAGCCCGCACATGACA‑TTGCGGGG‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACAATGCGAGTTTTGAGGTGCGGCGGTCCAGCAGTGGCAAATGCa                                                                  >  2:3683978/1‑150 (MQ=255)
                                                                                      ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg                                                                 <  1:582328/150‑1 (MQ=255)
                                                                                        aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa                                                               >  2:2518410/1‑150 (MQ=255)
                                                                                         ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTCCACCAACGGTAACGAGGTACCAACCATCCGAGGGTTGAAGTTCGGCGGTACATCATTGGCAATTGCAGaa                                                               >  2:3988615/1‑150 (MQ=255)
                                                                                           acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                                >  1:6398138/1‑114 (MQ=255)
                                                                                           acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGt                                                            >  1:1242870/1‑150 (MQ=255)
                                                                                               gTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCTGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                        >  1:1934238/1‑150 (MQ=255)
                                                                                               gTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCaa                                                                                                                               <  2:714032/79‑1 (MQ=255)
                                                                                                   cGGTGCGGGCTGACGCGTACAGGAAACACAC‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGTACGTTTTCTGCg                                                    <  1:6180227/150‑1 (MQ=255)
                                                                                                    ggTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                   <  1:2479478/150‑1 (MQ=255)
                                                                                                      tGCGGGCTGACGCGGCCAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgtgt                                                 <  2:3413163/150‑1 (MQ=255)
                                                                                                      ggCGGGCTGTCGCGTACTGGAAACACAC‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGTGGTAACAACCATGCGAGTGTTTAAGTTCGGCGGTACATCAGTGGCAAATGCAGATCGTTTTCTGCgtgt                                                 <  2:5134568/149‑1 (MQ=255)
                                                                                                          ggCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                             >  2:2960851/1‑150 (MQ=255)
                                                                                                           gCTGACGCGTACAGGAAACACAG‑AATACAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCACAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCg                                            <  2:4283055/150‑1 (MQ=255)
                                                                                                               aCGCGTACAGGAAACACAG‑AAAAATGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGTCCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGatat                                        >  1:1542381/1‑150 (MQ=255)
                                                                                                                 gcgTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTc                                      >  2:153821/1‑150 (MQ=255)
                                                                                                                  cgTACAGGAAACACAG‑AAAAAAGCCCGCACCTTACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAAAAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCTAATTCAGAACGTTTTGTGCGTGTTTGCGATATTCt                                     >  1:2654417/1‑150 (MQ=255)
                                                                                                                  cgTACAGCAAACACAG‑AACAATGCCCCCAGCTGACG‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCt                                     <  1:2942315/150‑1 (MQ=255)
                                                                                                                   gTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGttt                                                          >  2:2198839/1‑128 (MQ=255)
                                                                                                                   gTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGttt                                                          <  1:2198288/128‑1 (MQ=255)
                                                                                                                   gTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                    <  1:3946068/150‑1 (MQ=255)
                                                                                                                   gTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                    <  1:1982566/150‑1 (MQ=255)
                                                                                                                     aCAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                  <  1:5978681/150‑1 (MQ=255)
                                                                                                                      cAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaa                                 <  2:1591453/150‑1 (MQ=255)
                                                                                                                       aGGAAACACAG‑AAAAATGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGTAGTTCGGCGGTACAGCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                                >  2:6203349/1‑150 (MQ=255)
                                                                                                                        ggAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGGTTTCTGCGTGTTGCCGATATTCTGGAAAg                               >  2:1019834/1‑150 (MQ=255)
                                                                                                                           aaCACCG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCaa                            <  2:4950408/150‑1 (MQ=255)
                                                                                                                            acacAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAt                           >  2:530480/1‑150 (MQ=255)
                                                                                                                              acaG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                         <  1:2063963/150‑1 (MQ=255)
                                                                                                                              acaG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                         <  1:4340458/150‑1 (MQ=255)
                                                                                                                              acaG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                         >  2:6556197/1‑150 (MQ=255)
                                                                                                                              acaG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                         >  2:55395/1‑150 (MQ=255)
                                                                                                                               gaaaaaaaaaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaa                                 <  1:3157145/138‑1 (MQ=38)
                                                                                                                                aG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                       >  2:5873992/1‑150 (MQ=255)
                                                                                                                                aG‑AAAAAAGCCCGAACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                       <  1:495498/150‑1 (MQ=255)
                                                                                                                                   aaaaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                     >  2:1730430/1‑150 (MQ=255)
                                                                                                                                     aaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggca                   <  2:54061/150‑1 (MQ=255)
                                                                                                                                      aaaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcag                  >  2:6040029/1‑150 (MQ=255)
                                                                                                                                       aaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg                 <  2:1480886/150‑1 (MQ=255)
                                                                                                                                       aaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg                 <  1:5250861/150‑1 (MQ=255)
                                                                                                                                       aaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg                 >  2:3302720/1‑150 (MQ=255)
                                                                                                                                        aGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCTGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggg                >  2:86178/1‑150 (MQ=255)
                                                                                                                                         gCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg               >  2:2385618/1‑150 (MQ=255)
                                                                                                                                         gCCCGCA‑CTGTCAGGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggg                >  2:3238587/1‑150 (MQ=255)
                                                                                                                                           ccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGTAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTATGGAAAGCAATGCCAggcaggggca             >  1:2801094/1‑150 (MQ=255)
                                                                                                                                           ccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca             >  2:1865063/1‑150 (MQ=255)
                                                                                                                                           ccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAAATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTTCCGATATTCTGGAAAGCAATGCCAggcaggggca             >  1:5102858/1‑150 (MQ=255)
                                                                                                                                            cGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggccg            >  1:883077/1‑148 (MQ=255)
                                                                                                                                            cGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGACGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACCTTATCTGCGTGTTGCCGATATTATGGAAAGCAATGCCAggcaggggcag            >  2:4995184/1‑150 (MQ=255)
                                                                                                                                              cACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGt          <  1:3090262/150‑1 (MQ=255)
                                                                                                                                               aCCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg         >  2:3911697/1‑150 (MQ=255)
                                                                                                                                               aCCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg         <  2:2931537/150‑1 (MQ=255)
                                                                                                                                               aCCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg         <  1:4217544/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
AGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCG  >  NC_000913/144‑442

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