Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,706,193 C→T 8.3% G100D (GGC→GAC)  PP_4162 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,706,1930CT 8.3% 205.7 / 5.0 71G100D (GGC→GAC) PP_4162membrane protein
Reads supporting (aligned to +/- strand):  ref base C (37/28);  new base T (2/4);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGACCGCCCA  >  NC_002947/4706048‑4706337
                                                                                                                                                 |                                                                                                                                                
cAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGTGCAGGTCga                                                                                                                                               >  2:356761/1‑148 (MQ=255)
cAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGTGCAGGTCga                                                                                                                                               >  2:351226/1‑148 (MQ=255)
cAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACACTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGt                                                                                                                                               >  2:117543/1‑149 (MQ=255)
 aGCAACCAGTTGCCCACCTGCCCACAGTAAGGCTTGATCCACGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCCACCGTTGACCAGACGCTCGCCAACAAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTc                                                                                                                                              <  6:317712/149‑1 (MQ=255)
 aGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGTGCAGGTCgac                                                                                                                                              <  1:351226/149‑3 (MQ=255)
 aGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGTGCAGGTCgac                                                                                                                                              <  1:356761/149‑3 (MQ=255)
    aaCCAGTTGAGCACCTGAACAAACTAAGGCTTGATCGAGGTAAGCAATTCTGCGCCTTTCTGGTGCAGCGATTGCCTACGTTTGACCAGACGCTCGCCAACGAATGTAATTCCAGCCAGGCAGGCCGGTGCTTCGGGCAGGCCGTCTAc                                                                                                                                           <  2:195315/149‑1 (MQ=255)
         gTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGTCGACGGCCTGcc                                                                                                                                      <  4:129795/149‑2 (MQ=255)
         gTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGTCGACGGCCTGcc                                                                                                                                      <  4:135204/149‑2 (MQ=255)
         gTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCa                                                                                                                                      >  1:58386/1‑149 (MQ=255)
          ttGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCAc                                                                                                                                     <  6:350209/149‑1 (MQ=255)
              ccACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTca                                                                                                                                 >  3:268556/1‑149 (MQ=255)
                     cccAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGAc                                                                                                                                           >  3:206795/1‑132 (MQ=255)
                     cccAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGAc                                                                                                                                           <  4:206795/132‑1 (MQ=255)
                      ccAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGTCACTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGaa                                                                                                                         <  1:117543/149‑1 (MQ=255)
                        aaaGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACg                                                                                                                       >  4:294485/1‑149 (MQ=255)
                         cactAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGA‑CAGACGCTCGCCAACGAAAGGAAAGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGcc                                                                                                                     <  6:42437/146‑1 (MQ=255)
                              aGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTa                                                                                                                  <  3:244871/148‑1 (MQ=255)
                              aGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTACCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTa                                                                                                                  <  5:359442/148‑1 (MQ=255)
                              aGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTATGGTCCAGCGATTGCCATCCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTa                                                                                                                  <  3:223608/148‑1 (MQ=255)
                                gCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCg                                                                                                               >  1:205978/1‑149 (MQ=255)
                                          ggTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAAtgt                                                                                                     >  5:50502/1‑149 (MQ=255)
                                             aaGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGAt                                                                                                  >  8:29961/1‑149 (MQ=255)
                                                cAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCCCCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCCCGCACGAACGCCGTAGCGTCGCGAATGTGATCgg                                                                                               >  3:312362/1‑149 (MQ=255)
                                                cAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATc                                                                                                 <  2:189073/147‑1 (MQ=255)
                                                cAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATc                                                                                                 >  1:189073/1‑147 (MQ=255)
                                                cAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCgg                                                                                               >  3:249778/1‑149 (MQ=255)
                                                        cgcCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGtt                                                                                        <  6:50502/148‑1 (MQ=255)
                                                                tGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAg                                                                               >  2:113372/1‑149 (MQ=255)
                                                                   tCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTcgc                                                                                                           >  4:196192/1‑118 (MQ=255)
                                                                   tCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTcgc                                                                                                           <  3:196192/118‑1 (MQ=255)
                                                                           ttGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCg                                                                    <  7:29961/149‑1 (MQ=255)
                                                                              ccAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCgg                                                                   >  6:135975/1‑147 (MQ=255)
                                                                              ccAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCgg                                                                   >  1:138455/1‑147 (MQ=255)
                                                                                aaCCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGTCCGGTGCATCGGGCAGGCAGTCGACATGCACGTCACGCACGAACGCCGAAGCGTCGCGAAAGTGATCGGCCAGGTTGAACCCCAGCCACCCCAGCGGCAGt                                                               >  3:261311/1‑149 (MQ=255)
                                                                                aaCCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTcgc                                                                                                           >  5:55707/1‑105 (MQ=255)
                                                                                aaCCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTcgc                                                                                                           <  6:55707/105‑1 (MQ=255)
                                                                                 aCCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTg                                                              >  2:19397/1‑149 (MQ=255)
                                                                                 aCCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATAGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTg                                                              >  6:294982/1‑149 (MQ=255)
                                                                                  ccGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGc                                                             <  2:138455/149‑1 (MQ=255)
                                                                                    gTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCa                                                           <  1:358770/149‑1 (MQ=255)
                                                                                        aCCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGa                                                        <  5:135975/148‑1 (MQ=255)
                                                                                            gACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAgg                                                   <  5:240493/149‑1 (MQ=255)
                                                                                              cGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGAt                                                 >  5:172373/1‑149 (MQ=255)
                                                                                                   gCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAgg                                            >  2:283210/1‑149 (MQ=255)
                                                                                                       aCGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCAt                                        >  7:104977/1‑149 (MQ=255)
                                                                                                        cGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATg                                       <  3:6349/149‑1 (MQ=255)
                                                                                                        cGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATg                                       <  2:58386/149‑1 (MQ=255)
                                                                                                        cGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATg                                       <  1:19397/149‑1 (MQ=255)
                                                                                                           aaGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGt                                     >  3:114531/1‑148 (MQ=255)
                                                                                                              gAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCATCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTaag                                  >  4:19988/1‑148 (MQ=255)
                                                                                                              gAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcag                                  >  1:35469/1‑148 (MQ=255)
                                                                                                               aaTGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcagca                                <  7:61162/149‑1 (MQ=255)
                                                                                                               aaTGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcagca                                <  6:172373/149‑1 (MQ=255)
                                                                                                               aaTGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcag                                  >  2:36545/1‑147 (MQ=255)
                                                                                                                     aCCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCcgg                          >  5:245531/1‑149 (MQ=255)
                                                                                                                      cccAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcag                                  <  3:273046/140‑1 (MQ=255)
                                                                                                                      cccAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcag                                  >  4:273046/1‑140 (MQ=255)
                                                                                                                       ccagccagGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCcggcg                        >  2:229783/1‑149 (MQ=255)
                                                                                                                               gCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGt                <  2:35469/149‑1 (MQ=255)
                                                                                                                                  ggTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGtctc             <  3:173942/149‑1 (MQ=255)
                                                                                                                                   gTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTcg            >  5:7016/1‑149 (MQ=255)
                                                                                                                                   gTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTcg            >  2:371458/1‑149 (MQ=255)
                                                                                                                                   gTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTcg            >  1:161981/1‑149 (MQ=255)
                                                                                                                                      cATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGa         <  5:150536/149‑1 (MQ=255)
                                                                                                                                      cATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGa         <  3:294485/149‑1 (MQ=255)
                                                                                                                                       aTCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCCGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGAc        >  7:348276/1‑149 (MQ=255)
                                                                                                                                       aTCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACCCCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGAc        >  4:152404/1‑149 (MQ=255)
                                                                                                                                       aTCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCCGCGtcgagcgca        >  4:208626/1‑142 (MQ=255)
                                                                                                                                           ggCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGACCGcc    >  3:158632/1‑149 (MQ=255)
                                                                                                                                             cAGGCCGTCGCCCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGAACCAGGGCATGGTCAGCAGGCCGGCGGCCAGCGTCGAGCGACCGCCca  >  4:13914/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGACCGCCCA  >  NC_002947/4706048‑4706337

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