Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,289,437 C→A 7.7% L1050L (CTG→CTT)  PP_3761 ← two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,289,4370CA 7.7% 202.1 / 5.2 65L1050L (CTG→CTT) PP_3761two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base A (2/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

TGGCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGTGCAGGA  >  NC_002947/4289291‑4289582
                                                                                                                                                  |                                                                                                                                                 
tGGCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAg                                                                                                                                                 >  3:26234/1‑149 (MQ=255)
  gCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAggg                                                                                                                                               <  5:294899/149‑1 (MQ=255)
    ggATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGgcg                                                                                                                                             <  4:305598/149‑1 (MQ=255)
    ggATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTGGAGCACAAGGgcg                                                                                                                                             >  8:117948/1‑149 (MQ=255)
      aTCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCT‑CAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTg                                                                                                                                          <  8:177393/149‑1 (MQ=255)
       tCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTGGAGCACAAGGgcgcca                                                                                                                                          <  7:117948/149‑3 (MQ=255)
       gcAGGCGGGTGGCCTGGATCCCGTGCATCTGCGGCTTCATGACGTCCATCAGGAGCAGGTCGATTTCTTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTGGAGCACAAGGgcgcca                                                                                                                                          <  8:180474/148‑3 (MQ=255)
          ggCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCCGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGt                                                                                                                                        >  6:356802/1‑148 (MQ=255)
             gggTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAgg                                                                                                                                    <  5:228488/149‑1 (MQ=255)
                  gacTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTACAGGCGCTCGATGGCTTCGCGACCGTTGCTGCAGATTTCGACAATTGCGCCCTTGTGCTCCAGGGCGCTGGTCAGTGCaa                                                                                                                                >  6:122732/3‑148 (MQ=255)
                   ccTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCa                                                                                                                                 >  5:226668/1‑146 (MQ=255)
                        aaGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGGCAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATAtt                                                                                                                         <  1:342154/149‑1 (MQ=255)
                            cGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCTATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTTCGGCCGATTTCGACAATGTCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGc                                                                                                                     <  5:122732/149‑1 (MQ=255)
                            cGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGc                                                                                                                     <  2:197628/149‑1 (MQ=255)
                            cGTCCATCTCCGGCATCATCACGTCCATCAGAACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACg                                                                                                                      >  4:141335/1‑148 (MQ=255)
                             gTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCa                                                                                                                                                  >  8:217659/1‑119 (MQ=255)
                             gTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCa                                                                                                                                                  <  7:217659/119‑1 (MQ=255)
                             gTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCaaa                                                                                                                               <  6:113637/138‑1 (MQ=255)
                             gTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCaa                                                                                                                                >  3:240955/1‑137 (MQ=255)
                             gTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCaa                                                                                                                                <  4:240955/137‑1 (MQ=255)
                             gTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCCTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCaaa                                                                                                                               >  5:113637/1‑138 (MQ=255)
                                  tctcCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACAtcg                                                                                                               >  5:217837/1‑149 (MQ=255)
                                     ccGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACAtcgtc                                                                                                             <  4:26234/148‑1 (MQ=255)
                                      cGGCATCATCACGTCCATCAGCACCAGGTCGATATAGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCACTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGCTATTAAGCACAtcgtcg                                                                                                            >  6:53266/1‑148 (MQ=255)
                                            catcaCGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTccacca                                                                                                     >  5:151049/1‑149 (MQ=255)
                                              tcaCGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACcag                                                                                                    >  6:145753/1‑148 (MQ=255)
                                                 cGTCCATCAGCACCAGGTCGATATCGTCGTGCTTCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCTGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCagcagc                                                                                                <  5:53266/149‑1 (MQ=255)
                                                          gCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGAGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTccacc                                                                                                      >  5:283649/1‑134 (MQ=255)
                                                               aGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTc                                                                                  <  4:124455/149‑1 (MQ=255)
                                                                 gTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCg                                                                                 >  2:133229/1‑148 (MQ=255)
                                                                       atCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCt                                                                          >  5:68959/1‑149 (MQ=255)
                                                                        gcgtcgtGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGGACACCtt                                                                         <  8:188844/148‑1 (MQ=255)
                                                                          gtcgtGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTg                                                                        <  5:365137/148‑1 (MQ=255)
                                                                               gctgctCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCgg                                                                  >  8:18543/1‑149 (MQ=255)
                                                                                     ccAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTgcgc                                                            >  1:27670/1‑149 (MQ=255)
                                                                                      cAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTgcgc                                                            <  6:226668/148‑1 (MQ=255)
                                                                                        ggCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGcgcgcg                                                         <  5:145753/149‑1 (MQ=255)
                                                                                            cTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCt                                                     <  6:151049/149‑1 (MQ=255)
                                                                                             tCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCtt                                                    >  7:379433/1‑149 (MQ=255)
                                                                                               gATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGa                                                  <  5:216485/149‑1 (MQ=255)
                                                                                               gATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGa                                                  <  3:141335/149‑1 (MQ=255)
                                                                                               gATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGa                                                  <  6:68959/149‑1 (MQ=255)
                                                                                                aTGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGACCACCTTGTCGCGGCTGCGCGCGGTCTTGAg                                                 >  6:312885/1‑149 (MQ=255)
                                                                                                 tGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCACGCGGTCTTGAg                                                 >  8:141393/1‑148 (MQ=255)
                                                                                                   gCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCAt                                              >  8:6422/1‑149 (MQ=255)
                                                                                                     ttCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGa                                                  <  6:270999/143‑1 (MQ=255)
                                                                                                     ttCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGa                                                  >  5:270999/1‑143 (MQ=255)
                                                                                                     ttCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGc                                            >  3:385011/1‑149 (MQ=255)
                                                                                                          gcCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGa                                       >  2:227328/1‑149 (MQ=255)
                                                                                                                   ggCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTg                              <  1:227328/149‑1 (MQ=255)
                                                                                                                     ccGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGa                            <  7:6422/149‑1 (MQ=255)
                                                                                                                       gATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACa                          <  4:385011/149‑1 (MQ=255)
                                                                                                                       gATTTCGACAACGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGc                                                <  7:215713/127‑1 (MQ=255)
                                                                                                                       gATTTCGACAACGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGc                                                >  8:215713/1‑127 (MQ=255)
                                                                                                                        aTTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAg                         <  1:133229/149‑1 (MQ=255)
                                                                                                                              aCAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGa                   >  8:300962/1‑149 (MQ=255)
                                                                                                                                aaTGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGAtt                 >  6:267267/1‑149 (MQ=255)
                                                                                                                                aaTGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGAt                  >  2:362647/1‑148 (MQ=255)
                                                                                                                                 aTGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTc                <  8:379433/149‑1 (MQ=255)
                                                                                                                                 aTGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCACGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTc                <  7:141393/149‑1 (MQ=255)
                                                                                                                                  cagcgcCCTGGAGCACAAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTc                <  2:359003/146‑1 (MQ=255)
                                                                                                                                  cagcgcCCTGGAGCACAAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTc                >  1:359003/3‑148 (MQ=255)
                                                                                                                                         cTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCTTTGACGTTCGTTGGACAGCTGCGATTCGACCTtgt        >  4:224333/1‑149 (MQ=255)
                                                                                                                                          ttGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGACCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTtgtg       >  7:230611/1‑149 (MQ=255)
                                                                                                                                                tCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACATCTGCGATTCGCCCTTGTACAGGa  >  4:97162/1‑148 (MQ=255)
                                                                                                                                                  cAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTg                              >  8:132679/1‑118 (MQ=255)
                                                                                                                                                  cAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTg                              <  7:132679/118‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGGCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGTGCAGGA  >  NC_002947/4289291‑4289582

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