Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,291,085 G→T 5.8% P501Q (CCG→CAG)  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,291,0850GT 5.8% 269.0 / 6.2 86P501Q (CCG→CAG) PP_3761two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (42/39);  new base T (3/2);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

GGCGATATCGAGAATGTCGTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGTTCCGA  >  NC_002947/4290940‑4291230
                                                                                                                                                 |                                                                                                                                                     
ggCGATATCGAGAATGTCGTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTACTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAgcggcg                                                                                                                                                    >  7:67321/1‑149 (MQ=255)
        cGAGAATGTCGTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGTGCAGc                                                                                                                                            >  5:50313/1‑149 (MQ=255)
         gagaATGTCGTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGc                                                                                                                                            <  2:38949/148‑1 (MQ=255)
         gagaATGTCGTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCt                                                                                                                                           >  4:349919/1‑149 (MQ=255)
            aaTGTCGTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCTAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCTGCGCACGCCGCTCaa                                                                                                                                        <  3:251068/149‑2 (MQ=255)
            aaTGTCGTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCTGCGCACGCCGCTCaa                                                                                                                                        <  2:6552/149‑2 (MQ=255)
             aTGTCGTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATg                                                                                                                                       <  8:227670/149‑1 (MQ=255)
              tGTCGTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATgg                                                                                                                                      <  2:32896/149‑1 (MQ=255)
              tGTCGTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATgg                                                                                                                                      <  2:65141/149‑1 (MQ=255)
                  gTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGTGCAGCTCATGGGAc                                                                                                                                   <  6:50313/148‑1 (MQ=255)
                       tCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGTCCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTg                                                                                                                             >  1:290812/1‑149 (MQ=255)
                       tCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTg                                                                                                                             <  5:134906/149‑1 (MQ=255)
                       tCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTg                                                                                                                             >  2:164601/1‑149 (MQ=255)
                          ggTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGc                                                                                                                           >  7:355031/1‑148 (MQ=255)
                           gTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGc                                                                                                                           >  3:53944/1‑147 (MQ=255)
                                gCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAt                                                                                                                    >  5:317801/1‑149 (MQ=255)
                                 cAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAt‑‑t                                                                                                                 >  2:326707/1‑149 (MQ=255)
                                 cAGAGCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAt‑‑t                                                                                                                 <  8:81520/149‑1 (MQ=255)
                                   gATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCg                                                                                                               <  8:355031/149‑1 (MQ=255)
                                    aTCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCg                                                                                                               >  5:307843/1‑148 (MQ=255)
                                     tCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGa                                                                                                             <  1:164601/149‑1 (MQ=255)
                                     tCGCTGCCGGCCGAGTACATGGATTCGGCGAACGTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGGGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGa                                                                                                             <  4:230083/149‑1 (MQ=255)
                                      cGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGa                                                                                                             <  4:53944/148‑1 (MQ=255)
                                       gTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGAc                                                                                                            <  6:155903/148‑1 (MQ=255)
                                        ttGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCt                                                                                                                                           >  3:131156/1‑118 (MQ=255)
                                        ttGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCt                                                                                                                                           <  4:131156/118‑1 (MQ=255)
                                             ggccgAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACt                                                                                                     <  2:346773/149‑1 (MQ=255)
                                                cgAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑Gc                                                                                                <  3:349919/149‑1 (MQ=255)
                                                   gttgATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGTGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCg                                                                                             <  6:1699/146‑1 (MQ=255)
                                                           aTTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTg                                                                                      >  5:210519/1‑148 (MQ=255)
                                                            ttCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTg                                                                                      >  3:17646/1‑147 (MQ=255)
                                                            ttCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTg                                                                                      >  1:250972/1‑147 (MQ=255)
                                                              cGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAg                                                                                   <  4:245245/148‑1 (MQ=255)
                                                                gCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCACCAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGCCAGGTGGGCGAGGAAT‑‑TCGGACGTATCCTT‑‑GCTCGATCTCTGCAg                                                                                   >  8:75962/1‑146 (MQ=255)
                                                                gCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGgtc                                                                                >  1:230477/1‑149 (MQ=255)
                                                                gCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTGGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAg                                                                                   <  7:75962/146‑1 (MQ=255)
                                                                              ctcctcGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCt                                                                  >  1:222645/1‑149 (MQ=255)
                                                                              ctcctcGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTAATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGCTCGCTGAAGGTCGTCGGCACGGGCCt                                                                  >  1:228423/1‑149 (MQ=255)
                                                                                cctcGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCg                                                                                         >  6:51859/1‑124 (MQ=255)
                                                                                cctcGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCg                                                                                         <  5:51859/124‑1 (MQ=255)
                                                                                cctcGCTGAGGTTGCCTTCAGCGTTTGCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGGGCGCAGCGCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGc                                                                <  5:357259/149‑1 (MQ=255)
                                                                                 ctcGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTg                                                                 >  4:41129/1‑147 (MQ=255)
                                                                                  tcGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTg                                                                                                                             >  1:172315/1‑90 (MQ=255)
                                                                                  tcGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTg                                                                                                                             >  1:178118/1‑90 (MQ=255)
                                                                                  tcGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTg                                                                                                                             <  2:172315/90‑1 (MQ=255)
                                                                                  tcGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTg                                                                                                                             <  2:178118/90‑1 (MQ=255)
                                                                                    gCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCTTGCAGCt                                                            >  7:292799/1‑149 (MQ=255)
                                                                                    gCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGc                                                             <  4:137206/148‑1 (MQ=255)
                                                                                        aGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAg                                                        <  3:41129/149‑1 (MQ=255)
                                                                                           ttGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTg                                                     >  1:262546/1‑149 (MQ=255)
                                                                                              ccTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTgcgc                                                  <  6:210519/149‑1 (MQ=255)
                                                                                               cTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTGGAGCTgcagg                                                 <  6:317801/149‑4 (MQ=255)
                                                                                                      gtgttCTGCCAACAGCTGCGCCCCGCTCAGCGAGCTGCTGCGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGgctca                                          <  8:289365/146‑1 (MQ=255)
                                                                                                       ttttCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGgctcag                                         >  6:159694/1‑149 (MQ=255)
                                                                                                         ttCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTg                                                                                                                             <  5:120685/67‑1 (MQ=255)
                                                                                                         ttCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTg                                                                                                                             >  6:120685/1‑67 (MQ=255)
                                                                                                            tGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCAt                                    >  7:234446/1‑149 (MQ=255)
                                                                                                              ccAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAATAGGGGGCCTTTAACTTTCGGTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGgctcagct                                       >  1:315310/1‑148 (MQ=255)
                                                                                                              ccAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCg                                  >  3:272040/1‑149 (MQ=255)
                                                                                                              ccAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCg                                  >  6:218404/1‑149 (MQ=255)
                                                                                                                aaCAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGa                                                                                                             <  4:193203/74‑1 (MQ=255)
                                                                                                                aaCAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGa                                                                                                             >  3:193203/1‑74 (MQ=255)
                                                                                                                  cAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCt                              >  8:65727/1‑149 (MQ=255)
                                                                                                                    gCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTg                            <  2:262546/149‑1 (MQ=255)
                                                                                                                    gCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTg                            <  2:230477/149‑1 (MQ=255)
                                                                                                                          ccAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGATCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGGCGTCGGCACGGGCATGAAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTc                       >  8:256221/1‑148 (MQ=255)
                                                                                                                           cAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTc                       >  4:301918/1‑147 (MQ=255)
                                                                                                                             ggATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAgc                   <  4:17646/149‑1 (MQ=255)
                                                                                                                                  aGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCtt               <  7:65727/148‑1 (MQ=255)
                                                                                                                                   gCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTg                            <  2:363910/134‑1 (MQ=255)
                                                                                                                                   gCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTg                            >  1:363910/1‑134 (MQ=255)
                                                                                                                                   gCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCg             <  8:67321/149‑1 (MQ=255)
                                                                                                                                     gAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGt           <  8:234446/149‑1 (MQ=255)
                                                                                                                                      aGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGc                                                       >  5:101391/1‑104 (MQ=255)
                                                                                                                                      aGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGc                                                       <  6:101391/104‑1 (MQ=255)
                                                                                                                                       gCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCATGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGc         >  5:367411/1‑149 (MQ=255)
                                                                                                                                         ccatGAGCTGCGCACGCCGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGt       >  7:283718/4‑149 (MQ=255)
                                                                                                                                         ccatGAGCTGCGCACGCCGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGt       >  6:148958/4‑149 (MQ=255)
                                                                                                                                         ccatGAGCTGCGCACGCCGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGt       >  2:237042/4‑149 (MQ=255)
                                                                                                                                          gTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGtt      >  8:227230/1‑149 (MQ=255)
                                                                                                                                          gTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGtt      >  5:163304/1‑149 (MQ=255)
                                                                                                                                           ttGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGTTc     <  5:218404/149‑1 (MQ=255)
                                                                                                                                           ttGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGTTc     >  3:188793/1‑149 (MQ=255)
                                                                                                                                            cgAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGTTcc    <  7:7603/148‑1 (MQ=255)
                                                                                                                                             gAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGTTcc    >  8:95596/1‑148 (MQ=255)
                                                                                                                                              aGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGTTCCGa  >  8:357227/1‑149 (MQ=255)
                                                                                                                                              aGCGGCGGGCGGAGCTCCTGGGACATGTTGGCGCGGAAT‑‑TCGGACTGGTTCGT‑‑GCTCGCTCGCTGCAGGTCGTGGGCGCGGGCCTGCAGCTCCAGCTGCGCGTGGCTCAGCTCCTCGGGCGGGGGGGCCAGCGGGTCGGTGGGttgg    >  8:7603/1‑145 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GGCGATATCGAGAATGTCGTTGATCAGGTTGAGCAGATCGTTGCCGGCCGAGTAGATGGATTCGGCGAACTTGACCTGCTCCTCGCTGAGGTTGCCTTCAGCGTTTTCTGCCAACAGCTTGGCCAGGATCAGCGAGCTGTTGAGCGGCGTGCGCAGCTCATGGGACATGTTGGCGAGGAAT‑‑TCGGACTTATACTT‑‑GCTCGATCGCTGCAGGTCGTCGGCACGGGCCTGCAGCTCCAGCTGCGCCTGGCTCAGCTCATCGTTCTTGCGGTCCAGCGCTTCGGTGCGTTCCGA  >  NC_002947/4290940‑4291230

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