Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,067,301 A→T 6.4% intergenic (+117/‑57) 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,3010AT 6.4% 234.4 / 4.3 78intergenic (+117/‑57)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 = 9.39e-01

GCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGTGCTGCTCAACGGCTT  >  NC_002947/6067154‑6067449
                                                                                                                                                   |                                                                                                                                                    
gcCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCa                                                                                                                                                      >  8:223939/1‑148 (MQ=255)
   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                                                                                                                                          <  1:205873/147‑2 (MQ=255)
             gCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCTTTTGCTGGCAca                                                                                                                                          >  2:205873/1‑146 (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                                                                                                                                 >  6:254885/1‑148 (MQ=255)
                     aaaGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTtctc                                                                                                                                 <  5:254885/148‑1 (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                                                                            <  6:95135/149‑1 (MQ=255)
                                                                         cACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGt                                                                            >  7:138893/1‑149 (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:21223/147‑1 (MQ=255)
                                                                                                            ttACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGc                                           <  4:72515/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             >  4:118093/1‑149 (MQ=255)
                                                                                                                                        ccTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGtgctg             >  6:138964/1‑149 (MQ=255)
                                                                                                                                        ccTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGtgctg             >  1:249997/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)
                                                                                                                                                   |                                                                                                                                                    
GCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTAAAAACTCCATCATGGACCCTTCCTCTGGTATTAGCCTCACCTCGTTATTCGCCGACTTCGGCATGATCATCTTTGCCTTCCTGCTGGTGCTGCTCAACGGCTT  >  NC_002947/6067154‑6067449

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