Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,798,594 C→T 12.0% intergenic (‑185/+33) pvdH ← / ← PP_4224 diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,798,5940CT12.0% 535.5 / 49.1 205intergenic (‑185/+33)pvdH/PP_4224diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (100/78);  new base T (9/15);  total (112/93)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 208 total aligned reads displayed.

TGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGCCCC  >  NC_002947/4798445‑4798741
                                                                                                                                                     |                                                                                                                                                   
tGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTc                                                                                                                                                     <  2:5830934/150‑1 (MQ=255)
tGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGGAGGAGCGGGTTc                                                                                                                                                     >  1:3955824/1‑150 (MQ=255)
  gccctGGGTACCTCCTCGCGAGCAAGCGGAGGGGGGGTTATCGCTTTGGAAACGCGTGCGTGGGGGGAGGATTTAGCGGTTGGTCCGCGCTTTTATTTTTGTCTGTGCGGGCCTCTTCGCGGGTGCCCCCGCTCCTACGGGCGGGTTTAc                                                                                                                                                   <  2:8929254/146‑1 (MQ=19)
  gACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAc                                                                                                                                                   <  1:1718437/150‑1 (MQ=255)
  gACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAc                                                                                                                                                   <  1:3375617/150‑1 (MQ=255)
   aCATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGGTTTAcc                                                                                                                                                  <  1:5329798/150‑1 (MQ=38)
   aCATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGAGTAAACCCGCCCGTAGGAGCGGGTTCAcc                                                                                                                                                  <  2:2237208/150‑1 (MQ=255)
    cATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGGTTTAccc                                                                                                                                                 <  1:1423784/150‑1 (MQ=38)
     aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcg                                                                                                                                                >  2:2794014/1‑150 (MQ=255)
      tGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                               <  1:7460331/150‑1 (MQ=255)
      tGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                               <  1:3749776/150‑1 (MQ=255)
        ggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGa                                                                                                                                             >  1:79752/1‑150 (MQ=255)
         gTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGaa                                                                                                                                            >  1:5847025/1‑150 (MQ=255)
         gTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGaa                                                                                                                                            <  1:3899411/150‑1 (MQ=255)
           aaaTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCGTTGGAAACGCGTGAGTGGGGGGCGGATTTAGCGGTTGTTACGCGATGTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGGAGGAGCGGGTTCACCCGCGAagc                                                                                                                                          >  2:7418716/1‑149 (MQ=255)
                cTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCa                                                                                                                                     <  2:7225976/150‑1 (MQ=37)
                 tCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAg                                                                                                                                    >  1:6019118/1‑150 (MQ=255)
                 tCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAAACCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAg                                                                                                                                    >  1:2690208/1‑150 (MQ=255)
                        cAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGc                                                                                                                             >  2:1438442/1‑150 (MQ=255)
                        cAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGc                                                                                                                             <  2:4677595/150‑1 (MQ=255)
                         aaGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCg                                                                                                                            >  2:3420064/1‑150 (MQ=255)
                          aggggAGGGGGGGTTTTCGATTTGGAAACGCTTGTGTGGGGGGTGTTTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTCGGCGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGa                                                                                                                           <  1:8809814/147‑1 (MQ=37)
                            cGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATa                                                                                                                         <  1:2066942/150‑1 (MQ=255)
                              gATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGGAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATaaa                                                                                                                       >  2:1796557/1‑150 (MQ=255)
                               aTGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAAc                                                                                                                      <  2:5028177/150‑1 (MQ=255)
                                  ggCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATc                                                                                                                   >  1:8266761/1‑150 (MQ=255)
                                   gCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCa                                                                                                                  <  1:4536055/150‑1 (MQ=255)
                                     ggTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAAt                                                                                                                <  2:9134875/150‑1 (MQ=255)
                                     ggTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAAt                                                                                                                <  1:2898855/150‑1 (MQ=255)
                                          tggCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTAACGCTATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCa                                                                                                           <  1:7017221/148‑1 (MQ=255)
                                          tCGCTTTGTAACCGCTGGCGGGGGGGGTGTATTTAGCGTTTGTACCGCTATTTTATTTTGGTCTGTGCGGGCCTCTGCGGGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCTAAGAGGCCAGCACAGGCGATAAACATCAATGCCCa                                                                                                           <  1:7709495/150‑1 (MQ=25)
                                          tCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCa                                                                                                           <  2:1943669/150‑1 (MQ=255)
                                             cTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGa                                                                                                        <  2:2686015/150‑1 (MQ=255)
                                              tttGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATTCCCAGGAc                                                                                                       >  1:3621852/1‑150 (MQ=255)
                                              tttGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGAc                                                                                                       >  2:4752444/1‑150 (MQ=255)
                                               ttGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCCCGAAGAGGCCAGCACAGGCGATAAACATCAAATCCCAGGAcg                                                                                                      >  1:1016518/1‑150 (MQ=255)
                                                  gAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCCCCCGTAGGAGCGGGTTAACCCCCGAAGAGGCCAGCACACGCGATAAAAAAAAAATCCCAGGACgcgc                                                                                                   >  1:4180435/1‑150 (MQ=38)
                                                    aaCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCg                                                                                                 >  2:66041/1‑150 (MQ=255)
                                                    aaCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAAAGCCCAGGACGCGCCg                                                                                                 >  1:7294020/1‑150 (MQ=255)
                                                       gCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTa                                                                                              >  2:2125996/1‑150 (MQ=255)
                                                       gCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCACCACACGCGATAAAAATCAAATCCCACGACGCGCCCGTa                                                                                              >  1:968387/1‑150 (MQ=38)
                                                          tGAGTGGGGTGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCt                                                                                           <  2:7275501/150‑1 (MQ=255)
                                                           gAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTg                                                                                          >  2:3745355/1‑150 (MQ=255)
                                                             gTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGca                                                                                        >  2:4876589/1‑150 (MQ=255)
                                                             gTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACAACAATGCCCAGGACGCGCCCGTACCTGca                                                                                        >  1:4309873/1‑150 (MQ=255)
                                                              tGGGGGGATGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGcaa                                                                                       <  2:6815018/150‑1 (MQ=255)
                                                              tGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGCAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGcaa                                                                                       <  2:3940171/150‑1 (MQ=255)
                                                               ggggggAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGcaac                                                                                      >  2:77663/1‑150 (MQ=255)
                                                                gggggAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGcaaca                                                                                     >  2:4408337/1‑150 (MQ=255)
                                                                gggggAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCACGACGCGCCGGTACCTGcaaca                                                                                     >  2:1549707/1‑150 (MQ=255)
                                                                    gaggaTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAAc                                                                                 >  1:1506940/1‑150 (MQ=255)
                                                                          ttAGGGGTTGTTACGCTATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTTGGAGGGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCaa                                                                           <  1:872906/150‑1 (MQ=255)
                                                                          ttAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCaa                                                                           <  2:8541399/150‑1 (MQ=255)
                                                                             gCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGcc                                                                        <  2:6530624/150‑1 (MQ=255)
                                                                                gttgttACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGTTAATCATCAATGGCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTc                                                                     >  2:8427277/1‑150 (MQ=255)
                                                                                        gcgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                                                                                       >  2:8945359/1‑76 (MQ=255)
                                                                                        gcgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTCCCTGCAACAAAACCCGCAAGCCGTCACTACGGc                                                             >  2:2418003/1‑150 (MQ=255)
                                                                                         cgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCt                                                            <  2:3275508/150‑1 (MQ=255)
                                                                                            ttttattttTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCATGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTc                                                         <  2:8499531/150‑1 (MQ=255)
                                                                                            ttttattttTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGGAGGAGCGGGTGCACCCGCGACGAGGCCAGCACAGGCGCTAAACATCAATGCCCAGGACGCGCCGGTACCTGCACCAAAACCCGCACGCCGTCACTACGGCTggc                                                         >  2:8330728/1‑148 (MQ=38)
                                                                                                     tGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAg                                                <  2:3735985/150‑1 (MQ=255)
                                                                                                         tgtgCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGcac                                                                                                                                 >  1:6863393/1‑65 (MQ=255)
                                                                                                         tgtgCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTa                                            >  2:170340/1‑150 (MQ=255)
                                                                                                                  ctctTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCacac                                   >  1:2909471/1‑150 (MQ=255)
                                                                                                                   tctTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACAcc                                  >  2:6573964/1‑150 (MQ=255)
                                                                                                                   tacgggcgGGTTTACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACAcc                                  >  2:5593111/6‑150 (MQ=255)
                                                                                                                     ttCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGc                                >  1:3886993/1‑150 (MQ=255)
                                                                                                                        gcgGGTAAACCCGCCCGTTGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGc                             <  1:3088129/150‑1 (MQ=255)
                                                                                                                        gcgGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGc                             >  2:5950495/1‑150 (MQ=255)
                                                                                                                        gcgGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGc                             >  1:1154927/1‑150 (MQ=255)
                                                                                                                         cgtgTAATCGCGCCCGTCGGAGCGGGTTCACCCGCGAAGAGGCCAGCTCAGGCGATAAACATCATTGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCTGGCTACCGGCACACCGCTGCa                            <  2:8233178/147‑1 (MQ=255)
                                                                                                                         cgGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCa                            >  2:7248038/1‑150 (MQ=255)
                                                                                                                          gggTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGACGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc                           >  2:1760784/1‑150 (MQ=37)
                                                                                                                          gggTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc                           >  1:2850508/1‑150 (MQ=255)
                                                                                                                            gtgaaCCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCATAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGa                         >  1:819861/4‑150 (MQ=38)
                                                                                                                            gtgaaCCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGa                         >  1:9033345/4‑150 (MQ=38)
                                                                                                                            gTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACaaaa                                                                                  <  1:8809008/93‑1 (MQ=255)
                                                                                                                            gTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACaaaa                                                                                  >  2:8809008/1‑93 (MQ=255)
                                                                                                                            gTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGGCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGa                         >  2:7994320/1‑150 (MQ=255)
                                                                                                                             tgaaCCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGAt                        >  2:4060556/3‑150 (MQ=255)
                                                                                                                              gaaCCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATc                       >  2:8726443/2‑150 (MQ=38)
                                                                                                                              gaaCCCGCTCCTACGGGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATc                       >  2:4228843/2‑150 (MQ=255)
                                                                                                                               aaCCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACAATCg                      <  2:3886993/150‑1 (MQ=255)
                                                                                                                                aCCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGc                     >  1:1983733/1‑150 (MQ=255)
                                                                                                                                 cccgctccTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTc                                                         <  2:1188458/113‑1 (MQ=37)
                                                                                                                                 cccgctccTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGcc                    >  1:6609617/1‑150 (MQ=39)
                                                                                                                                 cccgctccTACGGGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGcc                    >  2:6658365/1‑150 (MQ=255)
                                                                                                                                 cccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGcc                    <  1:3734739/150‑1 (MQ=255)
                                                                                                                                  ccgcccgTAGGATAGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAATACGTAACTACGGCTGTCAAAGCGCATGCTACCGGCACACCGCTGCTCGATCGCCt                   >  2:5101076/1‑150 (MQ=255)
                                                                                                                                  ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGcc                    <  2:452139/149‑1 (MQ=255)
                                                                                                                                     cccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGca                >  2:2758263/1‑150 (MQ=255)
                                                                                                                                           ttcgCGGGTGAACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCg          >  1:7742130/4‑150 (MQ=255)
                                                                                                                                            gAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCATGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGc         <  2:4489225/150‑1 (MQ=255)
                                                                                                                                            gAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGa                         >  1:4435076/1‑134 (MQ=255)
                                                                                                                                            gAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGc         >  1:2346044/1‑150 (MQ=255)
                                                                                                                                            gAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGc         >  1:2346264/1‑150 (MQ=255)
                                                                                                                                             aGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGcc        >  2:6357046/1‑150 (MQ=255)
                                                                                                                                             aGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGcc        >  2:3319880/1‑150 (MQ=255)
                                                                                                                                                gggTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGc     <  1:3289027/150‑1 (MQ=255)
                                                                                                                                                   ttCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGcaac  >  1:7441344/1‑147 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGCCCC  >  NC_002947/4798445‑4798741

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: