Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 436,905 A→T 5.9% L128Q (CTG→CAG)  PP_0358 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947436,9050AT 5.9% 258.0 / 4.8 84L128Q (CTG→CAG) PP_0358membrane protein
Reads supporting (aligned to +/- strand):  ref base A (41/38);  new base T (2/3);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

GGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGT  >  NC_002947/436757‑437045
                                                                                                                                                    |                                                                                                                                             
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                                                                                                             >  6:28049/1‑147 (MQ=255)
                                    cccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGcc                                                                                                             >  1:195646/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                                                                                  <  3:58286/149‑1 (MQ=255)
                                                             aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGAt                                                                                  <  5:117085/149‑1 (MQ=255)
                                                                   aGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATa                                                                            <  3:142620/149‑1 (MQ=255)
                                                                    gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑g                                                                          <  5:215423/149‑1 (MQ=255)
                                                                    gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑g                                                                          <  6:166386/149‑1 (MQ=255)
                                                                    gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑g                                                                          <  5:28049/149‑1 (MQ=255)
                                                                     tGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑Gc                                                                         >  8:215810/1‑149 (MQ=255)
                                                                         aCCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCGCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAg                                                                     >  3:216056/1‑149 (MQ=255)
                                                                           ccACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGc                                                                   >  7:256468/1‑149 (MQ=255)
                                                                           ccACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGc                                                                   <  7:138556/149‑1 (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                                 >  1:284885/1‑148 (MQ=255)
                                                                                                              ccACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACca                                 >  6:98289/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)
                                                                                                                                                   cAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCgttgtt                                          >  4:266683/1‑102 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGT  >  NC_002947/436757‑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: ‑