Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,072,625 G→A 10.1% A111T (GCC→ACC) ‡ PP_2683 → 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_0029473,072,6250GA10.1% 159.6 / 7.2 60A111T (GCC→ACC) ‡PP_2683two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (3/3);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

TCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTG  >  NC_002947/3072477‑3072771
                                                                                                                                                    |                                                                                                                                                  
ggttggcGGCACGCATGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCg                                                                                                                                                    >  4:96982/7‑149 (MQ=255)
 cAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCg                                                                                                                                                    <  4:377283/148‑1 (MQ=255)
    gACGGCATGCGTGCGGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCg                                                                                                                                                >  6:103834/1‑149 (MQ=255)
      cGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGt                                                                                                                                               >  1:48592/1‑148 (MQ=255)
       ggCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGa                                                                                                                                             <  4:210290/149‑1 (MQ=255)
       ggCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGa                                                                                                                                             <  6:152894/149‑1 (MQ=255)
       ggCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGa                                                                                                                                             <  2:48592/149‑1 (MQ=255)
           tGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGcac                                                                                                                                         <  2:303094/149‑1 (MQ=255)
            gCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGca                                                                                                                                          >  8:26330/1‑147 (MQ=255)
            gCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCACGGGCGAGcacc                                                                                                                                        >  5:71848/1‑148 (MQ=255)
            gCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCATGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCTGTGCTCACGGGCGAGcacc                                                                                                                                        >  8:287702/1‑148 (MQ=255)
             cGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCACGGGCGAGcaccg                                                                                                                                       <  6:71848/149‑3 (MQ=255)
              gTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGc                                                                                                                                      <  2:171294/149‑1 (MQ=255)
               tGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCATGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCTGTGCTCACGGGCGAGCACCGcg                                                                                                                                     <  7:287702/149‑2 (MQ=255)
                       aCCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGc                                                                                                                             <  2:174509/149‑1 (MQ=255)
                               cTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAAc                                                                                                                     <  1:201090/149‑1 (MQ=255)
                                      gcATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTg                                                                                                              <  1:264569/149‑1 (MQ=255)
                                         tGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGc                                                                                                                                      >  2:16285/1‑122 (MQ=255)
                                         tGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGc                                                                                                                                      <  1:16285/122‑1 (MQ=255)
                                         tGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTgcgg                                                                                                           >  3:268523/1‑149 (MQ=255)
                                            tCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCACGGGCGAGCACCGc                                                                                                                                      >  5:25118/1‑119 (MQ=255)
                                            tCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCACGGGCGAGCACCGc                                                                                                                                      <  6:25118/119‑1 (MQ=255)
                                                cTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAg                                                                                                     >  4:125518/1‑148 (MQ=255)
                                                    gacgacCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACgg                                                                                                >  2:20946/1‑149 (MQ=255)
                                                    gacgacCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACgg                                                                                                >  5:102270/1‑149 (MQ=255)
                                                           cGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCAt                                                                                         >  7:251191/1‑149 (MQ=255)
                                                            gCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCAtc                                                                                        >  4:130708/1‑149 (MQ=255)
                                                            gCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCAtc                                                                                        >  3:14050/1‑149 (MQ=255)
                                                                 cGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGAt                                                                                   <  1:378495/149‑1 (MQ=255)
                                                                 cGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGAt                                                                                   <  4:78145/149‑1 (MQ=255)
                                                                    tGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGAtgtg                                                                                >  4:272936/1‑149 (MQ=255)
                                                                      cTGTTTTCCCTGACTGACCTGGCCGCAAACCTTTTTGATGGCGGTGCAGAGGAGTTGCAGGGCATTACCGCGGTGCTCGCCCGTGAGCACAGACTGCACGGCTACGAAACCCGTCTTCGGCGCAAGGACGGCAGCCATCTCGAtgttc                                                                               <  1:223292/148‑3 (MQ=255)
                                                                            tCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGaa                                                                        >  5:253630/1‑149 (MQ=255)
                                                                              ccTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAAc                                                                       <  3:125518/148‑1 (MQ=255)
                                                                               ctgactgaCCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACct                                                                     >  7:162812/1‑149 (MQ=255)
                                                                                      acCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCtgctgc                                                                <  6:102270/147‑1 (MQ=255)
                                                                                      acCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTgctgct                                                               <  5:103834/148‑1 (MQ=255)
                                                                                         tGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCaaa                                                           >  5:360349/1‑149 (MQ=255)
                                                                                          ggCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAAc                                                          >  2:336253/1‑149 (MQ=255)
                                                                                           gccgccAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCaaa                                                           <  3:130708/147‑1 (MQ=255)
                                                                                             cgccAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCg                                                       >  5:54290/1‑149 (MQ=255)
                                                                                              gccAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCgg                                                      >  3:133788/1‑149 (MQ=255)
                                                                                                   ccTGTTTGATGGCGGTGCAGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATg                                                 <  3:280203/149‑1 (MQ=255)
                                                                                                     tGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATg                                                 >  3:201752/1‑147 (MQ=255)
                                                                                                       tttGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGaa                                               <  3:96982/147‑1 (MQ=255)
                                                                                                       tttGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAg                                              <  7:26330/148‑1 (MQ=255)
                                                                                                           aTGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCtg                                         >  2:170196/1‑149 (MQ=255)
                                                                                                           aTGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCt                                          <  6:360349/148‑1 (MQ=255)
                                                                                                                 gTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAg                                   <  4:86447/149‑1 (MQ=255)
                                                                                                                  tGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACg                                                                                                 >  5:15476/1‑86 (MQ=255)
                                                                                                                  tGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACg                                                                                                 <  6:15476/86‑1 (MQ=255)
                                                                                                                     cGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGt                               <  4:268523/149‑1 (MQ=255)
                                                                                                                          agTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTc                          <  1:20946/149‑1 (MQ=255)
                                                                                                                            ttGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAg                                              <  7:151932/127‑1 (MQ=255)
                                                                                                                            ttGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAg                                              >  8:151932/1‑127 (MQ=255)
                                                                                                                               cAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCg                      <  1:336253/148‑1 (MQ=255)
                                                                                                                                  gCCATTACGGGGGTGCTCTGCGTTGAGAACAGCCTGGACTGCTAGGAAACCCGTTTGCGGAGAAAGGACTGCAGCCATCTCTATGTGCTGATGAACCTGCTGATCAAACCGGGCCATGATGGGCTGGTGTAGGGGTTTGTCGCCGACAt                  <  4:133788/149‑1 (MQ=255)
                                                                                                                                     aTTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGTTGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTAc               <  5:266165/149‑1 (MQ=255)
                                                                                                                                     aTTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGttt                             <  4:309222/135‑1 (MQ=255)
                                                                                                                                     aTTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGttt                             >  3:309222/1‑135 (MQ=255)
                                                                                                                                         ccGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGa            <  1:170196/148‑1 (MQ=255)
                                                                                                                                                  tCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTg  >  3:332513/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTTTGATGGCGGTGCCGAGGAGTTGCAGGCCATTACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTG  >  NC_002947/3072477‑3072771

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