Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 436,899 C→A 5.9% W130L (TGG→TTG)  PP_0358 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947436,8990CA 5.9% 272.2 / 5.9 85W130L (TGG→TTG) PP_0358membrane protein
Reads supporting (aligned to +/- strand):  ref base C (41/39);  new base A (2/3);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

CCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGT  >  NC_002947/436755‑437045
                                                                                                                                                |                                                                                                                                                   
ccGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTg                                                                                                                                                 <  4:41865/149‑1 (MQ=255)
 cGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACCCCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTg                                                                                                                                                 >  6:211128/1‑148 (MQ=255)
  gggCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCa                                                                                                                                               <  7:73034/149‑1 (MQ=255)
      agcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAggcg                                                                                                                                           >  2:199786/1‑149 (MQ=255)
        cgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggt                                                                                                                                         <  3:243798/149‑2 (MQ=255)
         gcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggtg                                                                                                                                        <  8:279613/149‑3 (MQ=255)
           agcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTCCTGCGCCACACACACGGCAGCCATGCCCACCGCCAGCCATTGCAGccggcg                                                                                                                                        >  3:216044/1‑147 (MQ=255)
               accgcatccTGCTTACGGATAAGAACAAACAACCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCCGgcg                                                                                                                                  >  8:191847/8‑149 (MQ=255)
               aCCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCGCCGCCTGCAATGGCTggcggtgggcat                                                                                                                                   <  5:143462/148‑8 (MQ=255)
               aCCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCGCCGCCTGCAATGGCTggcggtgggcat                                                                                                                                   >  6:143462/1‑141 (MQ=255)
                ccGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGgcgc                                                                                                                                 <  6:162551/149‑1 (MQ=255)
                 cGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCt                                                                                                                                <  5:211128/149‑1 (MQ=255)
                   gcgcCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCa                                                                                                                              <  2:239918/149‑1 (MQ=255)
                     gcCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCAc                                                                                                                             >  5:160745/1‑148 (MQ=255)
                       cTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCaa                                                                                                                          <  4:80607/149‑1 (MQ=255)
                          cTTACGGATCAGCCCATAGAATCCCAAGCTCAACGCCAGCACCACTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                        <  3:191605/148‑1 (MQ=255)
                           ttACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAAt                                                                                                                      >  5:208705/1‑149 (MQ=255)
                              cGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAg                                                                                                                   >  2:6328/1‑149 (MQ=255)
                               ggATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAg                                                                                                                   >  8:212596/1‑148 (MQ=255)
                                    aGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGcc                                                                                                             <  6:160745/149‑1 (MQ=255)
                                      cccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGcc                                                                                                             >  1:44957/1‑147 (MQ=255)
                                      cccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGcc                                                                                                             >  1:195646/1‑147 (MQ=255)
                                      cccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGcc                                                                                                             >  6:28049/1‑147 (MQ=255)
                                       ccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGTCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAACATGCCcag                                                                                                          >  2:185965/1‑149 (MQ=255)
                                       ccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                          >  2:308483/1‑149 (MQ=255)
                                       ccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                          >  6:215423/1‑149 (MQ=255)
                                         aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcagca                                                                                                        <  1:199786/149‑1 (MQ=255)
                                         aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                          >  8:76/1‑147 (MQ=255)
                                         aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                          >  4:58286/1‑147 (MQ=255)
                                          tAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGACCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCCGCATGCCcagcag                                                                                                       >  8:267066/1‑149 (MQ=255)
                                           aGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCATCAAGAATCAGCATGCCcagcag                                                                                                       >  6:215705/1‑148 (MQ=255)
                                           aGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCaa                                                                                                                          >  1:43668/1‑129 (MQ=255)
                                           aGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCaa                                                                                                                          <  2:43668/129‑1 (MQ=255)
                                             aaGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTa                                                                                                     >  3:170196/1‑148 (MQ=255)
                                                  gAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGAGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCac                                                                                                                   >  6:300138/1‑128 (MQ=255)
                                                       tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                          >  4:142620/1‑149 (MQ=255)
                                                       tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                          >  7:65578/1‑149 (MQ=255)
                                                         aaCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCgg                                                                                        <  1:6328/149‑1 (MQ=255)
                                                          aCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCggg                                                                                       <  2:44957/149‑1 (MQ=255)
                                                               aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGAt                                                                                  <  5:117085/149‑1 (MQ=255)
                                                               aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGAt                                                                                  <  3:58286/149‑1 (MQ=255)
                                                                     aGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATa                                                                            <  3:142620/149‑1 (MQ=255)
                                                                      gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑g                                                                          <  6:166386/149‑1 (MQ=255)
                                                                      gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑g                                                                          <  5:28049/149‑1 (MQ=255)
                                                                      gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑g                                                                          <  5:215423/149‑1 (MQ=255)
                                                                       tGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑Gc                                                                         >  8:215810/1‑149 (MQ=255)
                                                                           aCCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCGCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAg                                                                     >  3:216056/1‑149 (MQ=255)
                                                                             ccACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGc                                                                   <  7:138556/149‑1 (MQ=255)
                                                                             ccACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGc                                                                   >  7:256468/1‑149 (MQ=255)
                                                                                 ggCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGc                                                               <  8:65578/149‑1 (MQ=255)
                                                                                  gCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                          <  3:95049/122‑1 (MQ=255)
                                                                                  gCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                          >  4:95049/1‑122 (MQ=255)
                                                                                      gCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGa                                                          <  7:215810/149‑1 (MQ=255)
                                                                                       cTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGTCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAg                                                         <  8:48120/149‑1 (MQ=255)
                                                                                       cTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAg                                                         <  8:256468/149‑1 (MQ=255)
                                                                                          cccACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCa                                                       <  7:76/148‑1 (MQ=255)
                                                                                           ccACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAACATGCCCAGCAGTACATTGATCAGCGGGTTGa                                                                                   <  7:12900/120‑1 (MQ=255)
                                                                                           ccACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAACATGCCCAGCAGTACATTGATCAGCGGGTTGa                                                                                   >  8:12900/1‑120 (MQ=255)
                                                                                             aCTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACg                                                   <  1:308483/149‑1 (MQ=255)
                                                                                              cTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGcc                                                              <  3:779/137‑1 (MQ=255)
                                                                                              cTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGcc                                                              >  4:779/1‑137 (MQ=255)
                                                                                              cTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCACAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACgg                                                  <  6:208705/149‑1 (MQ=255)
                                                                                                tgggcagccTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGcc                                                <  1:17796/142‑1 (MQ=255)
                                                                                                 gCCACACCTGCTGCGCCACTCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGcc                                                <  4:216056/148‑1 (MQ=255)
                                                                                                 gCCACACCTGCTGCGCCACACCCACGGCAGGCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGcc                                                <  7:267066/148‑1 (MQ=255)
                                                                                                 gCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCg                                               >  3:286022/1‑149 (MQ=255)
                                                                                                  ccACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGTCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAACATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATATGCCCAGGCTGGCCTCGAGCATACGGCCg                                               <  1:185965/149‑1 (MQ=255)
                                                                                                    acacCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCgttg                                            <  4:170196/149‑1 (MQ=255)
                                                                                                       ccTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCATCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTc                                         <  5:215705/149‑1 (MQ=255)
                                                                                                        ctgctgCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCa                                        >  5:75460/1‑149 (MQ=255)
                                                                                                         tgctgcGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAc                                       >  1:41930/1‑149 (MQ=255)
                                                                                                         tgctgcGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAc                                       >  8:102976/1‑149 (MQ=255)
                                                                                                              cgcCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATc                                                                                                                     >  7:274052/1‑67 (MQ=255)
                                                                                                              cgcCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATc                                                                                                                     <  8:274052/67‑1 (MQ=255)
                                                                                                                ccACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACca                                 >  6:98289/1‑148 (MQ=255)
                                                                                                                ccACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACca                                 >  1:284885/1‑148 (MQ=255)
                                                                                                                 cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcaca                               >  1:108410/1‑149 (MQ=255)
                                                                                                                  acacCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATACCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCacac                              >  5:28854/1‑149 (MQ=255)
                                                                                                                      gcACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCTTACGGCCGTTGTTCACCGACCACACGTAg                          <  4:216044/148‑1 (MQ=255)
                                                                                                                          ggCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGt                        >  1:94379/1‑147 (MQ=255)
                                                                                                                              gCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGTCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCa                  <  7:191847/149‑1 (MQ=255)
                                                                                                                                 aTGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGtt               <  4:286022/149‑1 (MQ=255)
                                                                                                                                     ccACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCg           >  2:101190/1‑149 (MQ=255)
                                                                                                                                     ccACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACAGACCACACGTAGGTCAGCCAGTTGCCg           >  8:137618/1‑149 (MQ=255)
                                                                                                                                       gccGCCTGCAATGGCTGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGc         >  8:3546/2‑149 (MQ=255)
                                                                                                                                              gccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATACCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGt  <  6:28854/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGT  >  NC_002947/436755‑437045

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