Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,289,435 T→A 7.7% E1051V (GAG→GTG)  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,4350TA 7.7% 186.5 / 5.1 65E1051V (GAG→GTG) PP_3761two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base T (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 = 8.49e-01

CTGGCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGTGCAGGA  >  NC_002947/4289290‑4289582
                                                                                                                                                 |                                                                                                                                                   
cTGGCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTc                                                                                                                                                    >  6:365137/1‑147 (MQ=255)
 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                                                                                                                                <  4:240955/137‑1 (MQ=255)
                              gTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCaa                                                                                                                                >  3:240955/1‑137 (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                                                  <  3:141335/149‑1 (MQ=255)
                                                                                                gATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGa                                                  <  5:216485/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)
                                                                                                                                                 |                                                                                                                                                   
CTGGCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGTGCAGGA  >  NC_002947/4289290‑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: ‑