Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,247,811 G→A 7.5% S596S (AGC→AGT) ‡ PP_1983 ← PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,247,8110GA 7.5% 193.2 / 4.0 67S596S (AGC→AGT) ‡PP_1983PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (27/34);  new base A (3/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

CGCGGGCGGTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCGGTGATGGTGATGGCGT  >  NC_002947/2247663‑2247956
                                                                                                                                                      |                                                                                                                                                  
cgcgGGCGGTCAGCGGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAgcgcg                                                                                                                                                    >  7:334019/1‑149 (MQ=255)
   gggcggTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCgcgctg                                                                                                                                                 >  3:236640/1‑149 (MQ=255)
        gTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGcgctgcgctg                                                                                                                                            <  3:315636/149‑1 (MQ=255)
         tCAGCAGCAGGCAATCCACTTCGTTGAGTGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATCTCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCgctgcgctgc                                                                                                                                           <  7:75756/149‑1 (MQ=255)
         tCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCgctgcgctgc                                                                                                                                           >  1:65522/1‑149 (MQ=255)
           agcagcagGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCtt                                                                                                                                         <  4:351288/149‑1 (MQ=255)
           agcagcagGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCAGCGCgctgctg                                                                                                                                         <  7:34054/149‑2 (MQ=255)
                   gCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCgg                                                                                                                                 <  2:254208/149‑1 (MQ=255)
                       tCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTgt                                                                                                                                       <  6:238973/139‑1 (MQ=255)
                       tCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTgt                                                                                                                                       >  5:238973/1‑139 (MQ=255)
                       tCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGcc                                                                                                                               >  4:342836/1‑147 (MQ=255)
                          aCTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGAt                                                                                                                          <  8:48169/149‑1 (MQ=255)
                           cTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATg                                                                                                                         >  5:224466/1‑149 (MQ=255)
                           cTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGAGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATg                                                                                                                         >  2:289758/1‑149 (MQ=255)
                              cGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑c                                                                                                                       <  1:163017/147‑1 (MQ=255)
                               gTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGa                                                                                                                    >  1:254867/1‑149 (MQ=255)
                                ttGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGAt                                                                                                                   <  5:187992/149‑1 (MQ=255)
                                   aGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGAtctc                                                                                                                <  4:352198/149‑1 (MQ=255)
                                   aGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGAtct                                                                                                                 >  3:387810/1‑148 (MQ=255)
                                    gggTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTTCACAGGAAAGCTCGCTCAGACAAGGCCAGCAGCCCCGTGCGGTGGTTGTAGTCGGCCCCGATG‑CCGAtct                                                                                                                 >  1:126975/1‑148 (MQ=255)
                                       tGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATa                                                                                                             >  7:278531/1‑148 (MQ=255)
                                            ggTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCgt                                                                                                       <  6:224466/149‑1 (MQ=255)
                                             gTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCAGCGCgctgctggccga                                                                                                                                    >  5:235433/1‑114 (MQ=255)
                                             gTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCAGCGCgctgctggccga                                                                                                                                    <  6:235433/120‑7 (MQ=255)
                                                 cACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGc                                                                                                  <  2:254867/149‑1 (MQ=255)
                                                  aCCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑AATCGGCCAGCAGCGCGATGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGc                                                                                                  >  1:64020/1‑149 (MQ=255)
                                                    cGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGa                                                                                               >  2:313584/1‑149 (MQ=255)
                                                     gAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGAt                                                                                              <  4:389437/149‑1 (MQ=255)
                                                      aaaCCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATc                                                                                             <  3:342836/149‑1 (MQ=255)
                                                         ccTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACt                                                                                          <  4:387810/149‑1 (MQ=255)
                                                           tCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACtt                                                                                         >  7:106205/1‑148 (MQ=255)
                                                             gACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTtc                                                                                        <  2:65522/147‑1 (MQ=255)
                                                                   gAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTgg                                                                                                   <  7:362764/130‑1 (MQ=255)
                                                                   gAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTgg                                                                                                   >  8:362764/1‑130 (MQ=255)
                                                                   gAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGc                                                                                >  2:336485/1‑149 (MQ=255)
                                                                      gATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGt                                                                             <  7:277635/149‑1 (MQ=255)
                                                                       atatCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTa                                                                            <  5:143953/149‑1 (MQ=255)
                                                                          tCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGcc                                                                         <  8:106205/149‑1 (MQ=255)
                                                                          tCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGcc                                                                         <  8:278531/149‑1 (MQ=255)
                                                                             cctcGCTTG‑TTGCGCCCGAGGATCTGGGGGTTGTGTACCCGGCCGGGCTCGCGCA‑ACTCGGCCAGCATGGCGGTGCGCTGCTTGTGGTTGGCCCAGATGCCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCgg                                                                       <  1:289758/146‑1 (MQ=255)
                                                                             ccACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGTCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGt                                                                      <  2:281282/149‑1 (MQ=255)
                                                                             ccACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGt                                                                      <  6:302722/149‑1 (MQ=255)
                                                                             ccACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGt                                                                      >  7:341689/1‑149 (MQ=255)
                                                                                   ttg‑ttgCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAggcg                                                                >  8:69199/1‑149 (MQ=255)
                                                                                       ttgCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACTCCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAggcgg                                                               >  7:180616/1‑147 (MQ=255)
                                                                                         gcgcCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGCTG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAg                                                            >  3:134111/1‑148 (MQ=255)
                                                                                          cgcCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGaa                                                          <  1:313584/149‑1 (MQ=255)
                                                                                           gcCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAg                                                         >  8:19339/1‑149 (MQ=255)
                                                                                             ccGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAggcggc                                                              >  8:368823/1‑142 (MQ=255)
                                                                                             ccGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAggcggc                                                              <  7:368823/142‑1 (MQ=255)
                                                                                                   aTCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGtt                                                 <  6:105411/149‑1 (MQ=255)
                                                                                                    tctcGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGtt                                                 >  8:116672/1‑148 (MQ=255)
                                                                                                      tcGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAc                                              >  4:334075/1‑149 (MQ=255)
                                                                                                           gtggtgTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCa                                         <  3:334075/149‑1 (MQ=255)
                                                                                                              gtgtACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAgg                                       <  7:19339/148‑1 (MQ=255)
                                                                                                               tgtACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTa                                     <  4:236640/149‑1 (MQ=255)
                                                                                                                  aCCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCTGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCg                                  >  8:76335/1‑149 (MQ=255)
                                                                                                                  aCCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCg                                  >  7:6145/1‑149 (MQ=255)
                                                                                                                         cgcgCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCtg                           >  3:152095/1‑149 (MQ=255)
                                                                                                                          gcgcTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCtgt                          <  8:334019/149‑1 (MQ=255)
                                                                                                                                gcgcA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGtt                    <  6:213828/149‑1 (MQ=255)
                                                                                                                                  gcA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCg                  >  4:169392/1‑149 (MQ=255)
                                                                                                                                    a‑aCTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggt                >  8:26171/1‑149 (MQ=255)
                                                                                                                                      aCTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtg               <  7:168078/149‑1 (MQ=255)
                                                                                                                                        tCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgat             >  2:245613/1‑149 (MQ=255)
                                                                                                                                        tCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGGTCggtgat             >  4:93904/1‑149 (MQ=255)
                                                                                                                                          ggCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatg            <  8:6145/148‑1 (MQ=255)
                                                                                                                                          ggCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgat             <  7:69199/147‑1 (MQ=255)
                                                                                                                                            caagcagcGCAGCGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggt          >  3:64055/3‑148 (MQ=255)
                                                                                                                                             cagcagCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtga        <  4:134111/149‑1 (MQ=255)
                                                                                                                                             aagcagcGCAGCGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGTGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtga        >  4:319387/2‑149 (MQ=255)
                                                                                                                                               gcagcGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtgat       <  1:336485/148‑1 (MQ=255)
                                                                                                                                                   cgcgcTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCGGTGATGGTGATGGCGt  <  7:76335/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CGCGGGCGGTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTG‑TTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCA‑ACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATG‑CCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCGGTGATGGTGATGGCGT  >  NC_002947/2247663‑2247956

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