Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,908,416 A→G 5.6% N122S (AAC→AGC)  PP_3450 → TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,908,4160AG 5.6% 258.2 / 2.4 89N122S (AAC→AGC) PP_3450TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (42/41);  new base G (3/2);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGC  >  NC_002947/3908273‑3908564
                                                                                                                                               |                                                                                                                                                      
cGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACcgcg                                                                                                                                                   >  5:299912/1‑149 (MQ=255)
cGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACcgcg                                                                                                                                                   >  8:335720/1‑149 (MQ=255)
      aGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACcgc                                                                                                                                                    <  6:122426/142‑1 (MQ=255)
      aGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACcgc                                                                                                                                                    >  5:122426/1‑142 (MQ=255)
               gcctgccCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGcc                                                                                                                                    >  2:202721/1‑149 (MQ=255)
                    cccAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccac                                                                                                                               <  5:369493/149‑1 (MQ=255)
                    cccAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccac                                                                                                                               >  1:350153/1‑149 (MQ=255)
                     ccAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccacc                                                                                                                              <  6:299912/149‑1 (MQ=255)
                     ccAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGCGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCACCGGCCAACCAGACGCGGCTGGCCCTGGACAACCGCCCCATTGCCCTCGGccaccac                                                                                                                               >  4:126968/1‑148 (MQ=255)
                      cAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccacca                                                                                                                             >  5:30521/1‑149 (MQ=255)
                           cTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAgg                                                                                                                        <  8:305981/149‑1 (MQ=255)
                            tGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGTTGGCCATGGCCAGCcgcgtctg                                                                                                                                               <  7:203563/125‑5 (MQ=255)
                            tGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGTTGGCCATGGCCAGCcgcgtctg                                                                                                                                               >  8:203563/1‑121 (MQ=255)
                            tGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGc                                                                                                                       >  4:106499/1‑149 (MQ=255)
                            tGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGc                                                                                                                       >  8:30899/1‑149 (MQ=255)
                              gCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGccc                                                                                                                     <  6:30521/149‑1 (MQ=255)
                              gCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGcc                                                                                                                      >  3:79741/1‑148 (MQ=255)
                                aaCGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCaa                                                                                                                   >  8:137977/1‑149 (MQ=255)
                                aaCGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCACCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCAGCCACCCCCAGCAGGCCCaa                                                                                                                   >  3:86138/1‑149 (MQ=255)
                                   gCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgc                                                                                                                <  7:335720/149‑1 (MQ=255)
                                        gCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctggc‑‑t                                                                                                         >  2:168707/1‑149 (MQ=255)
                                        gCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctgg                                                                                                             <  4:34358/147‑1 (MQ=255)
                                        gCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctgg                                                                                                             <  3:201026/147‑1 (MQ=255)
                                         ccAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctggc‑‑tg                                                                                                        <  2:350153/149‑1 (MQ=255)
                                         ccAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctggc‑‑tg                                                                                                        <  4:79741/149‑1 (MQ=255)
                                          cAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTTGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGc                                                                                                       >  4:259536/1‑149 (MQ=255)
                                           aaTGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCa                                                                                                      <  1:187422/149‑1 (MQ=255)
                                             tGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGt                                                                                                    >  7:208671/1‑149 (MQ=255)
                                                cTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGc                                                                                                 <  1:168707/149‑1 (MQ=255)
                                                    atCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTTGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTAcc                                                                                             <  3:259536/149‑1 (MQ=255)
                                                      ccAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACAGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGc                                                                                           >  6:218357/1‑149 (MQ=255)
                                                         gCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCt                                                                                         >  7:294921/1‑148 (MQ=255)
                                                           aaGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctggc                                                                                                            <  2:175992/129‑1 (MQ=255)
                                                                   gaAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTgc                                                                               >  5:351232/1‑148 (MQ=255)
                                                                   gaAAAGGCCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTgc                                                                               <  6:351232/148‑1 (MQ=255)
                                                                           ccTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTgca                                                                              <  5:37932/141‑1 (MQ=255)
                                                                           ccTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTgca                                                                              >  6:37932/1‑141 (MQ=255)
                                                                           ccTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTgc                                                                               >  5:298377/1‑140 (MQ=255)
                                                                           ccTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTgc                                                                               <  6:298377/140‑1 (MQ=255)
                                                                              gACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGTCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGag                                                                    >  6:39287/1‑148 (MQ=255)
                                                                               aCTGCGCTGCTGGCGCACGCCCCCTCGCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCTGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagc                                                                   <  3:126968/148‑1 (MQ=255)
                                                                               aCTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCACCCGCGCCAAAGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGag                                                                    >  8:280738/1‑147 (MQ=255)
                                                                               aCTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGag                                                                    <  1:202721/147‑1 (MQ=255)
                                                                                 tGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGTTGGCCATGGCCAGCcgc                                                                                                                                                    <  3:357505/67‑1 (MQ=255)
                                                                                 tGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGTTGGCCATGGCCAGCcgc                                                                                                                                                    >  4:357505/1‑67 (MQ=255)
                                                                                 tGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagcagc                                                                >  1:165045/1‑149 (MQ=255)
                                                                                     ctgctgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGTCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCa                                                            <  5:39287/149‑1 (MQ=255)
                                                                                       gctgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAg                                                           >  3:353177/1‑148 (MQ=255)
                                                                                        ctgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAgcg                                                               <  2:259586/143‑1 (MQ=255)
                                                                                        ctgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAgcg                                                               >  1:259586/1‑143 (MQ=255)
                                                                                         tgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCaa                                                        >  1:142052/1‑149 (MQ=255)
                                                                                           gCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCCGGCCCAACGCTGCC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCACCAACt                                                      >  5:74364/1‑149 (MQ=255)
                                                                                           gCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACt                                                      <  7:30899/149‑1 (MQ=255)
                                                                                             tcccGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGAGCCATTGCCGTAGGCCACCCCCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCTAGCAGCGCCAGCAACttt                                                    <  5:214320/146‑1 (MQ=255)
                                                                                             ccACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACttt                                                    >  4:192926/1‑149 (MQ=255)
                                                                                             ccACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACttt                                                    >  3:231892/1‑149 (MQ=255)
                                                                                              cACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAAc                                                                                  <  7:318050/118‑1 (MQ=255)
                                                                                              cACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTg                                                   <  3:192926/149‑1 (MQ=255)
                                                                                              cACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCCGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCACGCCCCACGCTCGC‑‑AGCAGTCGCTACCGCCAGATCAAc                                                                                  >  8:318050/1‑118 (MQ=255)
                                                                                               aCGCCCCCTCCCCCGCAGCGGCCAACCCGACGCGGCTGGCCATGGCCAAACGCGCCCTTGCCGGCGGACCCCACCCGACGGCCCACCGCAGGC‑‑TGCCGTCGCAACCGCCTGCTCAACGCCAGCCACCCGCGCCCCCCCAGccaagtg                                                    >  6:214320/1‑141 (MQ=255)
                                                                                               aCGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTgg                                                  >  6:299464/1‑149 (MQ=255)
                                                                                               aCGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTgg                                                  >  8:80222/1‑149 (MQ=255)
                                                                                                cGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGc                                                 <  2:165045/149‑1 (MQ=255)
                                                                                                   ccccTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCTCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAAc                                              >  1:230932/1‑149 (MQ=255)
                                                                                                   ccccTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAAc                                              <  4:353177/149‑1 (MQ=255)
                                                                                                     ccTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTg                                            <  5:218357/149‑1 (MQ=255)
                                                                                                     ccTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTg                                            <  3:106499/149‑1 (MQ=255)
                                                                                                              cAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGaa                                   <  2:230932/149‑1 (MQ=255)
                                                                                                                gCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAgg                                 <  8:208671/149‑1 (MQ=255)
                                                                                                                gCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAgg                                 <  7:137977/149‑1 (MQ=255)
                                                                                                                 cGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGc                                >  2:341096/1‑149 (MQ=255)
                                                                                                                 cGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGc                                >  2:346430/1‑149 (MQ=255)
                                                                                                                        ccAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAg                                                           >  6:21725/1‑115 (MQ=255)
                                                                                                                        ccAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAg                                                           <  5:21725/115‑1 (MQ=255)
                                                                                                                        ccAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATg                         >  3:151745/1‑149 (MQ=255)
                                                                                                                        aatggcgcgGTTGGCCATGGCCAGCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATg                         >  3:33714/6‑149 (MQ=255)
                                                                                                                          aGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCa                       <  7:80222/149‑1 (MQ=255)
                                                                                                                                ggCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTg                 <  6:309581/149‑1 (MQ=255)
                                                                                                                                  cTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCGATGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTg                 >  1:285104/1‑149 (MQ=255)
                                                                                                                                   tggccatggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTg              >  8:264810/1‑149 (MQ=255)
                                                                                                                                       catggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTgcggc          <  5:299464/149‑1 (MQ=255)
                                                                                                                                        atggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTgcggcg         >  8:111361/1‑149 (MQ=255)
                                                                                                                                         tggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGgcgc        <  1:341096/149‑1 (MQ=255)
                                                                                                                                         tggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGgcgc        <  1:346430/149‑1 (MQ=255)
                                                                                                                                          ggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCaa                                                                                   <  6:175455/73‑1 (MQ=255)
                                                                                                                                          ggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCaa                                                                                   >  5:175455/1‑73 (MQ=255)
                                                                                                                                           gccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCt      >  5:218564/1‑149 (MQ=255)
                                                                                                                                             caACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCGCAAGGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGa    <  6:74364/149‑1 (MQ=255)
                                                                                                                                             caACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGa    <  5:16182/149‑1 (MQ=255)
                                                                                                                                              aaCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagc                                                                   <  6:116875/85‑1 (MQ=255)
                                                                                                                                              aaCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagc                                                                   >  5:116875/1‑85 (MQ=255)
                                                                                                                                              aaCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAg   >  7:116511/1‑149 (MQ=255)
                                                                                                                                               aCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGc  <  2:285104/149‑1 (MQ=255)
                                                                                                                                               aCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGc  <  4:151745/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
CGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGC‑‑TGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGC  >  NC_002947/3908273‑3908564

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: