Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,968,908 T→C 6.0% S153G (AGT→GGT)  cheR ← chemotaxis regulator, protein‑glutamate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,968,9080TC 6.0% 204.0 / 2.6 66S153G (AGT→GGT) cheRchemotaxis regulator, protein‑glutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base C (2/2);  total (32/34)
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.92e-01

AGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGACGTGCGTGATCCGCGAGCAGAGG  >  NC_000913/1968760‑1969055
                                                                                                                                                    |                                                                                                                                                   
aGCTCATGACGCACGCGTACCAGCCCTTCATGCGGACCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACt                                                                                                                                                     <  5:370734/149‑1 (MQ=255)
     cTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCa                                                                                                                                                 >  4:118533/1‑148 (MQ=255)
              gcgTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTc                                                                                                                                        >  1:30013/1‑148 (MQ=255)
               cgTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCa                                                                                                                                      <  1:300975/149‑1 (MQ=255)
               cgTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCa                                                                                                                                      <  2:293148/149‑1 (MQ=255)
                    cAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTc                                                                                                                                 >  5:292208/1‑149 (MQ=255)
                      gCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTc                                                                                                                                        >  7:15232/1‑140 (MQ=255)
                      gCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTc                                                                                                                                        <  8:15232/140‑1 (MQ=255)
                      gCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCg                                                                                                                               >  3:335943/1‑149 (MQ=255)
                      gCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCg                                                                                                                               >  2:160881/1‑149 (MQ=255)
                       tccTTCATGGGGCCGCGTCCGTCGCATGTAATAAGGGTGCAGTTGGTGCGGCGTCAGGTTTTTCAACTCTTCCTGGCGATAGATACCGCTTCTGGCTTTTTCCAGTACTTCTGTGTCGTTATCACTGGTAAACTGTTTCCAGAGTTCg                                                                                                                               <  2:87658/147‑1 (MQ=255)
                         cTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGgc                                                                                                                            <  2:30013/149‑1 (MQ=255)
                            cATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGgcgc                                                                                                                          <  1:160881/148‑1 (MQ=255)
                            cATGCGGCCCCGTCCCTCGCATGAAAAACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGcgcg                                                                                                                         >  8:277375/1‑149 (MQ=255)
                             aTGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCgg                                                                                                                        >  6:238117/1‑149 (MQ=255)
                                cGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGc                                                                                                                     <  6:292208/149‑1 (MQ=255)
                                     ccGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAAt                                                                                                                <  3:118533/149‑1 (MQ=255)
                                     ccGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAAt                                                                                                                <  5:156456/149‑1 (MQ=255)
                                          ccTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGcgcg                                                                                                                         >  6:74634/1‑135 (MQ=255)
                                            tCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGAGTCAGGGCGCGGTGCCCAATGTGTcag                                                                                                         >  5:20618/1‑149 (MQ=255)
                                               cATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTcagcca                                                                                                      <  5:238117/149‑1 (MQ=255)
                                                 tGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTcagccagc                                                                                                    >  7:10058/1‑149 (MQ=255)
                                                     aTACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAAtgt                                                                                                              <  7:358984/135‑1 (MQ=255)
                                                       aCCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATTCCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCTCGGTGCCCAATGTGTCAGCCAGCGTCATc                                                                                              >  2:234563/1‑149 (MQ=255)
                                                          gTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGc                                                                                            >  1:350063/1‑148 (MQ=255)
                                                           ttGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTTTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCaa                                                                                          <  2:376040/149‑1 (MQ=255)
                                                              cAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGTCTCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGc                                                                                       <  1:234563/149‑1 (MQ=255)
                                                               aGTTGCTGCGGCGTCAGGTTTTTCAACGCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGAGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCt                                                                                      <  6:20618/149‑1 (MQ=255)
                                                                   gctgcGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCGAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTa                                                                                   <  7:277375/148‑1 (MQ=255)
                                                                            tCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCg                                                                         >  1:247412/1‑149 (MQ=255)
                                                                             cAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGcgcg                                                                                                                         >  4:326373/1‑100 (MQ=255)
                                                                             cAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGcgcg                                                                                                                         <  3:326373/100‑1 (MQ=255)
                                                                             cAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGc                                                                        >  7:350990/1‑149 (MQ=255)
                                                                              aGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACCGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGcc                                                                       >  8:58174/1‑149 (MQ=255)
                                                                              aGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACCGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGcc                                                                       >  8:52224/1‑149 (MQ=255)
                                                                                   tttCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTc                                                                   <  8:10058/148‑1 (MQ=255)
                                                                                      cAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGaa                                                                <  4:335943/148‑1 (MQ=255)
                                                                                       aaCTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAgc                                                              >  6:341302/1‑149 (MQ=255)
                                                                                       aaCTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAgc                                                              >  6:56383/1‑149 (MQ=255)
                                                                                               aTGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATc                                                                                              >  1:274515/1‑109 (MQ=255)
                                                                                               aTGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATc                                                                                              <  2:274515/109‑1 (MQ=255)
                                                                                                    gatagataCCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCaa                                                                                                                                                >  7:319371/1‑54 (MQ=255)
                                                                                                    gatagataCCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCaa                                                                                                                                                <  8:319371/54‑1 (MQ=255)
                                                                                                        gataCCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATAcgc                                             <  2:247412/149‑1 (MQ=255)
                                                                                                          taCCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGa                                           >  1:44345/1‑149 (MQ=255)
                                                                                                             cGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACCGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATAc                                        <  7:52224/149‑1 (MQ=255)
                                                                                                             cGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACCGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATAc                                        <  7:58174/149‑1 (MQ=255)
                                                                                                               gagaaggCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTc                                      <  1:204549/144‑1 (MQ=255)
                                                                                                                ttCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCg                                     >  7:30965/1‑149 (MQ=255)
                                                                                                                     ggtatttCCAGCACTTCGGTGTCGATATCACTGCCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATTTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTTGAAGCCGCCGCGCTCCATACGCGATACTCGCCAg                                 >  6:209024/5‑148 (MQ=255)
                                                                                                                       tttttCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGAGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAAc                              >  3:88250/1‑149 (MQ=255)
                                                                                                                         tttCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAAc                              <  2:44345/147‑1 (MQ=255)
                                                                                                                          ttCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACgg                            <  2:350063/148‑1 (MQ=255)
                                                                                                                             cAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGAc                        <  8:350990/149‑1 (MQ=255)
                                                                                                                              aGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGaa                                                                >  3:94762/1‑108 (MQ=255)
                                                                                                                              aGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGaa                                                                <  4:94762/108‑1 (MQ=255)
                                                                                                                                 aCTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGAcgtgc                    <  5:341302/149‑1 (MQ=255)
                                                                                                                                   ttCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGAcgtgcgt                  >  5:282842/1‑149 (MQ=255)
                                                                                                                                   ttCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGAcgtgcgt                  >  2:155946/1‑149 (MQ=255)
                                                                                                                                         gtCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGACGTGCGTGATCcg            >  7:309728/1‑149 (MQ=255)
                                                                                                                                            gATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGACGTGCGTGATCCGCGa         <  5:56383/149‑1 (MQ=255)
                                                                                                                                                 cccTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGACGTGCGTGATCCGCGa         <  5:187681/142‑1 (MQ=255)
                                                                                                                                                 cACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGACGTGCGTGATCCGCGa         >  6:187681/1‑144 (MQ=255)
                                                                                                                                                  aCTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGACGTGCGTGATCCGCGAGCaga    <  8:30965/148‑1 (MQ=255)
                                                                                                                                                   cTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGACGTGCGTGATCCGCGa         <  5:262268/142‑1 (MQ=255)
                                                                                                                                                   cTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGACGTGCGTGATCCGCGa         >  6:262268/1‑142 (MQ=255)
                                                                                                                                                   cTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGACGTGCGTGATCCGCGAGCAGAgg  <  6:282842/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCGCAATGCTGTACGGCTCTTCGCCGGTCGAAGCCGCCGCGCTCCATACGCGATACTCGCCAGAACGGCGACGTGCGTGATCCGCGAGCAGAGG  >  NC_000913/1968760‑1969055

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