Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 293:1 +G 5.2% 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.2% 2060.8 / 101.2 446intergenic (+38/‑44)thrL/thrAthr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base . (195/228);  new base G (12/11);  total (207/239)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

Only 100 of 459 total aligned reads displayed.

AGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCG  >  NC_000913/144‑442
                                                                                                                                                      |                                                                                                                                                      
aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACa                                                                                                                                                         >  2:2199368/1‑150 (MQ=255)
aGGCAAAGCGCACAGACAGATAAAAAATACAGAGAACACAACATCCATGAAACGCATTAGCACCACCAAAACCACCACCATCACCATAACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTcaca                                                                                                                                                         >  2:1651365/1‑150 (MQ=255)
     tAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCg                                                                                                                                                   <  2:4413192/150‑1 (MQ=255)
      aGCGCACATACATATAAAAATTACAGAGTACACAACATCCATTAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAAGGGTGCGGGCTGACGCGTACAGGAAACACATAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                  <  1:505467/150‑1 (MQ=255)
      aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGACACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                  >  2:1815969/1‑150 (MQ=255)
        cgcAAAGACAGATAAAAATTACAGATTACACACACTACATGAAACGATTTAGCACCAGCAGTACCACCACCATAACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACTGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGc                                                                                                                                                <  1:7112780/150‑1 (MQ=255)
         gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑t                                                                                                                                              >  2:7759558/1‑150 (MQ=255)
          cacaGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCCCCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑tt                                                                                                                                             <  1:1145728/150‑1 (MQ=255)
            cagacagaTAAAAATTACAGAGGACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACCGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑tttt                                                                                                                                           <  2:4956859/150‑1 (MQ=255)
               acagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttt                                                                                                                                        >  1:3772488/1‑150 (MQ=255)
               acagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttt                                                                                                                                        <  1:5690489/150‑1 (MQ=255)
                  taTAAAAATTACAGATTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACTGGTAAGGGTGCGGGCTGACGCTTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTc                                                                                                                                     <  1:1826353/149‑1 (MQ=255)
                    tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCCCCATTACCCCCACCAACACCATTACCACAGGTACCGGTGCGGGCTGACGCGTACAGGAAACACACAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttttttctgtgt                                                                                                                                   <  1:2161777/150‑7 (MQ=37)
                    tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGa                                                                                                                                   <  2:5742779/150‑1 (MQ=255)
                    aaaaaaTGACAGATTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCCCTCGTAACGGTGCGGGCTCACGCGTACAGTAAACACCCAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGa                                                                                                                                   <  1:8618735/149‑1 (MQ=255)
                     aaaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCAACACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACATAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGAc                                                                                                                                  <  1:6204918/150‑1 (MQ=255)
                     aaaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACTGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGAc                                                                                                                                  <  1:8945556/150‑1 (MQ=255)
                         aTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCaaa                                                                                                                              <  2:6493068/150‑1 (MQ=255)
                               gagTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAATAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAAc                                                                                                                        >  1:345097/1‑150 (MQ=255)
                                  tACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAg                                                                                                                     <  1:7340564/150‑1 (MQ=255)
                                       aacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACAATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                                                                                <  1:5099193/150‑1 (MQ=255)
                                         caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                                                                              >  1:3143740/1‑150 (MQ=255)
                                           tCCATGAAACGCATTAGCACCACCATTACCACCACCATCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                                                                                            <  1:1678679/150‑1 (MQ=255)
                                            ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                           >  2:7308091/1‑150 (MQ=255)
                                              aTGAAACGCATAAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                         >  1:7531503/1‑150 (MQ=255)
                                               tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                        >  1:8577704/1‑150 (MQ=255)
                                                 aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                                                                                      <  1:8766955/150‑1 (MQ=255)
                                                 aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                                                                                      >  2:3853482/1‑150 (MQ=255)
                                                    cGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                                                                                    <  2:7136188/150‑1 (MQ=255)
                                                        ttAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                                                                               >  1:1713869/1‑150 (MQ=255)
                                                          aGCACCACCATTACTACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                                                                             >  2:1363250/1‑150 (MQ=255)
                                                          aGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                                                                             >  1:9176689/1‑150 (MQ=255)
                                                             atcaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTGTGTTCGACCACAGGTAACGAGGTAACACCCCTGCGCGGGTTGAAGTTc                                                                                          >  2:7676071/3‑150 (MQ=255)
                                                             accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                          <  1:5550461/150‑1 (MQ=255)
                                                               caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                        >  1:8707048/1‑150 (MQ=255)
                                                                accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                       <  1:3902905/150‑1 (MQ=255)
                                                                accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                       <  1:5444321/150‑1 (MQ=255)
                                                                 ccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcg                                                                                      <  2:3488468/150‑1 (MQ=255)
                                                                    ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCTTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                   <  1:7411994/150‑1 (MQ=255)
                                                                    ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                   <  2:8715671/150‑1 (MQ=255)
                                                                     tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                                                                                  <  2:3772488/150‑1 (MQ=255)
                                                                       ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                                                                                >  1:6913945/1‑150 (MQ=255)
                                                                        caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                               >  2:8073253/1‑150 (MQ=255)
                                                                        caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                               >  2:3250996/1‑150 (MQ=255)
                                                                         accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAATAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGTCCAATGGTAACGTGGTTACAACCTTTCGTGTTTTGTAGTTAGGCGGTACAtga                                                                              >  2:1597274/1‑148 (MQ=255)
                                                                           caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                            <  2:7531503/150‑1 (MQ=255)
                                                                           caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGACCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                            >  1:6416553/1‑150 (MQ=255)
                                                                            accatcaccatTACCACAGGTAACGGTGCAGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                                                                           >  2:6176613/1‑150 (MQ=255)
                                                                                tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa                                                                       >  1:1226655/1‑150 (MQ=255)
                                                                                 caccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaaa                                                                      >  2:5011634/1‑150 (MQ=255)
                                                                                  accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt                                                                     <  1:7743402/150‑1 (MQ=255)
                                                                                   ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg                                                                    <  1:2287708/150‑1 (MQ=255)
                                                                                    catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                   >  1:5783762/1‑150 (MQ=255)
                                                                                    catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                   >  1:1069603/1‑150 (MQ=255)
                                                                                     atTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa                                                                  <  1:4504181/150‑1 (MQ=255)
                                                                                      ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg                                                                 >  2:7997644/1‑150 (MQ=255)
                                                                                             aGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGttt                                                          >  2:8980490/1‑150 (MQ=255)
                                                                                               gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                        <  1:5291604/150‑1 (MQ=255)
                                                                                               gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCCAGTGGTGAAGGTCGGCGGTACCTCGGTGGCAAATGCAGAACGGTTTc                                                        >  2:8625633/1‑150 (MQ=255)
                                                                                                tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCt                                                       <  2:524434/150‑1 (MQ=255)
                                                                                                  aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCCAAGGTAACGAGGTAACAACCATGCGAGTGGTGAAGTTCGGCGGGACATCAGTGGCAAATGCAGAACGTTTTCTGc                                                     >  1:2448682/1‑150 (MQ=255)
                                                                                                    ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                   <  1:8156460/150‑1 (MQ=255)
                                                                                                    ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                   <  2:2565634/150‑1 (MQ=255)
                                                                                                    ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                   <  2:2975066/150‑1 (MQ=255)
                                                                                                     gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgtg                                                  >  2:5476153/1‑150 (MQ=255)
                                                                                                        cGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTg                                               <  2:2347097/150‑1 (MQ=255)
                                                                                                          ggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                             <  1:830634/150‑1 (MQ=255)
                                                                                                                 gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTc                                      >  1:4944049/1‑150 (MQ=255)
                                                                                                                   gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                    <  1:9318743/150‑1 (MQ=255)
                                                                                                                    tACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGGGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAAGTTTTCTGCGTGTTGCCGATATTCTgg                                   >  1:6018143/1‑150 (MQ=255)
                                                                                                                    tACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                                   >  1:402385/1‑150 (MQ=255)
                                                                                                                     aCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                  <  2:3143785/150‑1 (MQ=255)
                                                                                                                       aGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                                <  1:9096746/150‑1 (MQ=255)
                                                                                                                         gAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACCTCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                              >  2:4640355/1‑150 (MQ=255)
                                                                                                                         gAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                              >  2:2814173/1‑150 (MQ=255)
                                                                                                                         gAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                              >  2:4338622/1‑150 (MQ=255)
                                                                                                                         gAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                              >  2:1533051/1‑150 (MQ=255)
                                                                                                                         gAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTGGAAGTGCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                              >  2:2723566/1‑150 (MQ=255)
                                                                                                                         gAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAAAGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                              >  2:1219868/1‑150 (MQ=255)
                                                                                                                            acacAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAt                           <  1:6697739/150‑1 (MQ=255)
                                                                                                                              acaGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTCTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                         <  1:7976596/150‑1 (MQ=255)
                                                                                                                                aGAACAAAGCCCGCCCCTGACA‑GTGCGGGC‑TTTGTTTTTCGACCACAGGTACCGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCCAATGCAGACCGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                       <  1:7676071/150‑1 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCGAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                       <  1:2640443/150‑1 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                       <  2:6109046/150‑1 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                       <  2:9141337/150‑1 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAATGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                       >  2:6882585/1‑150 (MQ=255)
                                                                                                                                aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ATTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                       <  2:8253880/150‑1 (MQ=255)
                                                                                                                                  aaaaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                     <  1:4260903/150‑1 (MQ=255)
                                                                                                                                   aaaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGTAAAAAACAAGCGAGTGTTGAATTTAGGCGGTACATCAGTGGCAAATGCAGAACGATTTCTGCGTGTTTCCGATATTCTGGAAAGCAATGCcaggc                    >  1:4828758/1‑150 (MQ=255)
                                                                                                                                      aaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg                 >  2:4113144/1‑150 (MQ=255)
                                                                                                                                      aaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGTAAAAACCAAGCGAGTGTTGAAGTTCGGCGGTACATCAGTTGCAAATTCAGAACGTTTTGTGCGTGTTGCCGATATGATTGAAAGCAAAGCcaggcagg                 >  1:7853873/1‑150 (MQ=255)
                                                                                                                                       aGCCCGCA‑CTGTCAGGTGCGGGC‑‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg               >  2:8491850/1‑150 (MQ=255)
                                                                                                                                          ccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGGAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGGACATCAGTGGCCAATTCAGAACGTTTTATGCGTGTTGCAGATATTATGGAAAGCAATGCCAggcaggggca             >  1:8408683/1‑150 (MQ=255)
                                                                                                                                           cGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag            <  2:5166052/150‑1 (MQ=255)
                                                                                                                                           cGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag            <  1:5889365/150‑1 (MQ=255)
                                                                                                                                             cACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGt          <  1:9055229/150‑1 (MQ=255)
                                                                                                                                              aCCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg         <  1:3215531/150‑1 (MQ=255)
                                                                                                                                                cTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc       <  2:4778282/150‑1 (MQ=255)
                                                                                                                                                cTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc       >  1:7192357/1‑150 (MQ=255)
                                                                                                                                                 tGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGcc      >  1:1435836/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
AGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑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: ‑