Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,437,250 G→C 7.9% S137T (AGC→ACC)  PP_3932 → diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,437,2500GC 7.9% 242.2 / 9.3 76S137T (AGC→ACC) PP_3932diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (36/34);  new base C (3/3);  total (39/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 = 7.86e-01

GTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGCGGGCCTTG  >  NC_002947/4437104‑4437386
                                                                                                                                                  |                                                                                                                                        
gtgGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTCGCGGCGAGCa                                                                                                                                        >  8:141347/1‑149 (MQ=255)
   gTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATgg                                                                                                                                     <  4:320763/149‑1 (MQ=255)
   gTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATg                                                                                                                                      <  3:110891/148‑1 (MQ=255)
   gTTAAACCTGTGCAATCTTCCCGTGGTCGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGAACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCAAGGTGGCGGCGAGCATgg                                                                                                                                     >  3:90868/1‑149 (MQ=255)
    ttAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATgg                                                                                                                                     >  1:355904/1‑148 (MQ=255)
          cTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCgg                                                                                                                               >  1:247458/1‑148 (MQ=255)
           tgtgCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGt                                                                                                                              >  6:196073/1‑148 (MQ=255)
             tgCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTa                                                                                                                             <  2:103634/147‑1 (MQ=255)
             tgCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAAt                                                                                                                           >  1:148632/1‑149 (MQ=255)
             tgCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAAt                                                                                                                           <  6:354541/149‑1 (MQ=255)
              gCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATg                                                                                                                          <  5:26891/149‑1 (MQ=255)
               cAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCAAGGTGGCGGCGAGCATGGCAGCGGTAATg                                                                                                                          <  1:203294/148‑1 (MQ=255)
               cAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCACCGCTGCCATGCTCGCCGCCACCATGGCAGCggcggc                                                                                                                           <  2:283391/147‑5 (MQ=38)
               cAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCACCGCTGCCATGCTCGCCGCCACCATGGCAGCggcggc                                                                                                                           >  1:283391/1‑143 (MQ=38)
                aaTCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGc                                                                                                                        >  5:95594/1‑149 (MQ=255)
                aaTCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGc                                                                                                                        <  7:59983/149‑1 (MQ=255)
                 aTCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGcc                                                                                                                       >  7:196549/1‑149 (MQ=255)
                      ggCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCAcc                                                                                                                  >  4:268963/1‑149 (MQ=255)
                       gCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCAccc                                                                                                                 <  3:268963/149‑1 (MQ=255)
                          tggtggTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATCGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTgg                                                                                                              >  6:215101/1‑149 (MQ=255)
                            gtggtCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGtt                                                                                                            <  5:196073/149‑1 (MQ=255)
                            gtggtCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCATGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGtt                                                                                                            >  2:166623/1‑149 (MQ=255)
                             tggtCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTc                                                                                                           <  2:148632/149‑1 (MQ=255)
                               gtCACCCTCTGGCAGCTGCTGTGGCGATTGCCACGCATGCACGGGCGCTTGCGCAGCCGGCACGGGGTACTGAGCGTGTTCGGGGGCTGCGGCGCCGCTGCGATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCg                                                                                                          <  4:90868/148‑1 (MQ=255)
                                  aCCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGaa                                                                                                        >  7:315169/1‑147 (MQ=255)
                                  aCCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAAc                                                                                                       <  2:355904/148‑1 (MQ=255)
                                    ccTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGt                                                                                                    <  8:237614/149‑1 (MQ=255)
                                          gCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGtt                                                                                                <  8:196549/147‑1 (MQ=255)
                                                 ctgTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCAc                                                                                       >  6:287173/1‑149 (MQ=255)
                                                      gcgcTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCAc                                                                                                                   <  8:14403/116‑1 (MQ=255)
                                                      gcgcTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGc                                                                                  >  3:293068/1‑149 (MQ=255)
                                                        gcTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCACCGCTGCCATGCTCGCCGCCACCATGGCa                                                                                                                                   >  5:223040/1‑98 (MQ=37)
                                                        gcTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCACCGCTGCCATGCTCGCCGCCACCATGGCa                                                                                                                                   <  6:223040/98‑1 (MQ=37)
                                                            gCCANGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCt                                                                            >  5:263810/1‑149 (MQ=255)
                                                            gCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGc                                                                              <  5:287173/147‑1 (MQ=255)
                                                             ccACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTg                                                                           <  6:263810/149‑1 (MQ=255)
                                                             ccACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTg                                                                           <  2:247458/149‑1 (MQ=255)
                                                                    tGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGc                                                                    >  1:65556/1‑149 (MQ=255)
                                                                    tGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGc                                                                    <  7:141347/149‑1 (MQ=255)
                                                                     gCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGc                                                                    >  5:185469/1‑148 (MQ=255)
                                                                       aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCgca                                                                 >  3:182797/1‑147 (MQ=255)
                                                                           gcgcATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGc                                                                    >  6:287976/1‑142 (MQ=255)
                                                                           gcgcATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGc                                                                    <  5:287976/142‑1 (MQ=255)
                                                                             gcATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTgg                                                                          >  8:234026/1‑134 (MQ=255)
                                                                             gcATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTgg                                                                          <  7:234026/134‑1 (MQ=255)
                                                                             gcATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTtc                                                           >  3:9995/1‑149 (MQ=255)
                                                                             gcATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCATGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTtc                                                           <  1:166623/149‑1 (MQ=255)
                                                                              cATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCAACCTGGCTAGCCTGGTTCAGCGAACAGTtct                                                          >  4:325997/1‑149 (MQ=255)
                                                                                   gcAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGGTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGAcc                                                     >  1:251405/1‑149 (MQ=255)
                                                                                    cAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCg                                                    >  3:296890/1‑149 (MQ=255)
                                                                                     aGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCgg                                                   >  6:316923/1‑149 (MQ=255)
                                                                                     aGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCgg                                                   >  8:16600/1‑149 (MQ=255)
                                                                                         cGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTtct                                                          <  4:69001/138‑1 (MQ=255)
                                                                                         cGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTtct                                                          >  3:69001/1‑138 (MQ=255)
                                                                                            aCGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCgtgctg                                            <  4:9995/149‑1 (MQ=255)
                                                                                            aCGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCgtgctg                                            >  4:156726/1‑149 (MQ=255)
                                                                                                 gTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGTCgtgctggtgct                                       <  4:293068/149‑1 (MQ=255)
                                                                                                    cTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGtgctgc                                     <  5:215101/148‑1 (MQ=255)
                                                                                                     tCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATTGCAGCGGTAATGGCCGCACCCAGCTTCGAACAGTCGTTGCGAGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCg                                   >  6:162220/1‑149 (MQ=255)
                                                                                                       aGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCgg                                  <  3:156726/148‑1 (MQ=255)
                                                                                                          gtgtTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCg                                   <  5:276511/144‑1 (MQ=255)
                                                                                                          gtgtTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCg                                   >  6:276511/1‑144 (MQ=255)
                                                                                                            gtTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGtgctgc                            >  8:162151/1‑149 (MQ=255)
                                                                                                            gtTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGtgctgc                            <  4:296890/149‑1 (MQ=255)
                                                                                                             ttCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGtgctgc                            <  5:316923/148‑1 (MQ=255)
                                                                                                             ttCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTgctgct                           >  8:184416/1‑149 (MQ=255)
                                                                                                                 ccgccTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTa                         >  4:282543/1‑147 (MQ=255)
                                                                                                                  cgccTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAg                      >  2:295469/1‑149 (MQ=255)
                                                                                                                  cgccTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAg                      <  2:251405/149‑1 (MQ=255)
                                                                                                                         gccgccgcTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTc               >  4:1810/1‑149 (MQ=255)
                                                                                                                          ccgccgcTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCg              <  7:16600/149‑1 (MQ=255)
                                                                                                                               gcTGCCATGCTCGCCGCCACCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGc                                                  >  5:283760/1‑108 (MQ=38)
                                                                                                                               gcTGCCATGCTCGCCGCCACCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGc                                                  <  6:283760/108‑1 (MQ=38)
                                                                                                                                   ccATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGCggg       <  1:295469/147‑1 (MQ=255)
                                                                                                                                     aTGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGCGGGCCtt   >  2:216131/1‑149 (MQ=255)
                                                                                                                                      tggtggCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGCGGGCCTTg  <  8:315169/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                        
GTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGCGGGCCTTG  >  NC_002947/4437104‑4437386

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