Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,591,287 T→C 5.3% K297R (AAA→AGA)  hipA ← inactivating GltX kinase facilitating persister formation; toxin of HipAB TA pair; autokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,591,2870TC 5.3% 230.6 / 3.3 75K297R (AAA→AGA) hipAinactivating GltX kinase facilitating persister formation; toxin of HipAB TA pair; autokinase
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base C (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

GTATTCCCGTACCGCCAAGGACCGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCCTCCTGTGGCAAGCGA  >  NC_000913/1591141‑1591433
                                                                                                                                                  |                                                                                                                                                  
gtATTCCCGTACCGCCAAGGCCCGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAAttt                                                                                                                                                   >  1:308721/1‑148 (MQ=255)
       cGTACCGCCAAGGACCGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATataa                                                                                                                                            >  2:240333/1‑148 (MQ=255)
       cGTACCGCCAAGGACCGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATaaaa                                                                                                                                            >  2:16819/1‑148 (MQ=255)
             gCCAAGGACCGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATATAATCATAg                                                                                                                                      <  1:240333/148‑1 (MQ=255)
             gCCAAGGACCGGAAATGCTGAAATGATGTCATAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAg                                                                                                                                      >  5:70338/1‑148 (MQ=255)
                  ggACCGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATc                                                                                                                                 >  4:325432/1‑148 (MQ=255)
                  ggACCGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGCCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATc                                                                                                                                 >  1:366107/1‑148 (MQ=255)
                    aCCGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGTAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCtt                                                                                                                               <  4:308480/148‑1 (MQ=255)
                    aCCGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCtt                                                                                                                               <  8:208922/148‑1 (MQ=255)
                     ccGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCttt                                                                                                                              >  1:84583/1‑148 (MQ=255)
                      cGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTc                                                                                                                             <  8:348241/148‑1 (MQ=255)
                        gAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAg                                                                                                                           >  8:343221/1‑148 (MQ=255)
                        gAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGCCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAg                                                                                                                           <  2:366107/148‑1 (MQ=255)
                         aaaTGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAgc                                                                                                                          >  8:354055/1‑148 (MQ=255)
                            tGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGcc                                                                                                                       >  1:1808/1‑148 (MQ=255)
                                gAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATCTCATAAAATCATAGCGATCTTTCAGCGCCTCGc                                                                                                                   >  3:55719/1‑148 (MQ=255)
                                     gatgTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAACTCATAGCGATCTTTc                                                                                                                             >  5:363973/1‑133 (MQ=255)
                                     gatgTCGTAAAATGGCGTGAGGCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTc                                                                                                                             <  6:363973/133‑1 (MQ=255)
                                      atgTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGAc                                                                                                              <  7:357705/147‑1 (MQ=255)
                                            tAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATc                                                                                                        >  5:136816/1‑147 (MQ=255)
                                              aaaTGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGATTTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCaa                                                                                                      >  5:87248/1‑147 (MQ=255)
                                              aaaTGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCaa                                                                                                      <  6:100666/147‑1 (MQ=255)
                                               aaTGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCAATCaaa                                                                                                     >  6:291354/1‑147 (MQ=255)
                                                      tGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCAt                                                                                             <  7:343221/148‑1 (MQ=255)
                                                      tGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCAt                                                                                             >  4:342771/1‑148 (MQ=255)
                                                             aTAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCcg                                                                                       <  5:39732/147‑1 (MQ=255)
                                                             aTAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCcg                                                                                       >  6:39732/1‑147 (MQ=255)
                                                             aTAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTAAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCcgc                                                                                      >  6:159279/1‑148 (MQ=255)
                                                                aCTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCcgcgcg                                                                                   <  2:84583/148‑1 (MQ=255)
                                                                aCTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCcgcgcg                                                                                   <  3:295447/148‑1 (MQ=255)
                                                                    ccgccAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATg                                                                                <  1:16819/147‑1 (MQ=255)
                                                                     cgccAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGcc                                                                              <  3:325432/148‑1 (MQ=255)
                                                                     cgccAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGcc                                                                              <  1:100837/148‑1 (MQ=255)
                                                                      gccAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCt                                                                             <  6:70338/148‑1 (MQ=255)
                                                                        cAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAg                                                                                                 >  5:177054/1‑126 (MQ=255)
                                                                        cAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAg                                                                                                 <  6:177054/126‑1 (MQ=255)
                                                                           ccTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGcc                                                                        >  7:113640/1‑148 (MQ=255)
                                                                            cTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCt                                                                       >  6:42223/1‑148 (MQ=255)
                                                                              gAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATCTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTc                                                                      >  3:202286/1‑147 (MQ=255)
                                                                               aataaataCGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGAcccc                                                                                                           <  7:35873/109‑1 (MQ=255)
                                                                               aataaataCGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGAcccc                                                                                                           >  8:35873/1‑109 (MQ=255)
                                                                               aataaataCGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCa                                                                    >  6:259373/1‑148 (MQ=255)
                                                                                 taaataCGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATCTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATc                                                                  <  4:55719/148‑1 (MQ=255)
                                                                                  aaataAGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGATGGACACCATAAAAAAAGCCATGATCCGAGCGATCCCTGGGCCTCCATCt                                                                 >  2:235978/1‑148 (MQ=255)
                                                                                     taCGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGAt                                                              >  4:17919/1‑148 (MQ=255)
                                                                                      aCGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAAttt                                                                                                                                                   <  2:174372/62‑1 (MQ=255)
                                                                                      aCGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAAttt                                                                                                                                                   >  1:174372/1‑62 (MQ=255)
                                                                                             aGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATAttt                                                      >  5:246532/1‑148 (MQ=255)
                                                                                              gTTTTTTGCATGACCGTCCGGTGCGCCAATCAACCACTGGAAGACCTGGAATTTAATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTc                                                     <  7:7174/148‑1 (MQ=255)
                                                                                                 ttttGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCAcc                                                  >  3:110755/1‑148 (MQ=255)
                                                                                                   ttGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATCTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCg                                                 <  4:202286/147‑1 (MQ=255)
                                                                                                      cATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGa                                             <  5:291354/148‑1 (MQ=255)
                                                                                                         gACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAgg                                          >  2:205688/1‑148 (MQ=255)
                                                                                                          aCCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGt                                         <  3:17919/148‑1 (MQ=255)
                                                                                                               ccgtTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCcgcgcg                                                                                   >  5:351580/1‑101 (MQ=255)
                                                                                                               ccgtTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCcgcgcg                                                                                   <  6:351580/101‑1 (MQ=255)
                                                                                                                      gcCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTc                             <  5:42223/148‑1 (MQ=255)
                                                                                                                        cAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTTGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTg                           <  7:354055/148‑1 (MQ=255)
                                                                                                                        cAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTg                           <  1:235978/148‑1 (MQ=255)
                                                                                                                        cAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTg                           <  4:110755/148‑1 (MQ=255)
                                                                                                                         aaTCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGa                          <  2:308721/148‑1 (MQ=255)
                                                                                                                           tCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCcgcg                                                                                     >  8:380843/1‑87 (MQ=255)
                                                                                                                           tCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCcgcg                                                                                     <  7:380843/87‑1 (MQ=255)
                                                                                                                           tCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTTGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGaca                        <  5:259373/148‑1 (MQ=255)
                                                                                                                             aaCCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGAcaca                      <  1:205688/148‑1 (MQ=255)
                                                                                                                              aCCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGa                                  <  3:192233/135‑1 (MQ=255)
                                                                                                                              aCCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGa                                  >  4:192233/1‑135 (MQ=255)
                                                                                                                                cACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACatat                   <  6:246532/148‑1 (MQ=255)
                                                                                                                                  cTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATAtcc                 >  3:232774/1‑148 (MQ=255)
                                                                                                                                    ggAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATAtcctc               >  7:230344/1‑148 (MQ=255)
                                                                                                                                            tGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGTTGAAGGTAAACCGAATGTCTGACACATATCCTCCTGTGGc        >  5:254144/1‑147 (MQ=255)
                                                                                                                                            tGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCCTCCTGTGGc        >  7:390895/1‑147 (MQ=255)
                                                                                                                                            tGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCCTCCTGTGGc        >  8:128785/1‑147 (MQ=255)
                                                                                                                                            tGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCCTCCTGTGGc        >  1:276564/1‑147 (MQ=255)
                                                                                                                                             ggAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCa                                                           <  5:340813/95‑1 (MQ=255)
                                                                                                                                             ggAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCa                                                           >  6:340813/1‑95 (MQ=255)
                                                                                                                                             ggAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGa                                  >  2:295885/1‑120 (MQ=255)
                                                                                                                                             ggAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGa                                  <  1:295885/120‑1 (MQ=255)
                                                                                                                                                 tttCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCCTCCTGTGGCAAGCGa  <  6:136816/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTATTCCCGTACCGCCAAGGACCGGAAATGCTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCCTCCTGTGGCAAGCGA  >  NC_000913/1591141‑1591433

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑