Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,896,803 G→A 6.4% R41Q (CGG→CAG)  PP_2551 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,896,8030GA 6.4% 231.0 / 4.4 78R41Q (CGG→CAG) PP_2551LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base A (2/3);  total (36/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 = 7.77e-01

CCTTGAATTCCACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCG  >  NC_002947/2896655‑2896946
                                                                                                                                                    |                                                                                                                                               
ccTTGAATTCCACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCg                                                                                                                                                 >  2:58451/1‑149 (MQ=255)
 cTTGAATTCCACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCgg                                                                                                                                                <  1:176592/149‑1 (MQ=255)
          cACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGc                                                                                                                                        <  1:58451/148‑1 (MQ=255)
          cACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCt                                                                                                                                       >  3:246403/1‑149 (MQ=255)
           aCCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTc                                                                                                                                      >  3:274311/1‑149 (MQ=255)
             ggCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCTCCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAAGCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCAc                                                                                                                                    <  7:243986/148‑1 (MQ=255)
              gCAGTGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCAGCTGCCGGCTCagg                                                                                                                                   <  8:290810/149‑3 (MQ=255)
              gCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCAcc                                                                                                                                   <  8:253865/149‑1 (MQ=255)
                 gggAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGc                                                                                                                                <  7:290831/149‑1 (MQ=255)
                  ggAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCg                                                                                                                               >  6:24428/1‑149 (MQ=255)
                  ggAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCg                                                                                                                               <  7:40680/149‑1 (MQ=255)
                   gAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCTCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGt                                                                                                                              <  5:43611/149‑1 (MQ=255)
                   gAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGt                                                                                                                              >  6:348108/1‑149 (MQ=255)
                   gAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCAGCTGCCGACTCag                                                                                                                                    >  4:82341/1‑142 (MQ=255)
                   gAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCAGCTGCCGACTCag                                                                                                                                    <  3:82341/143‑2 (MQ=255)
                      ccTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTg                                                                                                                             >  6:209931/1‑147 (MQ=255)
                      ccTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGaa                                                                                                                           <  6:117443/149‑1 (MQ=255)
                          cATGAACCCCACCTTCGCCTCCCTCAGCCTCGGCCACCTGCGCACCCTCGACCACGTGTTGAAGCTTAAAAACCTCAGCCTCGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTTATCCGGCAGCTGGCTCACCTGCGTGAAGCCt                                                                                                                       <  3:120084/149‑1 (MQ=255)
                           aTGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCtt                                                                                                                      <  2:320389/149‑1 (MQ=255)
                             gAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTc                                                                                                                     >  7:294665/1‑148 (MQ=255)
                               aCCCCACCTTCGCCTCCCTCAGCCTGGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCATCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTTAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCga                                                                                                                   <  6:126749/148‑1 (MQ=255)
                               aCCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCgat                                                                                                                  >  5:27063/1‑149 (MQ=255)
                               aCCCCACCTTCGCCTACCTCAGCCTGGCCCACCTGCCCAACCTCGACCAACTGTTGCATCTTAAAAACCTCAGACACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCga                                                                                                                   <  6:132595/148‑1 (MQ=255)
                                ccccACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCgat                                                                                                                  >  7:137413/1‑148 (MQ=255)
                                  ccACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCgatgat                                                                                                               <  4:330852/149‑1 (MQ=255)
                                      cTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGAGGATCCGt                                                                                                           >  6:163225/1‑149 (MQ=255)
                                       ttCGCCTCCCTCAGCCTAGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCGGCGTGAAGCCTTCGATGATCCGtt                                                                                                          >  3:54033/1‑149 (MQ=255)
                                          gCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCt                                                                                                       >  3:368497/1‑149 (MQ=255)
                                          gCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCAATCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCt                                                                                                       >  3:368492/1‑149 (MQ=255)
                                               cctcAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTgc                                                                                                  >  3:365126/1‑149 (MQ=255)
                                                ctcAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAg                                                                                                                          >  2:212334/1‑124 (MQ=255)
                                                ctcAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAg                                                                                                                          <  1:212334/124‑1 (MQ=255)
                                                ctcAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTgcg                                                                                                 <  5:209931/149‑1 (MQ=255)
                                                    gCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCaa                                                                                             <  8:68345/149‑1 (MQ=255)
                                                     ccTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAg                                                                                            >  3:176446/1‑149 (MQ=255)
                                                       tCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAggg                                                                                          >  6:160236/1‑149 (MQ=255)
                                                          cccACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGt                                                                                       <  4:274311/149‑1 (MQ=255)
                                                           ccgcgcgcgcCCCCGCGACCACCTGGTGCCGCTTACACCCCTCAGCCCCGCGGCCGAGCGCCTGGGGGGCAGCCAATCGGCCCGGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTg                                                                                      <  3:223133/143‑1 (MQ=255)
                                                           ccACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTg                                                                                      <  8:294665/149‑1 (MQ=255)
                                                             acttGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGaa                                                                                                                           <  3:87889/107‑1 (MQ=255)
                                                             acttGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGaa                                                                                                                           >  4:87889/4‑110 (MQ=255)
                                                                   gcACCCTCGACCCCCTGTTGCAGCTTAAAACCCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGCGCCGGCAGCTGGCTCACCTGCGGGAAGCAGTCGAGGATCCTGTGCGGGTACGCCAAGGGCGTg                                                                                      >  4:223133/1‑141 (MQ=255)
                                                                     aCCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCt                                                                            >  7:135535/1‑149 (MQ=255)
                                                                             ccaccTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAAc                                                                    <  8:137413/149‑1 (MQ=255)
                                                                                ccTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTg                                                                                      >  2:260295/1‑128 (MQ=255)
                                                                                ccTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTg                                                                                      <  1:260295/128‑1 (MQ=255)
                                                                                ccTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCAGCGTGAAGACTTCGATGATGCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATg                                                                 >  4:310119/1‑149 (MQ=255)
                                                                                 cTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAATCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGc                                                                <  6:27063/149‑1 (MQ=255)
                                                                                      gCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAAc                                                                    <  3:103597/140‑1 (MQ=255)
                                                                                      gCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAAc                                                                    >  4:103597/1‑140 (MQ=255)
                                                                                      gCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAg                                                           >  8:44511/1‑149 (MQ=255)
                                                                                          cTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCt                                                       <  8:135535/149‑1 (MQ=255)
                                                                                             aaaaaCCTCAGCCTCGCTGCCGAGCTAATGGGTTTCAGCCAATCGGCCCTGAGCCGGAAGCTGGATCACCTGCGTGAAGCATTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTTGTTATGTGCTTAGCGAACATGCCGAAGCCCTgg                                                     <  4:368492/148‑1 (MQ=255)
                                                                                             aaaaaCCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTgg                                                     <  4:368497/148‑1 (MQ=255)
                                                                                                   cTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCAGGTCGAGcc                                              >  5:364004/1‑149 (MQ=255)
                                                                                                      aGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGc                                                                                                                                <  3:113391/64‑1 (MQ=255)
                                                                                                      aGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGc                                                                                                                                >  4:113391/1‑64 (MQ=255)
                                                                                                       gCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTg                                          <  5:348108/149‑1 (MQ=255)
                                                                                                        ccACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGcc                                                               >  4:5774/1‑127 (MQ=255)
                                                                                                        ccACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCGGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTgc                                         <  5:24428/149‑1 (MQ=255)
                                                                                                             gggCCGAGCGCCTGTGTGTCAGCCATTCGGCCCTGAGACGGAAGCTGGCTCACCTGAGTGAAGCCTTCGCTGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGTTCGAGCCTTTGCGCCAg                                    <  4:54033/148‑1 (MQ=255)
                                                                                                                     cgcCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCga                            <  4:176446/149‑1 (MQ=255)
                                                                                                                      tccTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCgaa                           <  4:365126/148‑1 (MQ=255)
                                                                                                                        cTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGc                                                                        >  6:328636/1‑102 (MQ=255)
                                                                                                                        cTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGc                                                                        <  5:328636/102‑1 (MQ=255)
                                                                                                                          gggTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCg                             <  3:270680/143‑1 (MQ=255)
                                                                                                                          gggTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGACGAAGCCCTGGTCGAGCGTTTGAGCCAGGTGCTCgaatatc                       >  4:270680/1‑145 (MQ=255)
                                                                                                                           ggTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGc                                                                >  7:204321/1‑107 (MQ=255)
                                                                                                                           ggTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGc                                                                <  8:204321/107‑1 (MQ=255)
                                                                                                                           ggTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCAGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAAc                       <  6:364004/148‑1 (MQ=255)
                                                                                                                               tCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCAc                  >  2:378849/1‑149 (MQ=255)
                                                                                                                                  gCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAg                                                                                            <  6:7363/72‑1 (MQ=255)
                                                                                                                                  gCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAg                                                                                            >  5:7363/1‑72 (MQ=255)
                                                                                                                                  gCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGt                                                    >  6:91887/1‑112 (MQ=255)
                                                                                                                                  gCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGt                                                    <  5:91887/112‑1 (MQ=255)
                                                                                                                                           aggtGAGCCAGCTGCCGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTg                                          <  5:284920/110‑1 (MQ=255)
                                                                                                                                           aggtGAGACAGCTGCCGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTg                                          >  6:284920/4‑113 (MQ=255)
                                                                                                                                              tGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGACCTGCACGCGCTGCGCCAGccc   >  4:217548/1‑149 (MQ=255)
                                                                                                                                                  ccGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCg  >  8:194525/1‑146 (MQ=255)
                                                                                                                                                  ccGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCg  <  7:194525/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCTTGAATTCCACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCG  >  NC_002947/2896655‑2896946

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: