Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,368,197 G→T 5.5% intergenic (+26/‑157) PP_1193 → / → PP_1194 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,368,1970GT 5.5% 244.3 / 3.4 82intergenic (+26/‑157)PP_1193/PP_1194membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/39);  new base T (2/3);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

ATGTTGGTGGTGAACCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTGG  >  NC_002947/1368049‑1368337
                                                                                                                                                     |                                                                                                                                            
aTGTTGGTGCTGAACCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAg                                                                                                                                              <  3:331046/149‑1 (MQ=255)
     ggtggtGAACCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACg                                                                                                                                        >  2:282458/1‑149 (MQ=255)
     catgctgAAGCGGGTGACGATCATGGCAATCACGTGGGTGGCGAACCCGGTAATGAACACGGTGTTCATGACGGCGGCGAGGCAGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGATGAGGCCCATACAACAACCAg                                                                                                                                              >  4:102997/6‑144 (MQ=21)
      gtggtgAACCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGc                                                                                                                                             >  5:374508/1‑144 (MQ=255)
      gtggtgAACCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGc                                                                                                                                             <  6:374508/144‑1 (MQ=255)
         gtgAACCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGt                                                                                                                                     >  7:281997/1‑149 (MQ=255)
             aCCTGGTGACGATCACGGCAACCATGGGGGTGGCGAACCCGGTAATGACCAGGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTg                                                                                                                                 <  2:296668/149‑1 (MQ=255)
              ccTGGTGACGATCACGGCAACCATGTGGGGGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGt                                                                                                                                <  7:32125/149‑1 (MQ=255)
               cTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTa                                                                                                                               <  4:99658/149‑1 (MQ=255)
                tGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTAt                                                                                                                              >  4:227621/1‑149 (MQ=255)
                tGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACAACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCAGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTAt                                                                                                                              >  6:218275/1‑149 (MQ=255)
                 ggTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCGTCCGCTGGGTTTTGTATg                                                                                                                             <  7:48593/149‑1 (MQ=255)
                 ggTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTAt                                                                                                                              >  8:279743/1‑148 (MQ=255)
                  gTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGTCCATACAAAACCCAGCGGACGGGTTTTGTa                                                                                                                               <  6:15056/146‑1 (MQ=255)
                  gTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGTCCATACAAAACCCAGCGGACGGGTTTTGTa                                                                                                                               >  5:15056/1‑146 (MQ=255)
                  gTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTAt                                                                                                                              <  5:93262/147‑1 (MQ=255)
                      cGATCACGGAAACCATGTGTTTTGCGAACCATGTAATTACAACGGTTTTCATGCCGGCTGCGAGGCCGGGTATGACAAAAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCATCGGACGGGTTTTGTATGAGCAg                                                                                                                        <  8:120578/149‑1 (MQ=255)
                         tCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGtt                                                                                                                     <  2:266331/149‑1 (MQ=255)
                         tCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGtt                                                                                                                     >  8:295530/1‑149 (MQ=255)
                            cGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGc                                                                                                                          >  8:166471/1‑141 (MQ=255)
                            cGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGc                                                                                                                          <  7:166471/141‑1 (MQ=255)
                            cGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTg                                                                                                                  >  1:222072/1‑149 (MQ=255)
                            cGGCAACCATGTGGGTGGCGAACCCGGTAATGAACACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGCCGGGTTTTGTATGAGCAGGTTCTg                                                                                                                  >  4:252933/1‑149 (MQ=255)
                                aaCCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACgg                                                                                                                                       <  8:284993/124‑1 (MQ=255)
                                aaCCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACgg                                                                                                                                       >  7:284993/1‑124 (MQ=255)
                                aaCCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCt                                                                                                              <  3:227621/149‑1 (MQ=255)
                                             gCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTa                                                                                                  >  6:9602/1‑148 (MQ=255)
                                                   ccGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTc                                                                                           <  2:383610/149‑1 (MQ=255)
                                                    cGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCa                                                                                          >  1:212662/1‑149 (MQ=255)
                                                         aTGACGACGGTGTTCATGCCGTCGGCGAGGCCGGGGATGACAAGAGCAGTATCTGATTGGCGCTGACCCGTAGTGGCCCATCCAAAACCCAGCGGAGGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTc                                                                                      <  7:403054/148‑1 (MQ=255)
                                                         aTGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTc                                                                                      <  6:257331/148‑1 (MQ=255)
                                                          tGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGTCCTGTCTTGCCTTTATCACGTCATTTCCa                                                                                    <  5:218275/149‑1 (MQ=255)
                                                            aCCACGGTGTTCATGCCGGCGGCGATGCCGGGGATGAAAACAGCAGAAGCTGATTTGCTCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCaa                                                                                   <  7:55699/148‑1 (MQ=255)
                                                              cACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAAt                                                                                  >  1:75067/1‑147 (MQ=255)
                                                                     ttCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGa                                                                         <  5:9602/149‑1 (MQ=255)
                                                                     ttCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGa                                                                         <  3:200843/149‑1 (MQ=255)
                                                                     ttCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGATTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGa                                                                         >  3:297907/1‑149 (MQ=255)
                                                                       cATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGAc                                                                        <  7:295530/148‑1 (MQ=255)
                                                                        aTGCCGGCGGCTAGGCCGGGGATGACAACTGCATAAGCTGATTTGCGCTGACCCGATGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACa                                                                       <  6:96737/148‑1 (MQ=255)
                                                                        aTGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGcc                                                                                                      >  7:238748/1‑117 (MQ=255)
                                                                        aTGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGcc                                                                                                      <  8:238748/117‑1 (MQ=255)
                                                                         tGCCGGCGGCGAGGCCGGGGATGAGCAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAg                                                                      >  7:63441/1‑149 (MQ=255)
                                                                         tGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGc                                                                     >  6:49760/1‑149 (MQ=255)
                                                                           ccGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCtt                                                                   <  4:314340/149‑1 (MQ=255)
                                                                             ggcggcGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCtttt                                                                 <  3:385015/149‑1 (MQ=255)
                                                                                 gcgAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCt                                                             <  2:222072/149‑1 (MQ=255)
                                                                                 gcgAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGGATGAGCAGGATCTGGCCTGTCTTGCCTTTATCACGTCATTTCGAATGAGCGCTGACTGCTTTTTACt                                                             >  6:63125/1‑149 (MQ=255)
                                                                                              tGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCggg                                                <  2:212662/149‑1 (MQ=255)
                                                                                                  caacaGCAGAAGCTGATTGGCGCTGACCCGAAGTTGCCCATACAAAACCCAGCGGACGGGttt                                                                                                                                   <  4:237831/63‑1 (MQ=255)
                                                                                                  caacaGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGata                                             >  5:384486/1‑149 (MQ=255)
                                                                                                   aacaGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATg                                                     <  5:424172/140‑1 (MQ=255)
                                                                                                   aacaGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATg                                                     >  6:424172/1‑140 (MQ=255)
                                                                                                   aacaGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGatat                                            >  8:46562/1‑149 (MQ=255)
                                                                                                   aacaGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCCGACGGGTTTTGTATGAGCATGTTCTGGCCTGTCATGCCTTTATCACGTCATTTACAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGatat                                            >  6:348399/1‑149 (MQ=255)
                                                                                                      agcagAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGGCTGTCTTTCCTTTATCACGTCATTTCCTATGAGCGCTGACAGATTTTTCGTCGTCTATGCCGGGATATGAt                                         >  4:247686/1‑149 (MQ=255)
                                                                                                        cagAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCTGGATATGAtt                                        <  6:223649/148‑1 (MQ=255)
                                                                                                        cagAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGAtt                                        >  5:223649/1‑148 (MQ=255)
                                                                                                        cagAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTg                                       >  6:384479/1‑149 (MQ=255)
                                                                                                        cagAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTg                                       >  8:337968/1‑149 (MQ=255)
                                                                                                         agAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGc                                      >  3:279426/1‑149 (MQ=255)
                                                                                                          gAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCg                                     <  1:282458/149‑1 (MQ=255)
                                                                                                          gAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGATTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCg                                     <  4:297907/149‑1 (MQ=255)
                                                                                                          gAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACCGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATAACGTCATTTCCAATGATCCCTGACAGCTTTTTCATCGTCCATGCCGGGATATCATTgag                                     >  4:68900/1‑147 (MQ=255)
                                                                                                              cTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGttt                                 <  3:68900/149‑1 (MQ=255)
                                                                                                                 agtGGCGCTGTCCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGTCTTTATCACGTGATTTCCAATGAGAGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGTTTGCTGTTTg                                <  4:79544/145‑1 (MQ=255)
                                                                                                                 aTTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGTGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGttt                                 <  7:279743/146‑1 (MQ=255)
                                                                                                                 aTTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGc                               <  7:337968/148‑1 (MQ=255)
                                                                                                                  ttGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGt                             <  8:63441/149‑1 (MQ=255)
                                                                                                                  ttGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGt                             <  3:93368/149‑1 (MQ=255)
                                                                                                                     gcgcTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTgg                           >  6:240226/1‑148 (MQ=255)
                                                                                                                      cgcTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGa                         <  5:49760/149‑1 (MQ=255)
                                                                                                                       gcTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGAt                        >  5:249705/1‑149 (MQ=255)
                                                                                                                           aCCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTaa                     <  7:46562/148‑1 (MQ=255)
                                                                                                                             ccGAAGTGGCCCATACAAAACCCAGCGGACGCGTTTTGTATGAGCAGGTTCTGGCCTGTCTAGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCAGGGATATGATTGCGGTTTGCGTGGGATTAATTc                  >  8:411200/1‑149 (MQ=255)
                                                                                                                               gAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTccc                >  3:245737/1‑149 (MQ=255)
                                                                                                                                aaGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGAGATGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCcg               >  2:103729/1‑149 (MQ=255)
                                                                                                                                 agttccgcATAGAAAACCCAGCGTACTGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTTCCTTTATGACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCcgc              <  3:252933/142‑1 (MQ=255)
                                                                                                                                     gctcATACAAAACCCGTCCGCTGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCt          >  4:56506/4‑149 (MQ=255)
                                                                                                                                     gCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCt          >  7:56183/1‑149 (MQ=255)
                                                                                                                                          tACAAAACCCGTCCGCTGGGTTTTGTATGGGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGAt                        >  4:11743/1‑130 (MQ=39)
                                                                                                                                          tACAAAACCCGTCCGCTGGGTTTTGTATGGGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGAt                        <  3:11743/130‑1 (MQ=39)
                                                                                                                                             aaaaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTgg  >  7:205749/1‑149 (MQ=255)
                                                                                                                                             aaaaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTgg  <  6:384486/149‑1 (MQ=255)
                                                                                                                                             aaaaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTgg  >  4:367726/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
ATGTTGGTGGTGAACCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAATGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTGG  >  NC_002947/1368049‑1368337

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