Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,761,927 A→T 7.5% N359Y (AAC→TAC) ‡ PP_1567 → capsid protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,761,9270AT 7.5% 243.1 / 7.9 80N359Y (AAC→TAC) ‡PP_1567capsid protein
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base T (4/2);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACG  >  NC_002947/1761786‑1762074
                                                                                                                                              |                                                                                                                                                    
gCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑cgc                                                                                                                                              <  7:81734/149‑1 (MQ=255)
  ggACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCAT‑cgcgc                                                                                                                                            <  8:241231/149‑2 (MQ=255)
  ggACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCgg                                                                                                                                            >  8:279504/1‑149 (MQ=255)
      aTCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGcc                                                                                                                                          >  2:139764/1‑147 (MQ=255)
       tCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAg                                                                                                                                       <  5:96585/149‑1 (MQ=255)
              gACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCAT‑cgcgc                                                                                                                                            >  6:23319/1‑136 (MQ=255)
               aCCAAGGACGCCAACGTCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGTAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTc                                                                                                                                 <  4:388348/147‑1 (MQ=255)
               aCCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTc                                                                                                                                 <  8:39665/147‑1 (MQ=255)
               aCCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTc                                                                                                                                 <  7:74222/147‑1 (MQ=255)
                 cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCAT‑cgcgcctacca                                                                                                                                      >  5:19636/1‑133 (MQ=255)
                 cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCAT‑cgcgcctacca                                                                                                                                      <  6:19636/140‑8 (MQ=255)
                 cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGAcc                                                                                                                             >  8:48428/1‑149 (MQ=255)
                 cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGAcc                                                                                                                             >  2:64926/1‑149 (MQ=255)
                 cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGAc                                                                                                                              >  5:72429/1‑148 (MQ=255)
                 cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGAc                                                                                                                              <  6:72429/148‑1 (MQ=255)
                   aGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACcgc                                                                                                                           >  1:407896/1‑149 (MQ=255)
                          cAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAg                                                                                                                                       >  1:294345/1‑131 (MQ=255)
                          cAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAg                                                                                                                                       <  2:294345/130‑1 (MQ=255)
                           aaCGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGc                                                                                                                   >  6:285337/1‑149 (MQ=255)
                           aaCGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGc                                                                                                                   <  2:231643/149‑1 (MQ=255)
                            aCGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCg                                                                                                                  <  5:285337/149‑1 (MQ=255)
                                 tCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGAGGCTGt                                                                                                             >  7:240700/1‑149 (MQ=255)
                                      cATCTGGGTGAATGTCCAGGATGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAg                                                                                                         <  5:207502/148‑1 (MQ=255)
                                             gTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCAc                                                                                                 >  3:35407/1‑149 (MQ=255)
                                             gTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCAc                                                                                                 <  8:265840/149‑1 (MQ=255)
                                                 aTGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGaa                                                                                             <  7:269846/149‑1 (MQ=255)
                                                    tCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGAc                                                                                          >  6:256888/1‑149 (MQ=255)
                                                    tCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGAc                                                                                          >  5:248985/1‑149 (MQ=255)
                                                    tCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGAc                                                                                          >  8:35283/1‑149 (MQ=255)
                                                      cAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACgg                                                                                        >  6:198093/1‑149 (MQ=255)
                                                       aGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGt                                                                                       <  2:407896/149‑1 (MQ=255)
                                                       aGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGt                                                                                       >  4:242117/1‑149 (MQ=255)
                                                        ggAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTg                                                                                      <  6:18775/149‑1 (MQ=255)
                                                         gAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGa                                                                                     >  5:223107/1‑149 (MQ=255)
                                                          aaGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGAc                                                                                    <  5:256888/149‑1 (MQ=255)
                                                                gcgTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGGTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGAGCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTc                                                                              <  7:249396/149‑1 (MQ=255)
                                                                            cTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCACGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGa                                                                                     >  6:32485/1‑131 (MQ=255)
                                                                            cTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCACGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGa                                                                                     <  5:32485/131‑1 (MQ=255)
                                                                               tgGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATg                                                                  <  5:114325/147‑1 (MQ=255)
                                                                               tgGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGt                                                                >  6:232098/1‑149 (MQ=255)
                                                                                 gAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGa                                                              >  7:46931/1‑149 (MQ=255)
                                                                                 gAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGa                                                              >  6:74649/1‑149 (MQ=255)
                                                                                         cGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCcgc                                                      >  6:271963/1‑149 (MQ=255)
                                                                                            tggtggACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGcc                                                   <  7:279504/149‑1 (MQ=255)
                                                                                               tggACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgag                                                >  7:70100/1‑149 (MQ=255)
                                                                                                ggACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgag                                                >  8:97971/1‑148 (MQ=255)
                                                                                                 gACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgagga                                              <  5:198093/149‑1 (MQ=255)
                                                                                                 gACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgag                                                >  2:152119/1‑147 (MQ=255)
                                                                                                  aCAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgaggag                                             >  6:342113/1‑149 (MQ=255)
                                                                                                     gCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAgcgc                                          <  6:223107/149‑1 (MQ=255)
                                                                                                         cgcgGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctgg                                      <  1:64926/149‑1 (MQ=255)
                                                                                                            ggTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctggcgc                                   <  5:271963/149‑1 (MQ=255)
                                                                                                             gTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctggcgct                                  <  1:152119/149‑1 (MQ=255)
                                                                                                                ccAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATGCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctggcgctgg                                >  1:60208/1‑149 (MQ=255)
                                                                                                                 cAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctggcgctggcg                              <  5:256679/149‑1 (MQ=255)
                                                                                                                  aGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCgg                             >  1:202968/1‑149 (MQ=255)
                                                                                                                  aGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGAGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctgg                                      >  3:52804/1‑140 (MQ=255)
                                                                                                                  aGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGAGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctgg                                      <  4:52804/140‑1 (MQ=255)
                                                                                                                         gAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCg                      >  8:81449/1‑149 (MQ=255)
                                                                                                                         gAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCg                      >  8:339999/1‑149 (MQ=255)
                                                                                                                         gAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGCAGCGCCTGGCGCTGGCGGTCTACCg                      >  1:321495/1‑149 (MQ=255)
                                                                                                                              ccTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCa                  <  7:97971/148‑1 (MQ=255)
                                                                                                                              ccTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCa                  <  8:240700/148‑1 (MQ=255)
                                                                                                                               cTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTAGACCGccagcg                <  6:115740/149‑5 (MQ=255)
                                                                                                                                tGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCagag               >  8:63108/1‑149 (MQ=255)
                                                                                                                                 ctgggccgcgATGTTCAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGt              >  3:376822/7‑149 (MQ=255)
                                                                                                                                  gTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCgcgc                                                    <  4:124878/110‑1 (MQ=255)
                                                                                                                                  gTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCgcgc                                                    >  3:124878/1‑110 (MQ=255)
                                                                                                                                     ggCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGtcgtt          <  8:70100/149‑1 (MQ=255)
                                                                                                                                      ccgcgATGTTCAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGtcgttc         >  6:256953/2‑149 (MQ=255)
                                                                                                                                       cgcgATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGtcgttcg        <  4:35407/149‑1 (MQ=255)
                                                                                                                                        gcgATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTg                                                               >  8:221018/1‑93 (MQ=255)
                                                                                                                                        gcgATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTg                                                               <  7:221018/93‑1 (MQ=255)
                                                                                                                                        gcgATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGAGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCagag               <  6:202951/141‑1 (MQ=255)
                                                                                                                                        gcgATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGAGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCagag               >  5:202951/1‑141 (MQ=255)
                                                                                                                                           aTGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCa    <  7:63108/149‑1 (MQ=255)
                                                                                                                                            tGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCa                  <  8:27712/134‑1 (MQ=255)
                                                                                                                                            tGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCa                  >  7:27712/1‑134 (MQ=255)
                                                                                                                                            tGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCAc   <  7:339999/149‑1 (MQ=255)
                                                                                                                                             gAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACg  >  1:391991/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACG  >  NC_002947/1761786‑1762074

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