Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,930,408 G→C 7.3% L189V (CTG→GTG)  cheZ ← protein phosphatase CheZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,930,4080GC 7.3% 248.1 / 6.7 79L189V (CTG→GTG) cheZprotein phosphatase CheZ
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base C (4/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CGGACCTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATCACG  >  NC_002947/4930265‑4930549
                                                                                                                                               |                                                                                                                                             
cGGACCTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAgt                                                                                                                                              <  6:338893/145‑1 (MQ=255)
cGGACCTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTgtgg                                                                                                                                          <  4:130529/149‑1 (MQ=255)
         aCCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACgcgc                                                                                                                                  <  7:31783/148‑1 (MQ=255)
          ccccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                >  5:167506/1‑149 (MQ=255)
          ccccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                <  1:96853/149‑1 (MQ=255)
          ccccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                >  3:274293/1‑149 (MQ=255)
          ccccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                >  7:15458/1‑149 (MQ=255)
           cccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                >  2:6627/1‑148 (MQ=255)
           cccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                >  7:19467/1‑148 (MQ=255)
           cccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                >  1:53303/1‑148 (MQ=255)
           cccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCtt                                                                                                                                >  5:124835/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                                                                                                   >  8:349285/1‑134 (MQ=255)
                                                      aTCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGt                                                                                                   <  7:349285/134‑1 (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)
                                                               aTGTTTTATACCGGCAAAGCGTTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAgt                                                                                                                                              <  4:235894/82‑1 (MQ=255)
                                                               aTGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAgt                                                                                                                                              >  3:235894/1‑82 (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     <  1:94985/146‑1 (MQ=255)
                                                                                                                                      gcgTCACCACACTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATc     >  2:94985/3‑148 (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)
                                                                                                                                               |                                                                                                                                             
CGGACCTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATCACG  >  NC_002947/4930265‑4930549

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