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% 1465.6 / 80.9 320intergenic (+38/‑44)thrL/thrAthr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base . (133/169);  new base G (10/8);  total (143/177)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

Only 100 of 323 total aligned reads displayed.

AGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAC  >  NC_000913/144‑440
                                                                                                                                                      |                                                                                                                                                    
aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACa                                                                                                                                                       <  2:1325086/150‑1 (MQ=255)
      aGCGCACAGACATATAAAAATTACAGACTACACAACACCCATGAAACGCATTAGCACCCCCATTACCACCACCATCACCATTCCCCCCGGTAACGGTGCGGGCTGACGCTTACAGTAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                <  1:1181882/150‑1 (MQ=255)
      aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCTTTAGCACCACCATTACCACCACCATCACCATTACCACCGTTATGGGTGGGGGCTGACGCGTACAGTAAACACATAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                <  1:1758814/150‑1 (MQ=255)
      aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                <  2:4542624/150‑1 (MQ=255)
      aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                <  2:2577854/150‑1 (MQ=255)
        cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGc                                                                                                                                              <  1:6496823/150‑1 (MQ=255)
        cgcACACACAGATAAAAATTACAGAGAACACAACAGCCATGAAACGCATTCGCACCACCATGACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACTGGAAACACACAAAAAAGCCCGCACCTGACA‑GTGCGGGc                                                                                                                                              <  1:1283126/150‑1 (MQ=255)
               acagaTAATAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTCCCACCACCATCACCATTCCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttttttc                                                                                                                                      <  1:6426524/150‑2 (MQ=255)
               acagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttt                                                                                                                                      >  2:6614365/1‑150 (MQ=255)
               acagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttttttc                                                                                                                                      <  2:6215671/150‑2 (MQ=255)
                cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAATAGCCCGCACCTGACA‑GTGCGGGC‑tttttttt                                                                                                                                     >  1:2563341/1‑150 (MQ=255)
                  gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTc                                                                                                                                   <  2:5820214/150‑1 (MQ=255)
                   aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCCTTCCCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCg                                                                                                                                  <  2:2074730/150‑1 (MQ=255)
                   aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCg                                                                                                                                  <  1:4097853/150‑1 (MQ=255)
                    tAAAAATTACAGAGTACACTACCTCCCTGAACCGCATGAGCACCACCAGTACCACCACCAGGACCATTACCCCAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGa                                                                                                                                 <  2:4258161/150‑1 (MQ=255)
                          ttACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAg                                                                                                                           >  1:1413030/1‑150 (MQ=255)
                           tACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAgg                                                                                                                          >  1:2412930/1‑150 (MQ=255)
                              agagTACACAAAATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTaa                                                                                                                       <  1:2135066/150‑1 (MQ=255)
                                  tACACAACATCCATTAAACGCATTAGCAACACCATTACCACAAACAACACAATAACCACAGGTAACGGTGCGGGCTGACGCGTAAATGAAACACAGAAAAAAGCACGCACCTGACA‑GAGAGGTC‑TTTTTTTATAGACCAAAGTTAACGAg                                                                                                                   >  1:5953490/1‑150 (MQ=255)
                                  tACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAAagag                                                                                                                   >  1:4680892/1‑150 (MQ=255)
                                    aacaacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGt                                                                                                                 <  1:5054063/149‑1 (MQ=255)
                                     acaccaTCCATGAAACGAATTAGCACCACCATTACCACCACCATCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTa                                                                                                                <  1:3402592/150‑1 (MQ=255)
                                         caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGGCA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                                                                            >  1:2840364/1‑150 (MQ=255)
                                          aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                                                                           <  1:899273/150‑1 (MQ=255)
                                              aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAATGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCagg                                                                                                       >  1:6364796/1‑148 (MQ=255)
                                              aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTCGACCAAAGGAAACGAGGTAACAACCAt                                                                                                        >  1:4277980/1‑150 (MQ=255)
                                               tGAAACGCATTAGCACCACCATTACCACCACCATCACCCTTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                      <  2:3056435/150‑1 (MQ=255)
                                               tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCCAAGGTAACGAGGTAACAACCATGc                                                                                                      >  2:5168635/1‑150 (MQ=255)
                                               tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                      >  1:2114119/1‑150 (MQ=255)
                                               tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                      >  1:4494203/1‑150 (MQ=255)
                                               agAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                      <  2:2747660/149‑1 (MQ=255)
                                                 aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                                                                                    >  1:3504919/1‑150 (MQ=255)
                                                   aCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACACAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                                                                                  <  1:1836045/150‑1 (MQ=255)
                                                     gCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtgt                                                                                                >  2:570250/1‑150 (MQ=255)
                                                     gCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtgt                                                                                                >  2:5031951/1‑150 (MQ=255)
                                                     gCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACAAAAGGTAACGAGGTAACAACCATGCGAgtgt                                                                                                >  2:5322357/1‑150 (MQ=255)
                                                      cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                               <  2:2509478/150‑1 (MQ=255)
                                                         tAGCACCCCCATTACCACCACCATCACCATTACCCCCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            <  1:1961299/150‑1 (MQ=255)
                                                         tAGCACCACCATTACCCCCACCCTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            <  1:3140129/150‑1 (MQ=255)
                                                         tAGCACCACCATTACCACCACCCTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            <  2:4680125/150‑1 (MQ=255)
                                                         tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            <  1:4830646/150‑1 (MQ=255)
                                                              ccaccaTTACCACCACCATCACCATTACCACAGTTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGACCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAATTTcg                                                                                       <  2:2224776/150‑1 (MQ=255)
                                                              ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                       >  1:1515101/1‑150 (MQ=255)
                                                               caccaTTACCACCACCCTCCCCATTACCACAGGTAACGGGGCGGGCGGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACC‑GTGCGGGC‑TGTTGTTTTCGACCAAAGGGACCGAGGGCACAACCATGCGCGGGGTGAAGTTcgg                                                                                      >  2:4152960/1‑150 (MQ=255)
                                                               caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                      <  1:662316/150‑1 (MQ=255)
                                                                 ccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcg                                                                                    >  2:2979342/1‑150 (MQ=255)
                                                                 ccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcg                                                                                    <  2:3669887/150‑1 (MQ=255)
                                                                 ccaTTACCACCACCATCACCATTACAACTGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcg                                                                                    <  1:6122652/150‑1 (MQ=255)
                                                                        caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                             <  1:3240635/150‑1 (MQ=255)
                                                                           caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                          <  1:302816/150‑1 (MQ=255)
                                                                             ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGCGTGGTGAAGTTCGGCGGTACATCAGTg                                                                         >  1:4033983/1‑150 (MQ=255)
                                                                               atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                                                                      >  1:3670060/1‑150 (MQ=255)
                                                                                tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGGTGAAGTTCGGCGGTACATCAGTGGCaa                                                                     >  1:2677955/1‑150 (MQ=255)
                                                                                   ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg                                                                  <  2:2446449/150‑1 (MQ=255)
                                                                                    catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                 <  2:6034971/150‑1 (MQ=255)
                                                                                    catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                 <  1:744453/150‑1 (MQ=255)
                                                                                    catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                 <  1:5157953/150‑1 (MQ=255)
                                                                                     atTACCACAGGTACCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa                                                                >  1:3916629/1‑150 (MQ=255)
                                                                                      ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg                                                               >  1:364704/1‑150 (MQ=255)
                                                                                       tACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa                                                              <  2:1171501/150‑1 (MQ=255)
                                                                                       gaCCCCAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa                                                              <  2:1973639/149‑1 (MQ=255)
                                                                                       gaCCACAGGTAACGGTGCGGGCTGCCGCGTACAGGAAACACAGAAACAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTATTTTCGACCAAAGGTAACGCGGTAACAACCATGCGAGTGTTGTAGTTCGGCGGTACATCTGTGGCAAATGCAGa                                                              <  2:4586939/149‑1 (MQ=255)
                                                                                        aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa                                                             <  1:1632451/150‑1 (MQ=255)
                                                                                            caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGtt                                                         <  2:3893692/150‑1 (MQ=255)
                                                                                               gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                      >  1:1006051/1‑150 (MQ=255)
                                                                                               gTAACGGGGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGAGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                      <  2:2344008/150‑1 (MQ=255)
                                                                                                tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCt                                                     <  1:3810620/150‑1 (MQ=255)
                                                                                                 aaCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTg                                                    <  2:4366029/150‑1 (MQ=255)
                                                                                                  aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGc                                                   <  2:5646277/150‑1 (MQ=255)
                                                                                                  aCGGTGCGGGCGGACGCGTACAGGAAACACAGAACAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGc                                                   <  2:4003242/150‑1 (MQ=255)
                                                                                                          ggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                           >  2:2837317/1‑150 (MQ=255)
                                                                                                            cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGa                                         <  2:1482253/150‑1 (MQ=255)
                                                                                                            cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGa                                         <  1:623465/150‑1 (MQ=255)
                                                                                                               aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGatat                                      >  2:798810/1‑150 (MQ=255)
                                                                                                                 gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTc                                    <  2:5863280/150‑1 (MQ=255)
                                                                                                                  cgTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCt                                   <  1:5534805/150‑1 (MQ=255)
                                                                                                                   gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTATCGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                  <  1:3733088/150‑1 (MQ=255)
                                                                                                                   gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                  <  2:2607678/150‑1 (MQ=255)
                                                                                                                      cAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaa                               <  1:6396242/150‑1 (MQ=255)
                                                                                                                       aGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                              <  1:848713/150‑1 (MQ=255)
                                                                                                                        ggAAACACAGAAAAAAGCCCGCACCTGCCA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAg                             <  1:2792250/150‑1 (MQ=255)
                                                                                                                           aaCACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCCGTGGCCAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCaa                          >  2:4478925/1‑150 (MQ=255)
                                                                                                                           aaCACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCaa                          <  1:2602899/150‑1 (MQ=255)
                                                                                                                           aaCACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCaa                          <  1:5422970/150‑1 (MQ=255)
                                                                                                                            acacacAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTATTTTCGACCAAACGTAACGAGGTAACAACCATGCGAGTGATGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAt                         <  2:4904475/150‑1 (MQ=255)
                                                                                                                            acacAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAt                         <  2:1593100/150‑1 (MQ=255)
                                                                                                                              acaGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                       >  2:2152305/1‑150 (MQ=255)
                                                                                                                                aaaaaaatGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTATGGAAAGAAATGCca                     >  1:2236570/3‑150 (MQ=255)
                                                                                                                                aaaaaaaaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAGGCca                     >  2:6531537/3‑150 (MQ=39)
                                                                                                                                aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                     <  1:4141550/150‑1 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                     <  2:3406068/150‑1 (MQ=255)
                                                                                                                                      aaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTATAGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTTCCGATATTATGGAAAGCAAAGCcaggcagg               >  1:2852757/1‑150 (MQ=255)
                                                                                                                                         cccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTCACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATCTGCTGGAAAGCAATGCCAggccggggc            >  2:5547857/1‑150 (MQ=255)
                                                                                                                                          ccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGAACAAAGGTAAAGAGGTAACAAACAATCGAGTTTTGAAGTTAGGCGGTAAAAAATTGGCAAAAGGAGAACGTTTTCTTCGTGTTGCCGATATTCTTGAAAGCAAAGCCAggcaggggcc           >  1:5996886/1‑149 (MQ=255)
                                                                                                                                           cGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag          <  1:4941607/150‑1 (MQ=255)
                                                                                                                                              aCCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg       <  1:1447566/150‑1 (MQ=255)
                                                                                                                                               ccTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg      <  1:5852624/150‑1 (MQ=255)
                                                                                                                                                  gACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa   <  1:4282626/150‑1 (MQ=255)
                                                                                                                                                  gACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa   >  1:5268024/1‑150 (MQ=255)
                                                                                                                                                   aCA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAc  >  2:353601/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAC  >  NC_000913/144‑440

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