Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,213,366 C→T 6.0% A427T (GCC→ACC)  yehU ← inner membrane putative sensory kinase in two‑component system with YehT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,213,3660CT 6.0% 209.2 / 2.4 66A427T (GCC→ACC) yehUinner membrane putative sensory kinase in two‑component system with YehT
Reads supporting (aligned to +/- strand):  ref base C (28/34);  new base T (2/2);  total (30/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 = 4.57e-01

TGGCGTTTTCCACTATCGGTTGCAGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGGGGATTCACCTGG  >  NC_000913/2213219‑2213514
                                                                                                                                                   |                                                                                                                                                    
tGGCGTTTTCCACTATCGGTTGCAGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGc                                                                                                                                                      >  8:307073/1‑148 (MQ=255)
   cGTTTTCCACTATCGGTTGCAGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCgaga                                                                                                                                                  <  3:316979/149‑1 (MQ=255)
   cGTTTTCCACTATCGGTTGCAGGGTAAACGCGGGCAATTGCTGATGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGAGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGTCgaga                                                                                                                                                  <  4:155664/149‑1 (MQ=255)
     ttttCCACTATCGGTTGCAGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGagag                                                                                                                                                 >  8:33969/1‑148 (MQ=255)
          cACTATCGGTTGCAGGGTAAACGCGGGCAATTGCTGCTGGGATAAATCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCTCCTTTTCAATTTGCAGATAAGCATTCACATGATCAATTTAGTCGCCGAGAGTAACa                                                                                                                                            >  7:333977/1‑148 (MQ=255)
           aCTATCGGTTGCAGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACaaa                                                                                                                                          >  7:111168/1‑149 (MQ=255)
             tATCGGTTGCAGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACt                                                                                                                                        >  8:256734/1‑149 (MQ=255)
             tATCGGTTGCAGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCACCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACt                                                                                                                                        >  3:263103/1‑149 (MQ=255)
                  gTTGCAGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCGTTTTCAATTTGCAGATAAGCGTTCACATGTTCAATTTCGTCGGCGAGAGGAACAGAGTCCg                                                                                                                                     >  4:228735/1‑147 (MQ=255)
                  gTTGCAGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCg                                                                                                                                     >  7:335333/1‑147 (MQ=255)
                       aGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCg                                                                                                                              <  8:46416/149‑1 (MQ=255)
                       aGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCg                                                                                                                              <  2:126793/149‑1 (MQ=255)
                       aGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGATGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGAGAAGCATTCACATGTTCATTTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCg                                                                                                                              >  2:305676/1‑149 (MQ=255)
                          gTAAACGGGGGCAATTGCTGCTTGGATAATTCTTGCGGAATAGCATTGTTGACCTGCAACAGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCtt                                                                                                                           <  8:333977/149‑1 (MQ=255)
                          gTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCtt                                                                                                                           <  7:256734/149‑1 (MQ=255)
                              aCGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCACCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAg                                                                                                                       <  4:263103/149‑1 (MQ=255)
                                   ggCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATCTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGttttt                                                                                                                  <  1:326205/149‑1 (MQ=255)
                                     cAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGc                                                                                                                <  2:85630/149‑1 (MQ=255)
                                      aaTTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGCCCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCg                                                                                                               >  4:111742/1‑149 (MQ=255)
                                      aaTTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCg                                                                                                               <  4:190938/149‑1 (MQ=255)
                                        ttGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGTGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCgg                                                                                                              >  2:142805/1‑148 (MQ=255)
                                                 ggATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTg                                                                                                     <  7:307073/148‑1 (MQ=255)
                                                  gATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTgg                                                                                                    >  6:257379/1‑148 (MQ=255)
                                                  gATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGa                                                                                                   <  4:295874/149‑1 (MQ=255)
                                                  gATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGa                                                                                                   <  8:73946/149‑1 (MQ=255)
                                                      aTTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGa                                                                                               <  5:257379/149‑1 (MQ=255)
                                                       ttCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCATCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGaaa                                                                                                 <  8:80378/146‑1 (MQ=255)
                                                       ttCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCATCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGaaa                                                                                                 >  7:80378/1‑146 (MQ=255)
                                                             cGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCAGTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCa                                                                                        <  1:250894/149‑1 (MQ=255)
                                                                  tAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAgc                                                                                   <  8:335333/149‑1 (MQ=255)
                                                                  tAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCAGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAgc                                                                                   >  7:65577/1‑149 (MQ=255)
                                                                        tGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAgctggctg                                                                             >  3:88953/1‑149 (MQ=255)
                                                                        tGTAGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAgctggctg                                                                             >  7:256765/1‑149 (MQ=255)
                                                                          ttGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGctggctgg                                                                            <  5:65811/148‑1 (MQ=255)
                                                                          ttGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCtggctggc                                                                           <  8:261903/149‑1 (MQ=255)
                                                                             aCCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCt                                                                         <  8:113378/148‑1 (MQ=255)
                                                                               cTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTATAAGTTTTTGCGGAAAAAAGTGGAAAGCTACTGCACCAGCTGGCTGGccactt                                                                      >  4:299859/1‑145 (MQ=255)
                                                                                  cAAACGCTACTGGAATCGCGCCTTTTAAATTTGCAGATAAGCATTCAAATGTTCTATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCg                                                                    <  5:93164/148‑1 (MQ=255)
                                                                                   aaCCGCGACTGGAAGCGCTCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCt                                                                  <  7:236091/149‑1 (MQ=255)
                                                                                   aaCCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGTGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTc                                                                     >  2:75605/1‑146 (MQ=255)
                                                                                   aaCCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGTGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTc                                                                     <  1:75605/146‑1 (MQ=255)
                                                                                   aaCCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCt                                                                  <  8:337505/149‑1 (MQ=255)
                                                                                    aCCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGTGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTg                                                                 <  1:142805/149‑1 (MQ=255)
                                                                                       gcgACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTcg                                                              >  3:84151/1‑149 (MQ=255)
                                                                                              gAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGa                                                       >  7:292627/1‑149 (MQ=255)
                                                                                                 gcgcgcCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCa                                                    >  4:240174/1‑149 (MQ=255)
                                                                                                       cTTTTCAATTTGCAGATTAGCATTCAGATGTTCAATTTCGTAGTCGCGAGTAAGAAACTCCGAAGGGCGCGTTAAGTTTTTGCGGTAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCtt                                              <  1:129787/149‑1 (MQ=255)
                                                                                                       cTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCtt                                              <  3:240174/149‑1 (MQ=255)
                                                                                                             aaTTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATgg                                        <  7:33969/149‑1 (MQ=255)
                                                                                                             aaTTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATg                                         <  8:141042/148‑1 (MQ=255)
                                                                                                                    aGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAg                                 >  6:286159/1‑149 (MQ=255)
                                                                                                                       tAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAgcgc                              <  1:319698/149‑1 (MQ=255)
                                                                                                                       tAAGCATTCACATGTTCAATTTCGTCGGCGAAAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAgcgc                              <  8:80379/149‑1 (MQ=255)
                                                                                                                          gCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTAGCGGAAAAAAGTCGAAAGATACTGCACCAGCTGGCTAGACAGTACTCTCTCTCGGCGGAACACCGCTTTAATTGTGTTAAGCTCAt                            >  1:337101/1‑148 (MQ=255)
                                                                                                                           aaTTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCTCGGCGGATAACCGCTTTAATGGTGTTAAGCGAAtt                           >  4:109704/2‑148 (MQ=255)
                                                                                                                                acaTGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAGCGCATTAAACa                      >  4:93029/1‑148 (MQ=255)
                                                                                                                                acaTGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAGCGCATTAAACa                      >  1:176371/1‑148 (MQ=255)
                                                                                                                                  aTGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGgtgt                                     <  5:69307/131‑1 (MQ=255)
                                                                                                                                  aTGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGgtgt                                     >  6:69307/1‑131 (MQ=255)
                                                                                                                                   tGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAGCGCATTAAACaaaa                   >  7:82059/1‑148 (MQ=255)
                                                                                                                                     ttCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGgt                                       <  4:261682/126‑1 (MQ=255)
                                                                                                                                     ttCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGgt                                       >  3:261682/1‑126 (MQ=255)
                                                                                                                                        aaTTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTaa                                           >  4:222625/1‑119 (MQ=255)
                                                                                                                                        aaTTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTaa                                           <  3:222625/119‑1 (MQ=255)
                                                                                                                                        aaTTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATggg             <  3:93029/149‑1 (MQ=255)
                                                                                                                                        aaTTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATgg              <  8:292627/148‑1 (MQ=255)
                                                                                                                                            acgtcgGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGAACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAAAGGTGTTAAGCGCATTAAACAAAAAATGGGGAt          >  7:146773/2‑148 (MQ=255)
                                                                                                                                                   cGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGGGGATTCACCTgg  <  3:228735/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGCGTTTTCCACTATCGGTTGCAGGGTAAACGCGGGCAATTGCTGCTGGGATAATTCTTGCGGAATAGCAATGTTGACCTGCAACCGCGACTGGAAGCGCGCCTTTTCAATTTGCAGATAAGCATTCACATGTTCAATTTCGTCGGCGAGAGTAACAAACTCCGAAGGCCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGGGGATTCACCTGG  >  NC_000913/2213219‑2213514

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