Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 188,322 T→C 7.1% L144P (CTT→CCT) ‡ PP_0165 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947188,3220TC 7.1% 197.9 / 4.2 70L144P (CTT→CCT) ‡PP_0165GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (31/34);  new base C (2/3);  total (33/37)
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.24e-01

GAGCCAGACCCTGGTGGTGTGCCGCTTTGGTTCGTGCGCCTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGCTATTGCCCGCCGCGAATTCCT  >  NC_002947/188177‑188469
                                                                                                                                                 |                                                                                                                                                   
gAGCCAGACCCTGGTGGTGTGCCGCTTTGGTTCGTGCGCCTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTg                                                                                                                                                  <  8:186970/149‑1 (MQ=255)
   ccAGACCCTGGTGGTGTGCCGATTTGGTTCGTGCGCCTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCACGCCGCGCGCGATGAAGTGATCAGCCACCCCATGTTCGCCATTGCCAAGCTTTGGc                                                                                                                                                >  8:220471/1‑148 (MQ=255)
       aCCCTGGTGGTGTGCCGCTTTGGTTCGTGCGCCTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAgc                                                                                                                                           <  3:64209/149‑1 (MQ=255)
       aCCCTGGTGGTGTGCCGCTTTGGTTCGTGCGCCTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAAGCTTGGCaatgg                                                                                                                                           <  2:79731/149‑4 (MQ=255)
         ccTGGTGGTGTGCCGCTTTGGTTCGTGCGCCTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAgcgc                                                                                                                                         <  1:105598/149‑1 (MQ=255)
           tggtggTGTGCCGCTTTGGTTCGTGCGCCTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCt                                                                                                                                       >  8:242391/1‑149 (MQ=255)
              tggtgTGCCGCTTTTGTTCGTGCGCCTGATCGGCCTGGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGTCTGGCAGCATGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAATCGCCTTgg                                                                                                                                    <  2:112472/149‑1 (MQ=255)
                             gTTCGTGCGCCTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGggttggt                                                                                                                     >  8:251064/1‑149 (MQ=255)
                                  tGCGCCTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCtg                                                                                                                <  7:95735/149‑1 (MQ=255)
                                      ccTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCtgttgt                                                                                                            <  1:167124/149‑1 (MQ=255)
                                       cTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTtgtg                                                                                                           <  7:251064/149‑1 (MQ=255)
                                         gATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAAGCTTGGCaatg                                                                                                                                            <  4:330679/114‑3 (MQ=255)
                                         gATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAAGCTTGGCaatg                                                                                                                                            >  3:330679/1‑112 (MQ=255)
                                          aTCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGTCAAAGCGCCTTTGGCAGCCTGGGTTGGTTGCTGTTgtgcg                                                                                                         <  7:220471/148‑1 (MQ=255)
                                           tCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTgtgcggt                                                                                                       >  1:286287/1‑149 (MQ=255)
                                              gCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGTCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTgtgcggtgcg                                                                                                    >  4:264498/1‑149 (MQ=255)
                                              gCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTgg                                                                                                                           >  4:245948/1‑126 (MQ=255)
                                              gCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTgg                                                                                                                           <  3:245948/126‑1 (MQ=255)
                                                cTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGc                                                                                                  >  7:147773/1‑149 (MQ=255)
                                                 tCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTgtgcggtgc                                                                                                     <  7:24629/145‑1 (MQ=255)
                                                 tCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTgtgcggtgc                                                                                                     >  8:24629/1‑145 (MQ=255)
                                                 tCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGc                                                                                                  >  8:268223/1‑148 (MQ=255)
                                                        cgccgGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAAGCTTGGCaatgg                                                                                                                                           >  7:94432/1‑97 (MQ=255)
                                                        cgccgGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAAGCTTGGCaatgg                                                                                                                                           <  8:94432/100‑4 (MQ=255)
                                                         gccgGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGt                                                                                         >  5:290303/1‑149 (MQ=255)
                                                         gccgGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCACGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGc                                                                                                                                                >  2:86172/1‑94 (MQ=255)
                                                              cggcgATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCg                                                                                    <  7:242391/149‑1 (MQ=255)
                                                              cggcgATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCg                                                                                    <  1:142768/149‑1 (MQ=255)
                                                                  gATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGc                                                                                                  >  4:81883/1‑131 (MQ=255)
                                                                  gATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGc                                                                                                  <  3:81883/131‑1 (MQ=255)
                                                                     gCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCttt                                                                             >  3:75197/1‑149 (MQ=255)
                                                                       aaTCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCt                                                                                                                             <  4:52455/99‑1 (MQ=255)
                                                                       aaTCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCt                                                                                                                             >  3:52455/1‑99 (MQ=255)
                                                                        aTCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCt                                                                          >  2:36088/1‑149 (MQ=255)
                                                                        aTCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCt                                                                          <  6:194166/149‑1 (MQ=255)
                                                                        aTCGTCAGTCGGGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCt                                                                          <  8:147773/149‑1 (MQ=255)
                                                                           gTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGtt                                                                       >  8:102037/1‑149 (MQ=255)
                                                                           gTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGtt                                                                       >  1:192894/1‑149 (MQ=255)
                                                                              aGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGTCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTgc                                                                     <  3:264498/148‑1 (MQ=255)
                                                                              aGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTgcg                                                                    >  7:165279/1‑149 (MQ=255)
                                                                                  gTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGcgccgc                                                                >  5:266064/1‑149 (MQ=255)
                                                                                   tggctggcAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGcgccgc                                                                >  2:285730/1‑148 (MQ=255)
                                                                                       tggcAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGTGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGtt                                                           >  3:39494/1‑149 (MQ=255)
                                                                                         gcagcaGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGtt                                                           <  2:286287/147‑1 (MQ=255)
                                                                                         gcagcaGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGa                                                         <  5:311652/149‑1 (MQ=255)
                                                                                            gcagGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGc                                                      <  1:36088/149‑1 (MQ=255)
                                                                                             cagGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGcc                                                     >  5:279639/1‑149 (MQ=255)
                                                                                                 ccGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGTGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTg                                                 <  4:39494/149‑1 (MQ=255)
                                                                                                 ccGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTg                                                 <  1:285730/149‑1 (MQ=255)
                                                                                                               gTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTaga                                   <  4:75197/149‑1 (MQ=255)
                                                                                                                    cAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAAt                              >  8:34642/1‑149 (MQ=255)
                                                                                                                        cACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGcc                                                                                                 >  6:116259/1‑78 (MQ=255)
                                                                                                                        cACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGcc                                                                                                 <  5:116259/78‑1 (MQ=255)
                                                                                                                        cACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCa                          <  6:266064/149‑1 (MQ=255)
                                                                                                                         tcccGATGTTCGCCATTGCCATGCTTTGTCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTa                                     <  6:258027/136‑1 (MQ=255)
                                                                                                                         aCCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGTCAGCCTGGCTTGCTTGCTGTAGTGCGGTGCTGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTa                                     >  5:258027/1‑137 (MQ=255)
                                                                                                                           ccGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGc                       <  2:192894/149‑1 (MQ=255)
                                                                                                                               tGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGCTAtt                   <  6:290303/149‑1 (MQ=255)
                                                                                                                                  tCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCACCGCCAGTTGAGGCCGCTGGACTACATg                                        >  4:66203/1‑125 (MQ=255)
                                                                                                                                  tCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCACCGCCAGTTGAGGCCGCTGGACTACATg                                        <  3:66203/125‑1 (MQ=255)
                                                                                                                                    gCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGCTATTGc                 <  4:33242/146‑1 (MQ=255)
                                                                                                                                    gCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGCTATTGc                 >  3:33242/1‑146 (MQ=255)
                                                                                                                                           ccAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGCt                      <  3:251674/134‑1 (MQ=255)
                                                                                                                                           ccAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGCt                      >  4:251674/1‑134 (MQ=255)
                                                                                                                                            cAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGCTATTGCCCGCCGCGAAt      >  3:316579/1‑149 (MQ=255)
                                                                                                                                             aaGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATgg                                       <  2:111028/115‑1 (MQ=255)
                                                                                                                                             aaGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATgg                                       >  1:111028/1‑115 (MQ=255)
                                                                                                                                             aaGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGCTATTGCCCGCCGCGAAtt     >  4:138976/1‑149 (MQ=255)
                                                                                                                                               gCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGCTATTGCCCGCCGCGAATTcc   <  6:279639/149‑1 (MQ=255)
                                                                                                                                                cTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGCTATTGCCCGCCGCGAATTCCt  <  8:165279/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GAGCCAGACCCTGGTGGTGTGCCGCTTTGGTTCGTGCGCCTGATCGGCCTCGAGGCCGCCGGCGGCGATGCAATCGTCAGTCGTGGCTGGCAGCAGGCCGCGCGCGTTGAAGTGATCAGCCACCCGATGTTCGCCATTGCCAAGCTTTGGCAAAGCGCCTTGGGCAGCCTGGGTTGGTTGCTGTTGTGCGGTGCGGCCAGTGCCGTGCTCGGTGCTTTGCTGTTGCGCCGCCAGTTGAGGCCGCTGGACTACATGGTAGAGCAATCTCACGCTATTGCCCGCCGCGAATTCCT  >  NC_002947/188177‑188469

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: