Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 436,902 T→G 6.2% Q129P (CAA→CCA)  PP_0358 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947436,9020TG 6.2% 202.5 / 4.5 65Q129P (CAA→CCA) PP_0358membrane protein
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base G (2/2);  total (31/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 = 9.51e-01

CCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGCCC  >  NC_002947/436755‑437049
                                                                                                                                                   |                                                                                                                                                    
ccGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCGACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTg                                                                                                                                                     >  6:175407/1‑149 (MQ=255)
  gggCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCa                                                                                                                                                   <  5:351590/149‑1 (MQ=255)
    gcagcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTgg                                                                                                                                                 <  3:99713/149‑1 (MQ=255)
           agcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCa                                                                                                                                          <  5:223681/149‑1 (MQ=255)
                ccGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCCACCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGACGgcgc                                                                                                                                     >  4:67071/1‑149 (MQ=255)
                ccGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGgcgc                                                                                                                                     <  5:15848/149‑1 (MQ=255)
                 cGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCt                                                                                                                                    <  1:193075/149‑1 (MQ=255)
                       cTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCaa                                                                                                                              <  8:210787/149‑1 (MQ=255)
                          cTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGaa                                                                                                                           >  7:262681/1‑149 (MQ=255)
                              cGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAg                                                                                                                       <  6:337811/149‑1 (MQ=255)
                                 aTCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCAt                                                                                                                    <  3:187535/149‑1 (MQ=255)
                                 aTCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCAt                                                                                                                    <  1:303214/149‑1 (MQ=255)
                                        cATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCGCACCAAGAATCAGCATGC‑Ccag                                                                                                             >  1:208985/1‑148 (MQ=255)
                                                  gAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGa                                                                                                  >  1:187231/1‑149 (MQ=255)
                                                  gAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGa                                                                                                  >  6:84271/1‑149 (MQ=255)
                                                    aGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAgg                                                                                                                                                 <  4:254133/101‑1 (MQ=255)
                                                    aGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAgg                                                                                                                                                 >  3:254133/1‑101 (MQ=255)
                                                       tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGc                                                                                             >  6:150392/1‑149 (MQ=255)
                                                          aCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCggg                                                                                          <  5:175407/149‑1 (MQ=255)
                                                               aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAt                                                                                     >  8:277708/1‑149 (MQ=255)
                                                               aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAt                                                                                     <  8:262681/149‑1 (MQ=255)
                                                               aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAt                                                                                     <  7:277708/149‑1 (MQ=255)
                                                                gCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACTTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATg                                                                                    <  1:301094/149‑1 (MQ=255)
                                                                gCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATg                                                                                    <  2:220372/149‑1 (MQ=255)
                                                                 cACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCATACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGc                                                                                                                      <  2:123993/115‑1 (MQ=255)
                                                                  aCCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCAACTGCAGCCATGCCCACAGCCAGCAATTGCAGGCGGCGCAGGCGCTCACGAAGAAACAGCATGC‑CCAGCAGTAAATTGATCAGCGGGTTGATGt                                                                                   >  6:54021/1‑148 (MQ=255)
                                                                     aGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATa                                                                               <  5:150392/149‑1 (MQ=255)
                                                                         aTACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGcc                                                                            >  5:217806/1‑148 (MQ=255)
                                                                             ccACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGc                                                                       >  1:359390/1‑149 (MQ=255)
                                                                             ccACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCACGCCCACCGCCTGCAATGGCTggcggtgg                                                                                                                                           >  2:105661/1‑79 (MQ=37)
                                                                              cACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAgg                                                                        <  8:317609/147‑1 (MQ=255)
                                                                              cACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAgg                                                                        >  7:317609/1‑147 (MQ=255)
                                                                                    aGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTc                                                                <  8:278110/149‑1 (MQ=255)
                                                                                        tGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAg                                                             >  4:4065/1‑148 (MQ=255)
                                                                                          cccACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAAGTAATAGCCCAGGCTGGCCTCGAGc                                                            <  2:208985/147‑1 (MQ=255)
                                                                                            cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAt                                                                                     <  1:356973/120‑1 (MQ=255)
                                                                                            cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAt                                                                                     >  2:356973/1‑120 (MQ=255)
                                                                                               ttGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGATGGCCTCGAGCATACGGc                                                     >  6:218252/1‑149 (MQ=255)
                                                                                                    acacCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgttg                                                <  5:84271/149‑1 (MQ=255)
                                                                                                       actgctgCGCCACACCCCCGGCATCCATGCCCACCGCCAGACATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATTC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAg                                                                              <  6:178660/115‑1 (MQ=255)
                                                                                                           ctgcGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAcc                                          <  2:187231/148‑1 (MQ=255)
                                                                                                              cgcCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGa                                        <  8:114595/147‑1 (MQ=255)
                                                                                                                ccACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCTGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGCGCATACGGCCGTTGTTCACCGACcac                                    <  1:8680/149‑1 (MQ=255)
                                                                                                                ccACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACca                                     >  1:238092/1‑148 (MQ=255)
                                                                                                                 cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcaca                                   >  7:262841/1‑149 (MQ=255)
                                                                                                                 cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcac                                    >  6:210786/1‑148 (MQ=255)
                                                                                                                 cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcac                                    <  5:210786/148‑1 (MQ=255)
                                                                                                                    acCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGt                                <  2:359390/149‑1 (MQ=255)
                                                                                                                      ccACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAg                              <  7:88151/149‑1 (MQ=255)
                                                                                                                       cACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGGTGTTCACCGCCCACACGTAgg                             >  4:85741/1‑149 (MQ=255)
                                                                                                                          ggCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCa                          <  6:217806/149‑1 (MQ=255)
                                                                                                                            cAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGcc                                                    >  4:88344/1‑121 (MQ=255)
                                                                                                                             aGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGc                        >  6:177221/1‑148 (MQ=255)
                                                                                                                               ccATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCGCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCa                      >  1:279316/1‑149 (MQ=255)
                                                                                                                                   gCCCACCGCCAGCCAATGCAGGCGGCGCAGGCGCTCAGCAAGAATGTGCATGC‑CCAGCAGGACAGTGATCAGCGGGTTGAAGTAATAGCCCAGGCTCGCCGCGAGCATACGGCCGTTGGTCACCGCCCACACGTAGGGCAGCCGGTTGc                 >  4:258280/1‑149 (MQ=255)
                                                                                                                                    tgcgccGCCTGCAATGGCTGGCGGTGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGcc                >  4:234620/5‑149 (MQ=255)
                                                                                                                                    cccACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGcc                >  5:310851/1‑149 (MQ=255)
                                                                                                                                      cACCGCCAGACATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCATCGGGTTGATGTAATAGCCCAGGATGGCCTCGAGCATACGGACGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCgg              >  8:108656/1‑149 (MQ=255)
                                                                                                                                       gccGCCTGCAATGGCTGGCGGTGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACTCGTAGGTCAGCCAGTTGCCGGc             <  7:227252/148‑1 (MQ=255)
                                                                                                                                          gccagccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGAt          <  5:177221/149‑1 (MQ=255)
                                                                                                                                            cagccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGAt          >  1:227413/1‑147 (MQ=255)
                                                                                                                                              gccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAgg                                                                        <  7:150264/83‑1 (MQ=255)
                                                                                                                                              gccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAgg                                                                        >  8:150264/1‑83 (MQ=255)
                                                                                                                                               ccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcaca                                   <  4:340856/119‑1 (MQ=255)
                                                                                                                                               ccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcaca                                   >  3:340856/1‑119 (MQ=255)
                                                                                                                                                  ttGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGccc  <  3:4065/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGCCC  >  NC_002947/436755‑437049

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