Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 188,613 T→C 5.9% L241P (CTG→CCG)  PP_0165 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947188,6130TC 5.9% 257.8 / 4.5 85L241P (CTG→CCG) PP_0165GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (44/36);  new base C (2/3);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

CTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGA  >  NC_002947/188468‑188760
                                                                                                                                                 |                                                                                                                                                           
cTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGAc                                                                                                                                                          >  6:173415/1‑149 (MQ=255)
cTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGAc                                                                                                                                                          >  1:291747/1‑149 (MQ=255)
 tCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGAcc                                                                                                                                                         >  2:1504/1‑149 (MQ=255)
      cTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCt                                                                                                                                                    >  1:149929/1‑149 (MQ=255)
      cTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCt                                                                                                                                                    <  3:13887/149‑1 (MQ=255)
       tGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAAGTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTg                                                                                                                                                   >  3:245143/1‑149 (MQ=255)
         ccAGAGTTGCCACGCACGCCGGAGCTACGCCTCTTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGt                                                                                                                                                  <  2:253286/148‑1 (MQ=255)
         ccAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGt                                                                                                                                                  <  6:15211/148‑1 (MQ=255)
         ccAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTc                                                                                                                                                 <  3:71167/149‑1 (MQ=255)
          cAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGTAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcc                                                                                                                                                <  2:329189/149‑2 (MQ=255)
          cAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcc                                                                                                                                                <  3:202032/149‑2 (MQ=255)
             agTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAc                                                                                                                                           <  2:291747/149‑1 (MQ=255)
                 gCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAccgc‑‑c                                                                                                                                   >  5:379513/1‑149 (MQ=255)
                      gcacgcCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTAt                                                                                                                              >  2:147178/1‑149 (MQ=255)
                      gcacgcCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTCAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTAt                                                                                                                              >  5:370237/1‑149 (MQ=255)
                            cGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGaa                                                                                                                      <  6:370237/149‑1 (MQ=255)
                                   aCGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGc                                                                                                               >  5:296871/1‑149 (MQ=255)
                                    cgccgcGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAg                                                                                                                <  7:6274/147‑1 (MQ=255)
                                      ccgcGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTg                                                                                                             >  8:74606/1‑148 (MQ=255)
                                       cgcgTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGaa                                                                                                           >  7:143767/1‑149 (MQ=255)
                                       cgcgTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGaa                                                                                                           >  1:329651/1‑149 (MQ=255)
                                            aGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCGAGTT‑‑TGAAATGCAGCGGAACAATc                                                                                                      >  8:162351/1‑149 (MQ=255)
                                             gTGCAAGCTATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCTGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATc                                                                                                      <  7:68337/148‑1 (MQ=255)
                                               gCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGaa                                                                                                           <  2:321717/141‑1 (MQ=255)
                                               gCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGaa                                                                                                           >  1:321717/1‑141 (MQ=255)
                                                cAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCggg                                                                                                   >  3:111699/1‑148 (MQ=255)
                                                 aaccGATGAACAAGATGGTCGAGAAACTCAAAGCCCTGTTAACCGAGCAGGCAGACCGCAGTGAGAGGCTGCGTGCGGAGTCCTACAAGGACAGCATGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGt                                                                                                  >  2:339483/4‑148 (MQ=255)
                                                 aaGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGt                                                                                                  >  4:70720/1‑148 (MQ=255)
                                                 aaGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGt                                                                                                  >  1:330800/1‑148 (MQ=255)
                                                  aGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTaa                                                                                                <  6:109416/149‑1 (MQ=255)
                                                  aGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTaa                                                                                                <  6:296871/149‑1 (MQ=255)
                                                       tGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAAc                                                                                           <  5:173415/149‑1 (MQ=255)
                                                              gATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GCACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAg                                                                                    >  4:86012/1‑149 (MQ=255)
                                                              gATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTA‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTgg                                                                                       >  7:163818/1‑146 (MQ=255)
                                                                     gagaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgcgc                                                                              <  8:163818/148‑1 (MQ=255)
                                                                      agaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGaa                                                                                                                      <  1:373611/107‑1 (MQ=255)
                                                                      agaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGaa                                                                                                                      >  2:373611/1‑107 (MQ=255)
                                                                      agaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgcgcgc                                                                            <  7:74606/149‑1 (MQ=255)
                                                                       gaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAc                                                                                  <  7:148356/142‑1 (MQ=255)
                                                                       gaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAc                                                                                  >  8:148356/1‑142 (MQ=255)
                                                                        aaaCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACgcgcgcgc                                                                          >  5:122159/1‑149 (MQ=255)
                                                                        aaaCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACgcgcgcgc                                                                          >  1:241263/1‑149 (MQ=255)
                                                                        aaaCTCAAAGCCCTGTTCACCGAGCAGGCCGAGAGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACgcgcgcgc                                                                          <  7:97414/149‑1 (MQ=255)
                                                                          aCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCg                                                                        <  8:124377/149‑1 (MQ=255)
                                                                              aaaGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGAGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTa                                                                    >  6:333375/1‑149 (MQ=255)
                                                                                  cccTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACctg                                                                <  6:379513/149‑1 (MQ=255)
                                                                                  cccTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGACCTACCAGGACAGCCTGACAGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAACACGCGCGCGCCGGTTACctg                                                                >  6:93923/1‑149 (MQ=255)
                                                                                   ccTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACctgc                                                               <  8:143767/149‑1 (MQ=255)
                                                                                       ttCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACctgctg                                                             <  4:111699/147‑1 (MQ=255)
                                                                                        tCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTtg                                                          >  1:10289/1‑149 (MQ=255)
                                                                                          aCCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTtgct                                                        >  4:162686/1‑149 (MQ=255)
                                                                                           ccGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACgcgcgcgc                                                                          <  2:196352/130‑1 (MQ=255)
                                                                                           ccGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTtgctg                                                       >  4:261894/1‑149 (MQ=255)
                                                                                            cGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTCTCTGGGACCGC‑‑CGCTATTT‑‑TGAAAAGCAGCTGAACAATCGCGTAAGCAACCTGGAAGACGCGCGCGCCGCTTACCTGCTGTtgct                                                        >  1:118390/1‑149 (MQ=255)
                                                                                               gCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCggg                                                   >  3:216051/1‑149 (MQ=255)
                                                                                                     cgtgcgcAGTGACCGGCTGCCTGCGCACTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑TAACCGC‑‑CGCTATTT‑‑TGAAATTCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAgg                                             <  2:149929/146‑1 (MQ=255)
                                                                                                      gAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAggg                                            <  7:162915/149‑1 (MQ=255)
                                                                                                         cgcCGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTCTTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGGTACCTGCTGTTGCTGCGGGGTCAGGGCct                                         <  3:86012/149‑1 (MQ=255)
                                                                                                         cgcAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCct                                         <  1:1504/149‑1 (MQ=255)
                                                                                                          gcAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGGTGTC‑‑GAACTCAGCATCCCTTTTTACTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAg                                              <  1:311596/149‑1 (MQ=255)
                                                                                                                aGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAACACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCtggctggc                                  <  5:285080/149‑1 (MQ=255)
                                                                                                                  cGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTc                                                <  7:336837/133‑1 (MQ=255)
                                                                                                                  cGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTc                                                >  8:336837/1‑133 (MQ=255)
                                                                                                                  cGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCt                                <  1:147178/149‑1 (MQ=255)
                                                                                                                          caggctgTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgccc                        >  8:201768/7‑149 (MQ=255)
                                                                                                                           gcgGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACct                                                                 <  6:354091/107‑1 (MQ=255)
                                                                                                                           gcgGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACct                                                                 >  5:354091/1‑107 (MQ=255)
                                                                                                                                gTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTg                  >  1:116080/1‑149 (MQ=255)
                                                                                                                                gTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGGTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTg                  <  7:162351/149‑1 (MQ=255)
                                                                                                                                 tCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgccc                        >  7:192516/1‑142 (MQ=255)
                                                                                                                                 tCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgccc                        <  8:192516/142‑1 (MQ=255)
                                                                                                                                    tACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcgg              >  7:93955/1‑149 (MQ=255)
                                                                                                                                     aCCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc             >  1:234535/1‑149 (MQ=255)
                                                                                                                                        tcgACAGCCCGGTCAGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTAGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAg          >  3:116473/3‑149 (MQ=255)
                                                                                                                                         ggACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGc         >  1:376853/1‑149 (MQ=255)
                                                                                                                                         cgACAGCCCGGTCAGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGc         >  2:359315/2‑149 (MQ=255)
                                                                                                                                          gACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCg        >  7:57248/1‑149 (MQ=255)
                                                                                                                                             aGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTAc     <  7:201768/149‑1 (MQ=255)
                                                                                                                                             aGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTAc     <  2:10289/149‑1 (MQ=255)
                                                                                                                                               ccTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCTACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCg   >  2:148406/1‑149 (MQ=255)
                                                                                                                                               ccTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTCCCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgc                          >  8:90498/1‑126 (MQ=255)
                                                                                                                                               ccTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgc                          <  7:90498/126‑1 (MQ=255)
                                                                                                                                               ccTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTCGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTAcc    >  2:318254/1‑148 (MQ=255)
                                                                                                                                               ccTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATAT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCg   <  4:245143/149‑1 (MQ=255)
                                                                                                                                                cTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGa  >  4:198007/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                           
CTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGA  >  NC_002947/188468‑188760

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: