Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,592,022 A→G 6.7% I52T (ATC→ACC)  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,592,0220AG 6.7% 179.5 / 3.5 60I52T (ATC→ACC) hipAinactivating GltX kinase facilitating persister formation; toxin of HipAB TA pair; autokinase
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base G (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GGCACCAACGCTGTCTCGCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTCGCCTACCCGCTGGTTGTTCATCCAAGTGACAAGT  >  NC_000913/1591877‑1592168
                                                                                                                                                 |                                                                                                                                                  
ggCACCAACGCTGTCTCGCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGatat                                                                                                                                                 <  4:9152/149‑1 (MQ=255)
    ccAACGCTGTCTCGCCCTATTTCTGACAATAAATCAAACGGTTGTCAGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTccc                                                                                                                                             <  8:375328/149‑1 (MQ=255)
      aaCGCTGTCTCGCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTcccc                                                                                                                                            >  2:109804/1‑148 (MQ=255)
               tcGCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGGTATTCCCCCTCTGCAATg                                                                                                                                  >  4:78394/1‑149 (MQ=255)
               tcGCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATg                                                                                                                                  >  4:77785/1‑149 (MQ=255)
               tcGCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATg                                                                                                                                  >  2:160068/1‑149 (MQ=255)
               tcGCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGAAACAGAGGAGGAATTCCCCCTCTGCAATg                                                                                                                                  >  4:329581/1‑149 (MQ=255)
                cGCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGTGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGATGTGATATTCCCCCTCTGCAATgg                                                                                                                                 >  5:73816/1‑149 (MQ=255)
                 gCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGATGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATgg                                                                                                                                 >  4:271247/1‑148 (MQ=255)
                 gCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATgg                                                                                                                                 <  7:87271/148‑1 (MQ=255)
                 gCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTACCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATgg                                                                                                                                 >  3:45339/1‑148 (MQ=255)
                  cccTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCa                                                                                                                               <  1:109804/149‑1 (MQ=255)
                   ccTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAg                                                                                                                              >  3:181291/1‑149 (MQ=255)
                    cTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCGATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGc                                                                                                                             <  2:317994/149‑1 (MQ=255)
                         tCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACTTTTAACGATCAGGTCACGTACAATAGGGCTATCGGTTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGaaa                                                                                                                         <  1:280/148‑1 (MQ=255)
                           tGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATgg                                                                                                                                 <  3:174531/138‑1 (MQ=255)
                           tGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATgg                                                                                                                                 >  4:174531/1‑138 (MQ=255)
                                aTAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAAc                                                                                                                 <  4:261587/149‑1 (MQ=255)
                                   aaTCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGGGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGt                                                                                                              <  5:266858/149‑1 (MQ=255)
                                             ttGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAAc                                                                                                    <  1:160068/149‑1 (MQ=255)
                                                tCTGGATTTGGGATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTTACAGGTTTTCGAAGAAGTTAAATACGGCATCAGAGGTGGTATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGc                                                                                                 <  3:329581/149‑1 (MQ=255)
                                                tCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGGTATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGc                                                                                                 <  3:78394/149‑1 (MQ=255)
                                                         ggCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAAc                                                                                         <  3:252225/148‑1 (MQ=255)
                                                          gCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCa                                                                                       <  7:334603/149‑1 (MQ=255)
                                                           cATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGaaa                                                                                                                         >  1:136685/1‑114 (MQ=255)
                                                           cATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGaaa                                                                                                                         <  2:136685/114‑1 (MQ=255)
                                                            aTGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCATGGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTACCCACt                                                                                     >  4:302630/1‑149 (MQ=255)
                                                              gATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTcc                                                                                   >  3:111446/1‑149 (MQ=255)
                                                               aTAACGTTTAACGATCCGGGCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCg                                                                                  <  4:231455/149‑1 (MQ=255)
                                                                  aCGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTcccc                                                                                                                                            <  6:84606/88‑1 (MQ=255)
                                                                  aCGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTcccc                                                                                                                                            >  5:84606/1‑88 (MQ=255)
                                                                  aCGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGt                                                                                <  4:111446/148‑1 (MQ=255)
                                                                    gTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGACGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCa                                                                             <  3:77785/149‑1 (MQ=255)
                                                                     tttAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGa                                                                                                                                                    <  8:134540/77‑1 (MQ=255)
                                                                     tttAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGa                                                                                                                                                    >  7:134540/1‑77 (MQ=255)
                                                                             tCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAg                                                                    <  8:338326/149‑1 (MQ=255)
                                                                             tCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAg                                                                    <  1:41354/149‑1 (MQ=255)
                                                                                    aCGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTgc                                                              <  7:218760/148‑1 (MQ=255)
                                                                                      gTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTgcgcg                                                           >  3:346969/1‑149 (MQ=255)
                                                                                      gTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTgcgcg                                                           >  2:348612/1‑149 (MQ=255)
                                                                                      gTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGAgcgcg                                                           >  2:109582/1‑149 (MQ=255)
                                                                                        aCAATCGTGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGcc                                                         >  8:223813/1‑149 (MQ=255)
                                                                                        aCAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGcc                                                         >  5:354774/1‑149 (MQ=255)
                                                                                            tCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTg                                                     >  8:265757/1‑149 (MQ=255)
                                                                                            tCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTg                                                     >  3:147903/1‑149 (MQ=255)
                                                                                                 cTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTaa                                                <  6:105638/149‑1 (MQ=255)
                                                                                                  tATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAAc                                               >  2:191739/1‑149 (MQ=255)
                                                                                                    tCGGGTAACAGGTTATCGAAGAAGTTTAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACtt                                             >  4:220770/1‑149 (MQ=255)
                                                                                                        gTAACAGGGTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATACCGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGCGCGCGCCGTTGGCTAACTTCGtt                                         <  4:181291/149‑1 (MQ=255)
                                                                                                             aGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTc                                    <  4:332454/149‑1 (MQ=255)
                                                                                                             aGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTc                                    <  7:223813/149‑1 (MQ=255)
                                                                                                              ggTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTCg                                   >  5:114428/1‑149 (MQ=255)
                                                                                                                     gaagaaGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTACCTCGCCTAccc                            >  1:2910/1‑149 (MQ=255)
                                                                                                                      aagaagTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTCGCCTACCCg                           >  5:196095/1‑149 (MQ=255)
                                                                                                                       agaagTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTCGCCTACCCGc                          <  6:354774/149‑1 (MQ=255)
                                                                                                                        gaagTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGtt                                         >  3:110505/1‑133 (MQ=255)
                                                                                                                        gaagTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGtt                                         <  4:110505/133‑1 (MQ=255)
                                                                                                                                       aTCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTCGCCTACCCGCTGGTTGTTCATCCAAg          <  3:220770/149‑1 (MQ=255)
                                                                                                                                           gagGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTCGCCTACCCGCTGGTTGTTCATCCAAGTGa       >  2:56004/1‑148 (MQ=255)
                                                                                                                                            agGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTCGCCTACCCGCTGGTTGTTCATCCAAGTGAc      >  7:30333/1‑148 (MQ=255)
                                                                                                                                             ggTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTCGCCTACCCGCTGGTTGTTCATCCAAGTGACaa    <  4:346969/149‑1 (MQ=255)
                                                                                                                                              gTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTCGCCTACCCGCTGGTTGTTCATCCAAGTGACAAg   >  1:29102/1‑149 (MQ=255)
                                                                                                                                               tGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTCGCCTACCCGCTGGTTGTTCATCCAAGTGACAAGt  >  5:90447/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCACCAACGCTGTCTCGCCCTATTTCTGACAATAAATCAAACGGTTGTCTGGATTTGGCATGATAACGTTTAACGATCCGGTCACGTACAATCGGGCTATCGGGTAACAGGTTATCGAAGAAGTTAAATACGGCATCAGAGGTGATATTCCCCCTCTGCAATGGCAGCGAAAGTGACAACGGTCTGGCATAACGGCTTGCTAACCACTCCGGTGCATACTTAAAGGTGTGCGCGCCGTTGGCTAACTTCGTTAACTCGCCTACCCGCTGGTTGTTCATCCAAGTGACAAGT  >  NC_000913/1591877‑1592168

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