Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,067,303 A→T 6.4% intergenic (+119/‑55) phoR → / → PP_5322 two‑component system sensor histidine kinase/metal ion transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,067,3030AT 6.4% 233.1 / 4.8 78intergenic (+119/‑55)phoR/PP_5322two‑component system sensor histidine kinase/metal ion transporter
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base T (2/3);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

AGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGTGCTGCTCAACGGCTTC  >  NC_002947/6067157‑6067450
                                                                                                                                                  |                                                                                                                                                   
aGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAg                                                                                                                                                   <  7:223939/149‑1 (MQ=255)
  cacaATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCt                                                                                                                                                 <  8:171629/149‑1 (MQ=255)
    caATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTgg                                                                                                                                               <  4:63219/149‑1 (MQ=255)
    caATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTgg                                                                                                                                               <  4:113278/149‑1 (MQ=255)
      aTTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGGCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCa                                                                                                                                             <  1:226886/149‑1 (MQ=255)
          gCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCTTTTGCTGGCAca                                                                                                                                           >  2:205873/1‑146 (MQ=255)
          gCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCTTTTGCTGGCAca                                                                                                                                           <  1:205873/147‑2 (MQ=255)
           ccAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCTTTTGCTGGCAcaggg                                                                                                                                        <  8:296791/149‑5 (MQ=255)
              aCGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACttttttt                                                                                                                                     >  1:179560/1‑149 (MQ=255)
                  aaaGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTtctct                                                                                                                                 <  1:121885/149‑1 (MQ=255)
                  aaaGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTtctct                                                                                                                                 >  4:92968/1‑149 (MQ=255)
                  aaaGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTtctc                                                                                                                                  <  5:254885/148‑1 (MQ=255)
                  aaaGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTtctc                                                                                                                                  >  6:254885/1‑148 (MQ=255)
                   aaGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTtctctc                                                                                                                                >  5:320034/1‑149 (MQ=255)
                    aGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCCACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTctctct                                                                                                                               <  4:216963/149‑1 (MQ=255)
                      gCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTAt                                                                                                                             <  1:113552/149‑1 (MQ=255)
                        ttGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTAtct                                                                                                                           >  2:56718/1‑149 (MQ=255)
                        ttGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCCACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCtatctatc                                                                                                                            >  3:13093/1‑148 (MQ=255)
                             aTTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCCCTGGGCCCGCAATGAGTCCGCTGCCCTGTGGTTTTACGGCCGCAGGCGCGAGAGTGGAGCCCCTGTGCCAGGAAAGGCGGGCACTTTGTTTGTCTGTTTCTCAAAg                                                                                                                      >  4:237009/1‑149 (MQ=255)
                                 aGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACg                                                                                                                   <  7:16150/148‑1 (MQ=255)
                                  ggCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTtctc                                                                                                                                  <  8:105926/132‑1 (MQ=255)
                                   gCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGAc                                                                                                                >  5:88423/1‑149 (MQ=255)
                                     tGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCTTTTGCTGGCAc                                                                                                                                            <  6:213492/119‑1 (MQ=255)
                                     tGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCTTTTGCTGGCAc                                                                                                                                            >  5:213492/1‑119 (MQ=255)
                                        ggcgCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGt                                                                                                           <  8:31819/149‑1 (MQ=255)
                                           gcAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTaa                                                                                                         >  5:95135/1‑148 (MQ=255)
                                              gCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATGTCAAAGACGGACCCCGTAAGAAc                                                                                                     >  3:222932/1‑149 (MQ=255)
                                               cccTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACt                                                                                                    >  3:290849/1‑149 (MQ=255)
                                               cccTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTATCTCTAACTCAAAGACGGACCCCCTAAAAACt                                                                                                    >  1:140118/1‑149 (MQ=255)
                                                ccTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTc                                                                                                   >  2:138392/1‑149 (MQ=255)
                                                 cTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTc                                                                                                   >  5:219589/1‑148 (MQ=255)
                                                         gCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGa                                                                                          >  7:37261/1‑149 (MQ=255)
                                                          cATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGAc                                                                                         <  6:320034/149‑1 (MQ=255)
                                                                   tGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACACTTCctct                                                                                >  5:185117/1‑149 (MQ=255)
                                                                     ccACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCctct                                                                                <  8:77440/147‑1 (MQ=255)
                                                                     ccACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTgg                                                                              >  3:72515/1‑149 (MQ=255)
                                                                     ccACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTg                                                                               <  1:305449/148‑1 (MQ=255)
                                                                     ccACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTg                                                                               <  1:56718/148‑1 (MQ=255)
                                                                     ccACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTATCTCTATCTCAAAGACGGACCCCGTAAAAAAACCATCATGCACCCTTCAACTgg                                                                              >  2:295988/1‑149 (MQ=255)
                                                                      cACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGt                                                                             >  7:138893/1‑149 (MQ=255)
                                                                      cACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGt                                                                             <  6:95135/149‑1 (MQ=255)
                                                                       aCTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTa                                                                            >  7:197407/1‑149 (MQ=255)
                                                                        cTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTc                                                                                    <  3:159420/140‑1 (MQ=255)
                                                                        cTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTc                                                                                    >  4:159420/1‑140 (MQ=255)
                                                                           ggCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTAtt                                                                          <  8:170419/147‑1 (MQ=255)
                                                                           ggCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTGTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAg                                                                        >  7:170970/1‑149 (MQ=255)
                                                                            gCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGc                                                                       >  3:21223/1‑149 (MQ=255)
                                                                                gCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCCTAAAAACTCCATCATGGACCCTTCCTCTgg                                                                              >  2:29957/1‑138 (MQ=255)
                                                                                gCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCCTAAAAACTCCATCATGGACCCTTCCTCTgg                                                                              <  1:29957/138‑1 (MQ=255)
                                                                                       gTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGtt                                                            <  3:188545/149‑1 (MQ=255)
                                                                                        tCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTa                                                           >  5:205923/1‑149 (MQ=255)
                                                                                        tCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTa                                                           <  6:185117/149‑1 (MQ=255)
                                                                                                  tttGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTa                                                                         <  2:16607/125‑1 (MQ=255)
                                                                                                  tttGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTa                                                                         >  1:16607/1‑125 (MQ=255)
                                                                                                     gCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTc                                               <  1:138392/148‑1 (MQ=255)
                                                                                                      cTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCg                                              <  2:179560/148‑1 (MQ=255)
                                                                                                      cTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCg                                              <  8:138893/148‑1 (MQ=255)
                                                                                                       ttttACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCgg                                             >  4:255506/1‑148 (MQ=255)
                                                                                                         ttACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGc                                            <  4:72515/147‑1 (MQ=255)
                                                                                                         ttACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGc                                            <  4:21223/147‑1 (MQ=255)
                                                                                                          tACGGCCCCAGCCGCGACATTGGTTACCCTGTGCCAGAAAATGCTGGCTCTTTTTTTCTCTCTATGTCAAAGACGGACCACGTAAATACTACTTCATGGACACTTCCTCTGGTATTAGCCTCACATCGTTATTCGCCGACTTCGTCAt                                          <  4:13093/148‑1 (MQ=255)
                                                                                                            cGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGa                                        <  6:219589/148‑1 (MQ=255)
                                                                                                             ggCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGat                                       >  5:200216/1‑148 (MQ=255)
                                                                                                               ccccAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTATTTCTCTCTATCTCAAAGTCGGACGCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGatcat                                    <  8:37261/149‑1 (MQ=255)
                                                                                                                 ccAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCg                                              <  3:178301/137‑1 (MQ=255)
                                                                                                                 ccAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCg                                              >  4:178301/1‑137 (MQ=255)
                                                                                                                   aGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCttt                                <  3:92968/149‑1 (MQ=255)
                                                                                                                                gATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTgg                   >  8:123778/1‑149 (MQ=255)
                                                                                                                                gATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGACTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTgg                   >  1:233777/1‑149 (MQ=255)
                                                                                                                                    cccTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGtgct               >  6:107722/1‑149 (MQ=255)
                                                                                                                                     ccTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGtgctg              >  1:249997/1‑149 (MQ=255)
                                                                                                                                     ccTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGtgctg              >  6:138964/1‑149 (MQ=255)
                                                                                                                                     ccTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGtgctg              >  4:118093/1‑149 (MQ=255)
                                                                                                                                            cAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCTAAGACGGACCCCCTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGGCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGTGCTGCTCAACg       <  2:140118/149‑1 (MQ=255)
                                                                                                                                            cAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCTTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGTGCTGCTCAACg       <  4:290849/149‑1 (MQ=255)
                                                                                                                                            cAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCcatca                                                                                              <  8:160095/62‑1 (MQ=255)
                                                                                                                                            cAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCcatca                                                                                              >  7:160095/1‑62 (MQ=255)
                                                                                                                                            cAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGTGCTGCTCAACg       <  3:255506/149‑1 (MQ=255)
                                                                                                                                                aaaaGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGTGCTGCTCAACGGctt   <  8:197407/149‑1 (MQ=255)
                                                                                                                                                 aaaGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGTGCTGCTCAACGGcttc  >  4:201124/1‑149 (MQ=255)
                                                                                                                                                 aaaGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGTGCTGCTCAACGGcttc  >  5:124347/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGTGCTGCTCAACGGCTTC  >  NC_002947/6067157‑6067450

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: