Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 188,618 G→A 5.5% G243R (GGG→AGG)  PP_0165 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947188,6180GA 5.5% 267.0 / 3.4 90G243R (GGG→AGG) PP_0165GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (42/43);  new base A (2/3);  total (44/46)
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.97e-02

CTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTT  >  NC_002947/188474‑188766
                                                                                                                                                |                                                                                                                                                            
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                                                                       >  5:354091/1‑107 (MQ=255)
                                                                                                                     gcgGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACct                                                                       <  6:354091/107‑1 (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)
                                                                                                                                            gACCGGGCTGCC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGt                                                                  <  3:29212/89‑1 (MQ=255)
                                                                                                                                            gACCGGGCTGCC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGt                                                                  >  4:29212/1‑89 (MQ=255)
                                                                                                                                             aCCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGAc       <  2:329651/147‑1 (MQ=255)
                                                                                                                                             aCCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGAc       <  4:216051/147‑1 (MQ=255)
                                                                                                                                             aCCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGAc       <  5:93923/147‑1 (MQ=255)
                                                                                                                                             aCCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCa     <  8:53858/149‑1 (MQ=255)
                                                                                                                                              ccGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAg    <  3:162686/149‑1 (MQ=255)
                                                                                                                                              ccGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAg    <  5:333375/149‑1 (MQ=255)
                                                                                                                                                gggCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAgtt  <  8:57248/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                            
CTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAACCGC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTT  >  NC_002947/188474‑188766

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: