Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,930,410 G→C 6.4% T188R (ACA→AGA) ‡ cheZ ← protein phosphatase CheZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,930,4100GC 6.4% 247.2 / 5.9 78T188R (ACA→AGA) ‡cheZprotein phosphatase CheZ
Reads supporting (aligned to +/- strand):  ref base G (35/38);  new base C (3/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

CGGACCTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATCACGAAATTCCG  >  NC_002947/4930265‑4930557
                                                                                                                                                 |                                                                                                                                                   
cGGACCTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTgtgg                                                                                                                                                  <  4:130529/149‑1 (MQ=255)
         aCCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACgcgc                                                                                                                                          <  7:31783/148‑1 (MQ=255)
          ccccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                        >  5:167506/1‑149 (MQ=255)
          ccccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                        >  7:15458/1‑149 (MQ=255)
          ccccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                        >  3:274293/1‑149 (MQ=255)
          ccccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                        <  1:96853/149‑1 (MQ=255)
           cccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                        >  2:6627/1‑148 (MQ=255)
           cccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                        >  5:124835/1‑148 (MQ=255)
           cccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                        >  7:19467/1‑148 (MQ=255)
           cccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                        >  1:53303/1‑148 (MQ=255)
                    gATGTTTTTCTCGATCTTTTTCTGCACGGAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTgg                                                                                                                              >  1:214546/1‑149 (MQ=255)
                      tGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGcc                                                                                                                            <  4:272916/149‑1 (MQ=255)
                        tttttCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCg                                                                                                                           <  6:167506/148‑1 (MQ=255)
                             ctcGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAg                                                                                                                      <  2:95790/148‑1 (MQ=255)
                                gATCTTTTTCTGCACGCAATTGATGGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTcc                                                                                                                  >  8:173806/1‑149 (MQ=255)
                                gATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGa                                                                                                                                      >  2:193839/1‑129 (MQ=255)
                                gATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGa                                                                                                                                      <  1:193839/129‑1 (MQ=255)
                                gATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTc                                                                                                                   <  5:278321/148‑1 (MQ=255)
                                  tCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGTGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCt                                                                                                                 >  7:94892/1‑148 (MQ=255)
                                      tttCTGCACGCAATTGATCGTTATCATCTTTTATACCGGCAAAGCTGTCGAACTGGCTGGACATCAGCACCAGATTGAGCAGGTTGCTTTCAACCTCGGTCACGAGTGTGGTGACGCGCTTGATCAACTGGCCGATCAGGTGCTGATAg                                                                                                            >  4:226275/1‑149 (MQ=255)
                                      tttCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGAt                                                                                                              <  6:124835/147‑1 (MQ=255)
                                            cACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTgg                                                                                                      <  5:35530/149‑1 (MQ=255)
                                               gCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGc                                                                                                    >  7:345203/1‑148 (MQ=255)
                                                cAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCATGTCCTGATAGTCCTGGGCCa                                                                                                  <  4:234434/149‑1 (MQ=255)
                                                 aaTTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAg                                                                                                 >  4:368354/1‑149 (MQ=255)
                                                      aTCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGt                                                                                                           <  7:349285/134‑1 (MQ=255)
                                                      aTCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGt                                                                                                           >  8:349285/1‑134 (MQ=255)
                                                       tCGTGGTCATGTTTTATACCTGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTTGCCGGTCAGGTCCTGATATTCCTGGGCCAGCAAAAt                                                                                           <  3:226275/149‑1 (MQ=255)
                                                        cGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACTAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATg                                                                                          >  8:100976/1‑149 (MQ=255)
                                                          tGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTc                                                                                        >  2:297360/1‑149 (MQ=255)
                                                          tGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTc                                                                                        <  2:214546/149‑1 (MQ=255)
                                                             tCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGtt                                                                                     <  8:924/149‑1 (MQ=255)
                                                              cATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGtt                                                                                     >  4:85904/1‑148 (MQ=255)
                                                               aTGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGTTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTg                                                                                    <  7:173806/148‑1 (MQ=255)
                                                               aTGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTg                                                                                    <  3:44724/148‑1 (MQ=255)
                                                               aTGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGa                                                                                   <  8:15458/149‑1 (MQ=255)
                                                               aTGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCACTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGtt                                                                                     >  2:141936/1‑147 (MQ=255)
                                                                 gTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAg                                                                                  <  8:345203/148‑1 (MQ=255)
                                                                  ttttATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGt                                                                                >  6:306501/1‑149 (MQ=255)
                                                                  ttttATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGt                                                                                >  6:166061/1‑149 (MQ=255)
                                                                     tataCCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGcc                                                                             >  3:155152/1‑149 (MQ=255)
                                                                      ataCCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGcc                                                                             <  2:53303/148‑1 (MQ=255)
                                                                           ggCAAAGTGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGaaa                                                                       <  8:94892/149‑1 (MQ=255)
                                                                                 gCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGc                                                                 <  1:6627/149‑1 (MQ=255)
                                                                                        aCCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTAcc                                                          <  1:297360/149‑1 (MQ=255)
                                                                                          ctggctggCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATATTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCt                                                        >  1:344684/1‑149 (MQ=255)
                                                                                            ggctggcCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCg                                                      >  1:172561/1‑149 (MQ=255)
                                                                                             gctggcCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCGCCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCg                                                      <  4:274293/148‑1 (MQ=255)
                                                                                               tggcCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGgcg                                                   >  5:250915/1‑149 (MQ=255)
                                                                                                 gcCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCt                                                 <  3:85904/149‑1 (MQ=255)
                                                                                                   cATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCt                                                 >  4:41456/1‑147 (MQ=255)
                                                                                                   cATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCt                                                 <  3:41456/147‑1 (MQ=255)
                                                                                                    aTCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCAGCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAATCCTTGCGGTTACCCTCGGCGCtgtg                                              <  1:141936/149‑1 (MQ=255)
                                                                                                        gCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATg                                                                                          >  1:283781/1‑101 (MQ=255)
                                                                                                        gCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATg                                                                                          <  2:283781/101‑1 (MQ=255)
                                                                                                         cACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCa                                         <  4:155152/149‑1 (MQ=255)
                                                                                                            cAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCtgt                                               >  4:278672/1‑140 (MQ=255)
                                                                                                            cAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCtgt                                               <  3:278672/140‑1 (MQ=255)
                                                                                                                   aGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTc                               <  8:19467/149‑1 (MQ=255)
                                                                                                                    gCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTc                               >  6:18966/1‑148 (MQ=255)
                                                                                                                        gTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAgg                                                                                                                     <  5:229700/58‑1 (MQ=255)
                                                                                                                        gTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAgg                                                                                                                     >  6:229700/1‑58 (MQ=255)
                                                                                                                           gCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGc                        <  6:250915/148‑1 (MQ=255)
                                                                                                                            cTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCt                       <  5:166061/148‑1 (MQ=255)
                                                                                                                               tCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGa                    <  2:172561/148‑1 (MQ=255)
                                                                                                                                    gcgcgTCACCACACTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACt                                  >  7:248605/5‑129 (MQ=255)
                                                                                                                                    gcgcgTCACCACACTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACt                                  <  8:248605/125‑1 (MQ=255)
                                                                                                                                      gcgTCACCACACTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCGTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACACTCGGCGCTGTGCGTCAGTAAGCTGTCGACTCGGTTGCGCAGATCa            >  8:342907/3‑149 (MQ=255)
                                                                                                                                      gcgTCACCACACTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATc             >  2:94985/3‑148 (MQ=255)
                                                                                                                                      gcgTCACCACACTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATc             <  1:94985/146‑1 (MQ=255)
                                                                                                                                       ggTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATCAc           >  7:297753/1‑149 (MQ=255)
                                                                                                                                        gTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGc                        >  1:71195/1‑135 (MQ=255)
                                                                                                                                        gTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGc                        <  2:71195/135‑1 (MQ=255)
                                                                                                                                        gTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATCAc           >  8:55796/1‑148 (MQ=255)
                                                                                                                                        gTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATCACg          <  7:100976/149‑1 (MQ=255)
                                                                                                                                            ccAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGgcg                                                   <  2:260240/104‑1 (MQ=255)
                                                                                                                                            ccAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGgcg                                                   >  1:260240/1‑104 (MQ=255)
                                                                                                                                             cAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAg                                                                                  <  2:327621/72‑1 (MQ=255)
                                                                                                                                             cAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAg                                                                                  >  1:327621/1‑72 (MQ=255)
                                                                                                                                                tgtgGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATCACGAAATTCCg  <  5:18966/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGGACCTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATCACGAAATTCCG  >  NC_002947/4930265‑4930557

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