Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 35,009 A→G 6.6% intergenic (+27/+31) czrSA → / ← PP_0031 metal‑sensor histidine kinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294735,0090AG 6.6% 219.2 / 4.4 76intergenic (+27/+31)czrSA/PP_0031metal‑sensor histidine kinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base G (3/2);  total (40/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 = 9.74e-01

GATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTTC  >  NC_002947/34865‑35157
                                                                                                                                                |                                                                                                                                                     
gATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATc                                                                                                                                                    >  2:269574/1‑148 (MQ=255)
gATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGCTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTATCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGTCCCGTCaaac                                                                                                                                                    >  2:273604/1‑146 (MQ=255)
  tCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGAc                                                                                                                                                 >  2:34102/1‑149 (MQ=255)
    aTGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACgg                                                                                                                                               <  3:260188/149‑1 (MQ=255)
    aTGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCCGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACgg                                                                                                                                               <  1:273604/149‑1 (MQ=255)
      gAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGc                                                                                                                                             >  2:349114/1‑149 (MQ=255)
       aGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGc                                                                                                                                             <  4:54478/148‑1 (MQ=255)
         tCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCtt                                                                                                                                          >  3:170995/1‑149 (MQ=255)
                cAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGc                                                                                                                                             <  5:291144/139‑1 (MQ=255)
                cAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGc                                                                                                                                             >  6:291144/1‑139 (MQ=255)
                cAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTtttgtttg                                                                                                                                   >  8:3290/1‑149 (MQ=255)
                 aGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGATACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGc                                                                                                                                  <  4:201596/149‑1 (MQ=255)
                 aGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTtttgttt                                                                                                                                    >  3:231497/1‑147 (MQ=255)
                 aGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGc                                                                                                                                  >  7:66842/1‑149 (MQ=255)
                     aTCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTc                                                                                                                               >  5:331026/1‑148 (MQ=255)
                     aTCGGCGAGCGTGGAGGTGACCGAAGAAAACTTAACCTGGTTCCACCTTCACTTTACCGAACCCCAGCAAATGAAGCCGGCAATGACCCCCGCTTAAATTGCGCAGACACGAATAGGCCCGTCAATCGACGGGCCATTAGTTTCCCTc                                                                                                                               >  1:47650/1‑148 (MQ=255)
                         gCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCttt                                                                                                                                         >  8:29704/1‑134 (MQ=255)
                         gCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCttt                                                                                                                                         <  7:29704/134‑1 (MQ=255)
                           gAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCAt                                                                                                                        >  6:154355/1‑149 (MQ=255)
                             gCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCa                                                                                                                      <  5:208468/149‑1 (MQ=255)
                                tGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTg                                                                                                                   >  7:111042/1‑149 (MQ=255)
                                     gTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGaa                                                                                                              <  6:309945/149‑1 (MQ=255)
                                      tGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAg                                                                                                             <  5:311332/149‑1 (MQ=255)
                                       gACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGc                                                                                                                                  >  6:56889/1‑127 (MQ=255)
                                       gACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGc                                                                                                                                  <  5:56889/127‑1 (MQ=255)
                                       gACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAg                                                                                                             >  2:317600/1‑148 (MQ=255)
                                         ccGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGAc                                                                                                          <  6:226767/149‑1 (MQ=255)
                                          cGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTg                                                                                                                   <  2:301899/139‑1 (MQ=255)
                                                 aaCTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGt                                                                                                  >  3:54782/1‑149 (MQ=255)
                                                  aCTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTc                                                                                                 >  4:348207/1‑149 (MQ=255)
                                                    ttCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCa                                                                                               >  7:223476/1‑149 (MQ=255)
                                                      cACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCAt                                                                                              <  8:223476/148‑1 (MQ=255)
                                                        ccTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCAt                                                                                              <  1:331163/146‑1 (MQ=255)
                                                        ccTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCAt                                                                                              >  2:331163/1‑146 (MQ=255)
                                                         cTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGc                                                                                          >  3:134186/1‑149 (MQ=255)
                                                            gTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGAc                                                                                       >  8:95790/1‑149 (MQ=255)
                                                                 aCCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACcac                                                                                                       >  4:295134/1‑128 (MQ=255)
                                                                 aCCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACcac                                                                                                       <  3:295134/128‑1 (MQ=255)
                                                                   cTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGaa                                                                                >  6:172277/1‑149 (MQ=255)
                                                                       aCTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCAc                                                                            <  1:269574/149‑1 (MQ=255)
                                                                        cTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCAcc                                                                           <  8:66842/149‑1 (MQ=255)
                                                                          ttCCGGAAGCCCAGCAAATTCAGCCGGCAATGACGGCAGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGa                                                                                        <  6:20745/134‑1 (MQ=255)
                                                                          ttCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGa                                                                                        >  5:20745/1‑134 (MQ=255)
                                                                           tCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCg                                                                          <  5:154355/147‑1 (MQ=255)
                                                                            ccGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGAc                                                                        <  8:111042/148‑1 (MQ=255)
                                                                              ggAAGCCCAGCAAATGAAGCCGGCAATGACAGCCGATTAAATAGCGCACACACGATTAGGCCCGTCAATCGACGGGCCTTTTCATTCCCTCAATGCATCAGTGCCGAACTACCACAGGGTCCATGGCCGACCGCAGAATCACCCACta                                                                      >  8:337180/1‑147 (MQ=255)
                                                                               gAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTtttg                                                                                                                                       >  8:123810/1‑82 (MQ=255)
                                                                               gAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTtttg                                                                                                                                       <  7:123810/82‑1 (MQ=255)
                                                                                 aGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACAACAGGGTCCATGGGCGACCGCAGAATCACCGACTaaat                                                                   >  7:285710/1‑148 (MQ=255)
                                                                                   cccAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTaataat                                                                <  6:139019/149‑1 (MQ=255)
                                                                                    ccAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtc                                                               >  4:210436/1‑149 (MQ=255)
                                                                                     cAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtct                                                              <  5:172277/149‑1 (MQ=255)
                                                                                     cAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtct                                                              >  5:34782/1‑149 (MQ=255)
                                                                                      aGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtctt                                                             <  1:34102/149‑1 (MQ=255)
                                                                                               aGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGACCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGt                                                    >  7:209519/1‑148 (MQ=255)
                                                                                                 ccGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAg                                                 >  3:351521/1‑149 (MQ=255)
                                                                                                 ccGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAg                                                 >  4:90893/1‑149 (MQ=255)
                                                                                                    gCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTg                                                                                                                   <  5:123431/81‑1 (MQ=255)
                                                                                                    gCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTg                                                                                                                   >  6:123431/1‑81 (MQ=255)
                                                                                                     cAATGACCGCATCTTATATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGAtt                                             <  8:285710/149‑1 (MQ=255)
                                                                                                      aaTGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTg                                            <  7:337180/149‑1 (MQ=255)
                                                                                                              ccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCGTCAATGCATCAGTGCCGAAGGACGACAGGGTCCATGGGCGACCGCAGAATCTCCGACTTAATTATCTTCTGTTCGTCAGGATTGGCGttttt                                    <  2:47650/149‑1 (MQ=255)
                                                                                                              ccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGttttt                                    <  6:331026/149‑1 (MQ=255)
                                                                                                               cgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTc                                   <  7:3290/149‑1 (MQ=255)
                                                                                                                gcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTc                                   >  2:53245/1‑148 (MQ=255)
                                                                                                                          cgcAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTga                        <  8:55063/149‑1 (MQ=255)
                                                                                                                           gcAGTCACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTga                        <  1:53245/149‑1 (MQ=255)
                                                                                                                             aGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTgag                       >  7:339116/1‑147 (MQ=255)
                                                                                                                                cacGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTa                  <  7:95790/149‑1 (MQ=255)
                                                                                                                                aaacaaaaGGCCCGTCGATTGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTa                  >  6:159682/8‑149 (MQ=255)
                                                                                                                                   gATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACt                >  7:143147/1‑148 (MQ=255)
                                                                                                                                   gATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTg               <  1:349114/149‑1 (MQ=255)
                                                                                                                                   caaaaGGCCCGTCGATTGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTg               >  8:183564/5‑149 (MQ=255)
                                                                                                                                    aTTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGa              <  1:317600/149‑1 (MQ=255)
                                                                                                                                         gCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGaa                     <  3:189557/137‑1 (MQ=255)
                                                                                                                                         gCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGaa                     >  4:189557/1‑137 (MQ=255)
                                                                                                                                           ccGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTc        >  8:305439/1‑148 (MQ=255)
                                                                                                                                                aaTCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGgtgt                          >  5:7947/1‑125 (MQ=255)
                                                                                                                                                aaTCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGgtgt                          <  6:7947/125‑1 (MQ=255)
                                                                                                                                                aaTCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTtc  >  5:82055/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTTC  >  NC_002947/34865‑35157

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: