Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,346,117 G→C 5.6% A128G (GCC→GGC)  PP_2942 ← two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,346,1170GC 5.6% 233.4 / 4.4 71A128G (GCC→GGC) PP_2942two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base C (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

AGCATGTCAGCGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGC  >  NC_002947/3345970‑3346260
                                                                                                                                                     |                                                                                                                                               
aGCATGTCAGCGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGc                                                                                                                                                >  5:312412/1‑149 (MQ=255)
  cATGTCAGCGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGCCGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCt                                                                                                                                              >  2:267577/1‑149 (MQ=255)
  cATGTCAGCGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCt                                                                                                                                              >  4:118081/1‑149 (MQ=255)
    tGTCAGCGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTACTCCAG‑‑GGCCTgg                                                                                                                                            <  5:78495/149‑1 (MQ=255)
       cAGCGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGgcg                                                                                                                                          >  7:98931/1‑148 (MQ=255)
          cGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAg                                                                                                                                      <  3:51850/149‑1 (MQ=255)
            gTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGt                                                                                                                                    <  8:98931/149‑1 (MQ=255)
             tCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTc                                                                                                                                   <  2:239489/149‑1 (MQ=255)
                 ctgcAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCAtt                                                                                                                                <  8:244451/148‑1 (MQ=255)
                  tgcAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCt                                                                                                                                              >  2:203112/1‑133 (MQ=255)
                  tgcAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCt                                                                                                                                              <  1:203112/133‑1 (MQ=255)
                  tgcAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAg                                                                                                                              >  4:226046/1‑149 (MQ=255)
                   gcAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCAt                                                                                                                                 >  2:324495/1‑145 (MQ=255)
                   gcAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCAt                                                                                                                                 <  1:324495/145‑1 (MQ=255)
                    cAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCTTTAGc                                                                                                                             <  3:101993/148‑1 (MQ=255)
                       cTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTACAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAAtt                                                                                                                         >  5:130537/1‑149 (MQ=255)
                        tCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATtc                                                                                                                        >  3:36397/1‑149 (MQ=255)
                           gcgGTACGGGCTGCAACGCGTTTCTCGAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATtcgtc                                                                                                                     <  7:262494/149‑1 (MQ=255)
                              gTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATtcgtcgtc                                                                                                                  <  3:222542/149‑1 (MQ=255)
                                  gggCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTc                                                                                                              >  4:275154/1‑149 (MQ=255)
                                  gggATGCAACGCGTTTCTCAAAGTCGTTATTTAGCAGTCTCAACTGGGCTTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCTTAAAGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTc                                                                                                              <  8:204810/149‑1 (MQ=255)
                                                 tctcCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGAt                                                                                                <  8:217948/148‑1 (MQ=255)
                                                   tcCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGt                                                                                              >  5:193947/1‑148 (MQ=255)
                                                     cAAGTCGTTATTGAGCAGTCTCAAATGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAACGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAAc                                                                                           >  7:306675/1‑149 (MQ=255)
                                                        gTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGt                                                                                        <  6:193947/149‑1 (MQ=255)
                                                         tCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTg                                                                                       <  2:243091/149‑1 (MQ=255)
                                                           gTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTgg                                                                                      <  4:288839/148‑1 (MQ=255)
                                                           gTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGa                                                                                     >  5:75329/1‑149 (MQ=255)
                                                           gTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGa                                                                                     >  6:74023/1‑149 (MQ=255)
                                                            ttATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAGTCGGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTgg                                                                                      <  1:267577/149‑1 (MQ=255)
                                                              aTTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCGGGGCTTACTGATGTAACGGTGGATCt                                                                                  >  4:56923/1‑149 (MQ=255)
                                                                  aGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCAt                                                                                                                                 <  3:291587/98‑1 (MQ=255)
                                                                  aGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCAt                                                                                                                                 >  4:291587/1‑98 (MQ=255)
                                                                       tctcAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCAt                                                                         >  1:249315/1‑149 (MQ=255)
                                                                        ctcAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCAtt                                                                        <  5:254151/149‑1 (MQ=255)
                                                                           aaCTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGa                                                                      <  3:275154/148‑1 (MQ=255)
                                                                               gggCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCa                                                                                                            >  5:142432/1‑106 (MQ=255)
                                                                               gggCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCa                                                                                                            <  6:142432/106‑1 (MQ=255)
                                                                                ggCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGcc                                                                 <  3:128321/148‑1 (MQ=255)
                                                                                   gTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGCCGGCCTtga                                                             >  8:56770/1‑149 (MQ=255)
                                                                                   gTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTtga                                                             <  4:36397/149‑1 (MQ=255)
                                                                                   gTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTtga                                                             <  3:226046/149‑1 (MQ=255)
                                                                                    ttCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGCTCTGACCGTCATTGACGGCCTtgat                                                            >  5:330810/1‑149 (MQ=255)
                                                                                                 ccAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTgg                                                <  6:75329/148‑1 (MQ=255)
                                                                                                     cagcTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCgg                                           >  7:140906/1‑149 (MQ=255)
                                                                                                         tCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGAGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTAATGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCa                                         >  6:291054/1‑147 (MQ=255)
                                                                                                             ggCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTAACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCg                                   >  1:332700/1‑149 (MQ=255)
                                                                                                                    gcgcTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCa                            <  6:312412/149‑1 (MQ=255)
                                                                                                                     ctcTCGCGTTCGGAATGCTGGTTCTCCAG‑‑GGCCTGGCGCAACGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTTCGCTTGTTTCGGCATAGCCGGTCAGCaa                           <  5:291054/147‑1 (MQ=255)
                                                                                                                      gcTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACTGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGt                                              <  6:285230/129‑1 (MQ=255)
                                                                                                                      gcTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGt                                              >  5:285230/1‑129 (MQ=255)
                                                                                                                        tCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGc                             >  1:193466/1‑144 (MQ=255)
                                                                                                                        tCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGc                             <  4:146025/144‑1 (MQ=255)
                                                                                                                        tCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGc                             >  3:146025/1‑144 (MQ=255)
                                                                                                                          gcgTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATcc                      >  4:172049/1‑149 (MQ=255)
                                                                                                                            gTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCTCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCt                     <  8:306675/148‑1 (MQ=255)
                                                                                                                               cGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGAc                  <  6:130537/148‑1 (MQ=255)
                                                                                                                                     gCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGAt           >  8:58227/1‑149 (MQ=255)
                                                                                                                                       tGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATgg         >  3:305993/1‑149 (MQ=255)
                                                                                                                                         ttGCGCCAG‑‑GCCCTGGAGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGt       >  4:263149/2‑149 (MQ=255)
                                                                                                                                         ttGCGCCAG‑‑GCCCTGGAGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGt       >  4:258680/2‑149 (MQ=255)
                                                                                                                                         gTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATgg         <  8:140906/147‑1 (MQ=255)
                                                                                                                                         gTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATgg         <  2:193466/147‑1 (MQ=255)
                                                                                                                                         gTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGt       >  5:45435/1‑149 (MQ=255)
                                                                                                                                         gTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGt       <  7:58227/149‑1 (MQ=255)
                                                                                                                                            cgccAG‑‑GCCCTGGAGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGa     <  3:263149/146‑1 (MQ=255)
                                                                                                                                            cgccAG‑‑GCCCTGGAGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGa     <  3:258680/146‑1 (MQ=255)
                                                                                                                                             tCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGAt    >  7:220720/1‑148 (MQ=255)
                                                                                                                                              ccAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGc  >  8:99177/1‑149 (MQ=255)
                                                                                                                                              ccAG‑‑GGCCTGGCACAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGa     <  2:128393/146‑1 (MQ=255)
                                                                                                                                              ccAG‑‑GGCCTGGCACAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGa     >  1:128393/1‑146 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AGCATGTCAGCGGTCTGCTGCAGCTCCGCGGTACGGGCTGCAACGCGTTTCTCCAAGTCGTTATTGAGCAGTCTCAACTGGGCGTTCTGCCGACGTGCCAGCAGCTCCAGGCGTTGGCGCTCGCGTTCGGAATGCTGGTGCTCCAG‑‑GGCCTGGCGCAAGGTCATTAGCAATTCGTCGTCGTTCCAGGGCTTACTGATGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGC  >  NC_002947/3345970‑3346260

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: