Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,234,559 A→G 6.7% L169P (CTG→CCG)  PP_3711 ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,234,5590AG 6.7% 171.9 / 3.0 60L169P (CTG→CCG) PP_3711diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base G (2/2);  total (30/30)
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.10e-01

CGGGCACCCCGCACGAATTGCTGCTCCAGATAGTCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTGCAGCAGTTCGAGCTCCTGGCGGTGGCGTTCGGAAATGTCGCG  >  NC_002947/4234412‑4234707
                                                                                                                                                     |                                                                                                                                                    
cGGGCACCCCGCACGAATTGCTGCTCCAGATAGTCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCa                                                                                                                                                      <  2:238460/148‑1 (MQ=255)
    cACCCCGCACGAATTGCTGCTCCAGATAGTCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTcc                                                                                                                                                 <  6:25078/149‑1 (MQ=255)
     aCCCCGCACGAATTGCTGCTCCAGATAGTCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCt                                                                                                                                                >  4:331359/1‑149 (MQ=255)
       cccGCACGAATTGCTGCTCCAGATAGTCACCGGCTAGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTg                                                                                                                                               >  1:50118/1‑149 (MQ=255)
            aCGAATTGCTGCTCCAGATAGTCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGAt                                                                                                                                         <  3:342245/149‑1 (MQ=255)
             cGAATTGCTGCTCCAGATAGTCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATg                                                                                                                                        <  1:248362/149‑1 (MQ=255)
                       ctCCAGATAGTCACCGGC‑AGCCGGTGCCAGGCATTCCCTTGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGACGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTcgc                                                                                                                                <  6:119921/147‑1 (MQ=255)
                       ctCCAGATAGTCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCgg                                                                                                                              >  7:309004/1‑149 (MQ=255)
                       ctCCAGATAGTCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCgg                                                                                                                              >  6:162735/1‑149 (MQ=255)
                              ctggcACCGGC‑TGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGtt                                                                                                                       >  7:339588/5‑149 (MQ=255)
                                 tCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGt                                                                                                                        <  6:369975/145‑1 (MQ=255)
                                 tCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGt                                                                                                                        >  5:369975/1‑145 (MQ=255)
                                 tCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGc                                                                                                                    <  6:259474/149‑1 (MQ=255)
                                   aCCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTg                                                                                                                      <  1:134750/145‑1 (MQ=255)
                                   aCCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTg                                                                                                                      >  2:134750/1‑145 (MQ=255)
                                   aCCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCaa                                                                                                                  >  6:179537/1‑149 (MQ=255)
                                       gC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGAc                                                                                                              <  2:50118/149‑1 (MQ=255)
                                          aGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGACAAGTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGAcc                                                                                                             >  1:111398/1‑149 (MQ=255)
                                            ccGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCTCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGaa                                                                                                          <  8:37902/149‑1 (MQ=255)
                                            ccGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGaa                                                                                                          >  7:18066/1‑149 (MQ=255)
                                              ggTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATc                                                                                                        <  8:309004/149‑1 (MQ=255)
                                              ggTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATc                                                                                                        <  3:331359/149‑1 (MQ=255)
                                                tGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCga                                                                                                      <  2:295023/149‑1 (MQ=255)
                                                tGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCga                                                                                                      <  8:18066/149‑1 (MQ=255)
                                                       cATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCgg                                                                                                                              <  6:211043/118‑1 (MQ=255)
                                                       cATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCgg                                                                                                                              >  5:211043/1‑118 (MQ=255)
                                                            ccTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGa                                                                                            >  3:185765/1‑147 (MQ=255)
                                                                ggCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCtt                                                                                      <  4:210148/149‑1 (MQ=255)
                                                                ggCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCtt                                                                                      <  4:224802/149‑1 (MQ=255)
                                                                ggCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCtt                                                                                      <  1:226179/149‑1 (MQ=255)
                                                                     cTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCGGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGAcc                                                                                                             <  5:231519/121‑1 (MQ=255)
                                                                     cTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCGGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGAcc                                                                                                             >  6:231519/1‑121 (MQ=255)
                                                                      tGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCAt                                                                                >  1:332586/1‑149 (MQ=255)
                                                                        gCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCAttt                                                                              <  2:111398/149‑1 (MQ=255)
                                                                           gATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCGGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCt                                                                           >  2:213808/1‑149 (MQ=255)
                                                                                  ctgtgctgAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTg                                                                                     >  8:202574/1‑132 (MQ=255)
                                                                                  ctgtgctgAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTg                                                                                     <  7:202574/132‑1 (MQ=255)
                                                                                          aGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCg                                                             <  2:332586/148‑1 (MQ=255)
                                                                                            gCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGc                                                           >  5:76113/1‑148 (MQ=255)
                                                                                             cccGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGc                                                                                   <  4:221037/123‑1 (MQ=255)
                                                                                             cccGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGc                                                                                   >  3:221037/1‑123 (MQ=255)
                                                                                              ccgccgAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTg                                                                                     >  2:208701/1‑120 (MQ=255)
                                                                                              ccgccgAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAAAGACCGAATCGAGAAACTGCGAAAGCTTg                                                                                     <  1:208701/120‑1 (MQ=255)
                                                                                               cgccgAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGc                                                       >  6:219022/1‑149 (MQ=255)
                                                                                                     aCAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTg                                                 >  3:259553/1‑149 (MQ=255)
                                                                                                       aGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCGGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGc                                               <  1:213808/149‑1 (MQ=255)
                                                                                                       aGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGc                                                                                   <  6:27560/113‑1 (MQ=255)
                                                                                                       aGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGc                                                                                   >  5:27560/1‑113 (MQ=255)
                                                                                                        gAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTgg                                                                         >  7:315174/1‑122 (MQ=255)
                                                                                                        gAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTgg                                                                         <  8:315174/122‑1 (MQ=255)
                                                                                                             ttCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGAACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTg                                            >  1:87857/1‑147 (MQ=255)
                                                                                                                   ctgctgATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTGCAGCAGtt                                    >  1:285699/1‑149 (MQ=255)
                                                                                                                       tgATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTGCAGCAGTTCGAg                                >  1:227927/1‑149 (MQ=255)
                                                                                                                            ggCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTGCAGCAGTTCGAGCTCCt                           >  3:266522/1‑149 (MQ=255)
                                                                                                                                 aCAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTGCAGCAGTTCGAGCTCCTGGCgg                      <  2:285699/149‑1 (MQ=255)
                                                                                                                                         cgcTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTGCAGCAGTTCGAGCTCCTGGCGGTGGCGTTc              <  8:339588/149‑1 (MQ=255)
                                                                                                                                          gcTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCGGCAGCAGTTCGAGCTCCTGGCGGTGGCGTTCg             >  8:336182/1‑149 (MQ=255)
                                                                                                                                              gCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTGCAGCAGTTCGAGCTCCTGGCGGTGGCGTTCGGaa          <  4:185765/148‑1 (MQ=255)
                                                                                                                                              gCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCAGTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTGCAGCAGTTCGAGCTCCTGGCGGTGGCGGTCGGaat         >  8:155483/1‑148 (MQ=255)
                                                                                                                                                    cAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAAt                                                                                                         >  5:52695/1‑47 (MQ=255)
                                                                                                                                                    cAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAAt                                                                                                         <  6:52695/47‑1 (MQ=255)
                                                                                                                                                     aGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTGCAGCAGTTCGAGCTCCTGGCGGTGGCGTTCGGAAATGTcgcg  >  7:102952/1‑149 (MQ=255)
                                                                                                                                                     aGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTGCAGCAGTTCGAGCTCCTGGCGGTGGCGTTCGGAAATGTcgcg  >  6:368489/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGGCACCCCGCACGAATTGCTGCTCCAGATAGTCACCGGC‑AGCCGGTGCCAGGCATTCCCTGGGCAACTGGCCGATCATGCTGTGCTGAGGCCCGCCGAACAGAC‑GTTCGGCCTGCTGATTGGCCAACAGGATGCGCTGGCTTTCCAGGTCCTGGACGATGACGCTCGCGGGAATGTTGGCAATGACCGAATCGAGAAACTGCGAAAGCTTGGCCATTTCCTGGCTGTACTGCTCGGCCAGCTCTTTGGCCTGCAGCAGTTCGAGCTCCTGGCGGTGGCGTTCGGAAATGTCGCG  >  NC_002947/4234412‑4234707

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