Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 723,608 T→G 5.8% A269A (GCA→GCC kdpD ← fused sensory histidine kinase in two‑component regulatory system with KdpE: signal sensing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913723,6080TG 5.8% 247.8 / 4.1 90A269A (GCA→GCCkdpDfused sensory histidine kinase in two‑component regulatory system with KdpE: signal sensing protein
Reads supporting (aligned to +/- strand):  ref base T (28/56);  new base G (2/4);  total (30/60)
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.94e-01

AGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGGCGCAGTGCC  >  NC_000913/723459‑723752
                                                                                                                                                     |                                                                                                                                                
cgtgtTGCCGGCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCt                                                                                                                                                  <  1:989911/149‑1 (MQ=255)
   gtTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTgcg                                                                                                                                               >  2:398286/1‑150 (MQ=255)
   gtTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTgcg                                                                                                                                               >  2:1055430/1‑150 (MQ=255)
      gCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCgg                                                                                                                                            >  2:317941/1‑150 (MQ=255)
      gCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCgg                                                                                                                                            >  2:1285112/1‑150 (MQ=255)
      gCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCgg                                                                                                                                            <  2:1572402/150‑1 (MQ=255)
        cGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAACCGGTGCAGGGCAGGGGTTTCAACATACACAGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGAc                                                                                                                                          >  2:733961/1‑150 (MQ=255)
         gTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTATCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGAcc                                                                                                                                         <  2:1188691/150‑1 (MQ=255)
          gctcCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCAGCGGCGCggctgg                                                                                                                                        <  1:855594/149‑5 (MQ=255)
          gctcCGCGCCCAGTGCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTGTTTTTTCCGGGACGCGGTGCAGGGCAGGGGTTGCAACATCCACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCAGCGGCGCggctgg                                                                                                                                        <  1:281007/149‑5 (MQ=255)
           ctcCGCGCCCCGTTCCTGCGACAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGAAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCAGCGGCGCggctggc                                                                                                                                       <  2:517361/150‑6 (MQ=255)
                cgcCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTCAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGttttt                                                                                                                                  <  1:1142893/150‑1 (MQ=255)
                 gcCCAGTTCCTGCGGCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTc                                                                                                                                 <  1:333285/150‑1 (MQ=255)
                 gcCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTc                                                                                                                                 <  1:1307057/150‑1 (MQ=255)
                 gcCCAGTGCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTc                                                                                                                                 <  1:1699023/150‑1 (MQ=255)
                  ccccGTCCCCGCGCCAGACGTAATGCGCTGAGAATTGCCCGACGTGTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGGGCCAGACGCTACCCAGACGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCg                                                                                                                                <  2:325511/150‑1 (MQ=255)
                  cccAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCg                                                                                                                                >  2:1620828/1‑150 (MQ=255)
                  cccAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCg                                                                                                                                >  2:444105/1‑150 (MQ=255)
                      gTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCCGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGc                                                                                                                            >  2:407890/1‑150 (MQ=255)
                      gTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAACCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCTCCGCTGCGCGGACAAGTTTTTCGCtcc                                                                                                                            >  1:1515467/1‑148 (MQ=255)
                       ttCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGcc                                                                                                                           <  2:781094/150‑1 (MQ=255)
                          cTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGgt                                                                                                                        <  2:1633196/150‑1 (MQ=255)
                           tGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGgtg                                                                                                                       >  2:439492/1‑150 (MQ=255)
                                 cgACTTATGGCGCTGAGAATTGGCCGACCTTTTTTTTCCGGTAAGCGGTGGAGTGCCGGGGTTTCAACATAGACCGCGTGCCAGACGCTACCCAGCCGTGACGGCAGCCGCGCCGCTGCGCGGACCAGATTTTCGCTGCCGGTGTTATgg                                                                                                                 <  2:852196/149‑1 (MQ=255)
                                   aCGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGcc                                                                                                               <  2:1353065/150‑1 (MQ=255)
                                       aaGGCGCTGATAATTGCCCGACGGTTTTTTTCCGGTAAGCGGTGCCGGGCAGGGGGTTCAACCTTCCCCGGGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATg                                                                                                           <  1:1294404/150‑1 (MQ=255)
                                         ggCGCTGAGAATTGCCCGACGTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCAt                                                                                                        <  1:573029/150‑1 (MQ=255)
                                               gagaATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAg                                                                                                   <  1:585209/150‑1 (MQ=255)
                                               gagaATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAg                                                                                                   <  2:1529406/150‑1 (MQ=255)
                                                  aaTTGCCCGACGTTTTTTTTCCGGTAAGCGGGGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAt                                                                                                <  2:1104911/150‑1 (MQ=255)
                                                   atggCCCGCCGGGTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGCTGCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCAGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAtc                                                                                               <  1:65055/147‑1 (MQ=255)
                                                   agtGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGGGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAtc                                                                                               <  2:589505/148‑1 (MQ=255)
                                                   aTTGCCGGCCTTTTTTTTCCCGGAATGCGGGGAAGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAtcg                                                                                              <  1:1593290/150‑1 (MQ=255)
                                                      ccccGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAtcgcg                                                                                            <  2:121846/149‑1 (MQ=255)
                                                         cGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAtcgcgtcg                                                                                         <  2:939237/150‑1 (MQ=255)
                                                            cGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTcgcgc                                                                                      >  1:1556717/1‑150 (MQ=255)
                                                               tttttttCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCgtg                                                                                   <  2:1571460/150‑1 (MQ=255)
                                                                ttttttCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCgtgt                                                                                  <  1:444105/150‑1 (MQ=255)
                                                                ttttttCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCgtgt                                                                                  <  1:127081/150‑1 (MQ=255)
                                                                 tttttCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTATCGCTGCCGCTGTTATGGCCGCTGCATAAAAGCATCGCGTCGCGAGtgtg                                                                                 >  2:1129252/1‑150 (MQ=255)
                                                                    ttCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCca                                                                              <  1:975025/150‑1 (MQ=255)
                                                                       cGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCacac                                                                           <  1:1622482/150‑1 (MQ=255)
                                                                       cGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCacac                                                                           <  2:1682136/150‑1 (MQ=255)
                                                                       cGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTATCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCacac                                                                           >  2:1469001/1‑150 (MQ=255)
                                                                        ggTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGACGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCCCACt                                                                          >  1:886462/1‑150 (MQ=255)
                                                                          tAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACttt                                                                        <  2:1464386/150‑1 (MQ=255)
                                                                          tAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAGGATCGCGTCGCGCGTGTGCCACACTTtc                                                                       <  1:670767/150‑1 (MQ=255)
                                                                           aaGCGGTGCAGGGCAAGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTtc                                                                       <  2:1602539/150‑1 (MQ=255)
                                                                                gTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTc                                                                  >  1:1132606/1‑150 (MQ=255)
                                                                                   cagggcagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGCCCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGcc                                                               >  2:1002894/1‑150 (MQ=255)
                                                                                    agggcagggGTTTCACCATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTGATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCa                                                              <  2:42611/150‑1 (MQ=255)
                                                                                     gggcagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAg                                                             <  1:488886/150‑1 (MQ=255)
                                                                                     gggcagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGATTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAg                                                             >  1:644538/1‑150 (MQ=255)
                                                                                      ggcagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAgg                                                            >  1:1551129/1‑150 (MQ=255)
                                                                                       gcagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGa                                                           <  1:1176850/150‑1 (MQ=255)
                                                                                        cagggGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGATATGGCCGGTGCTTAAAAGGATCGCGTCGCGCGTGTGCCACACGTTCTCGGCGCCAggcg                                                          >  2:275004/1‑148 (MQ=255)
                                                                                             gTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGGTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGcccc                                                     <  1:72367/150‑1 (MQ=255)
                                                                                               ttCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCg                                                   <  2:1485036/150‑1 (MQ=255)
                                                                                                tCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGCTCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGCTGCCCCCGc                                                  >  2:1170751/1‑150 (MQ=255)
                                                                                                 cAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACGTTCTGTTCGCCAGGATGCCCCCGcc                                                 >  1:65098/1‑150 (MQ=255)
                                                                                                  aaCATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCAGAAAAGGATCGCGTCGCGCGGGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCa                                                >  1:1177816/1‑150 (MQ=255)
                                                                                                      tACTCCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGgcg                                            <  2:1321172/150‑1 (MQ=255)
                                                                                                      tACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGgcg                                            >  1:661130/1‑150 (MQ=255)
                                                                                                         acCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGgcgcgc                                         <  2:1037237/150‑1 (MQ=255)
                                                                                                          ccGCGTGCCAGACGCTACCCAGCCGGGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCGGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCa                                        <  1:25943/150‑1 (MQ=255)
                                                                                                              gTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTg                                    <  1:957107/150‑1 (MQ=255)
                                                                                                                gCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCt                                  <  1:544002/150‑1 (MQ=255)
                                                                                                                 ccAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCtc                                 <  2:63063/150‑1 (MQ=255)
                                                                                                                 ccAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCtc                                 <  1:1385753/150‑1 (MQ=255)
                                                                                                                 ccAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCtc                                 >  2:1431526/1‑150 (MQ=255)
                                                                                                                  cAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCTGTGTTATGGCCGATGCATAAAAGTATCTCGTCGCGCGTTTGCCACACTTTCTCTTCTCCAGGATGCCCCCGCCAGGCGCGCATTTGCtca                                <  1:1225653/150‑1 (MQ=255)
                                                                                                                   aGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCtcat                               <  1:567167/150‑1 (MQ=255)
                                                                                                                    gACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCtcatc                              <  1:1451257/150‑1 (MQ=255)
                                                                                                                        cTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAAcg                          <  1:282479/150‑1 (MQ=255)
                                                                                                                        atACACAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTATCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCAAGGATGCCCCCGACATTAGCTCATTTGCTCATCAAcg                          <  2:102221/149‑1 (MQ=255)
                                                                                                                         tACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAAcgc                         <  1:1316883/150‑1 (MQ=255)
                                                                                                                         tACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAAcgc                         <  1:946190/150‑1 (MQ=255)
                                                                                                                         tACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAAcgc                         <  1:1415313/150‑1 (MQ=255)
                                                                                                                             ccgCCGCGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGAGCCCGGCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATc                     <  2:1816217/148‑1 (MQ=255)
                                                                                                                             cAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATc                     <  2:1057672/150‑1 (MQ=255)
                                                                                                                                   ggACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGt               <  1:1528518/149‑1 (MQ=255)
                                                                                                                                      cGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACg            >  2:1763965/1‑150 (MQ=255)
                                                                                                                                      actggtccgcgcagcGGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACg            >  2:770235/7‑150 (MQ=255)
                                                                                                                                         cAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGgcg         >  2:953257/1‑150 (MQ=255)
                                                                                                                                         cAGCCGCGCCGCTGCGCGGACCAGATTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGgcg         >  2:1773682/1‑150 (MQ=255)
                                                                                                                                          gtccgcgcagcGGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCCCACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCTGTACGgcgc        >  2:1362996/3‑150 (MQ=255)
                                                                                                                                            ccgcgccgcTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACCCTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCAGTTGCTCATCAACGCGATCGGCAGTACGGCGCAg      >  1:1119837/1‑150 (MQ=255)
                                                                                                                                            ccgcgccgcTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGGCGCAg      <  2:801430/150‑1 (MQ=255)
                                                                                                                                            ccgcgccgcTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGAGGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGGCGCAg      <  1:771403/150‑1 (MQ=255)
                                                                                                                                              gcgccgcTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGGTGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGGCGCAGTg    <  2:1132606/150‑1 (MQ=255)
                                                                                                                                                gccgcGGCGCGGACCAGTCGTTCGCTGCCGGTGTTATGGCAGAAGCATAAAAGGATCGCGGCGCGCGTGTGCCCCACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGGCGCAGTGcc  <  2:259848/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AGTGTTGCCGTCTCCGCGCCCAGTTCCTGCGCCAGACGTAAGGCGCTGAGAATTGCCCGACGTTTTTTTTCCGGTAAGCGGTGCAGGGCAGGGGTTTCAACATACACCGCGTGCCAGACGCTACCCAGCCGTGACGCCAGCCGCGCCGCTGCGCGGACCAGTTTTTCGCTGCCGGTGTTATGGCCGATGCATAAAAGGATCGCGTCGCGCGTGTGCCACACTTTCTCTTCGCCAGGATGCCCCCGCCAGGCGCGCATTTGCTCATCAACGCGATCGGCAGTACGGCGCAGTGCC  >  NC_000913/723459‑723752

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: