Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,706,187 G→T 5.6% P102H (CCC→CAC)  PP_4162 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,706,1870GT 5.6% 219.6 / 3.4 71P102H (CCC→CAC) PP_4162membrane protein
Reads supporting (aligned to +/- strand):  ref base G (35/32);  new base T (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

CTGCGCGCCAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGACCGCC  >  NC_002947/4706040‑4706335
                                                                                                                                                   |                                                                                                                                                    
ctgcgcgcGAGAAACCAGTTGCCCACCTGCCCAAAGTAATGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCggg                                                                                                                                                     <  5:295217/149‑1 (MQ=255)
   cgcgcCAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGc                                                                                                                                                    <  4:20126/147‑1 (MQ=255)
   cgcgcCAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAg                                                                                                                                                  <  4:205281/149‑1 (MQ=255)
        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                                                                                                                                         <  4:206795/132‑1 (MQ=255)
                             cccAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGAc                                                                                                                                         >  3:206795/1‑132 (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                                                                 >  1:138455/1‑147 (MQ=255)
                                                                                      ccAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCgg                                                                 >  6:135975/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                              <  6:172373/149‑1 (MQ=255)
                                                                                                                       aaTGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcagca                              <  7:61162/149‑1 (MQ=255)
                                                                                                                       aaTGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcag                                >  2:36545/1‑147 (MQ=255)
                                                                                                                             aCCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCcgg                        >  5:245531/1‑149 (MQ=255)
                                                                                                                              cccAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcag                                >  4:273046/1‑140 (MQ=255)
                                                                                                                              cccAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcag                                <  3:273046/140‑1 (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       <  3:294485/149‑1 (MQ=255)
                                                                                                                                              cATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGa       <  5:150536/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)
                                                                                                                                                   |                                                                                                                                                    
CTGCGCGCCAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGCCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGACCGCC  >  NC_002947/4706040‑4706335

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: ‑