Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,334,844 T→G 8.9% Q538P (CAG→CCG) ‡ PP_2052 ← bifunctional sugar‑phosphatase/mannitol‑1‑phosphate 5‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,334,8440TG 8.9% 196.8 / 6.5 65Q538P (CAG→CCG) ‡PP_2052bifunctional sugar‑phosphatase/mannitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (32/27);  new base G (2/4);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGAGGTTTTGCAGCGCTCGAGAAATGCCTCGGCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCG  >  NC_002947/2334697‑2334990
                                                                                                                                                    |                                                                                                                                                  
aaCGAGGTTTTGCAGCGCTCGAGAAATGCCTCGGCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTg                                                                                                                                                   >  5:273072/1‑149 (MQ=255)
            cAGCGCTCGAGAAATGCCTCGGCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGc                                                                                                                                       <  6:273072/149‑1 (MQ=255)
                gcTCGAGAAATGCCTCGGCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCTGGGTCGGCCAGGGCATgg                                                                                                                                   <  6:309/149‑1 (MQ=255)
                 cTCGAGAAATGCCTCGGCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCTGGGTCGGCCAGGGCATggg                                                                                                                                  <  5:288468/149‑2 (MQ=255)
                  tCGAGAAATGCCTCGGCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCTGGGTCGGCCAGGGCATgg                                                                                                                                   <  1:274720/147‑1 (MQ=255)
                  tCGAGAAATGCCTCGGCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCTGGGTCGGCCAGGGCATgg                                                                                                                                   >  2:274720/1‑147 (MQ=255)
                  tCGAGAAATGCCTCGGCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGcc                                                                                                                                 >  5:209017/1‑149 (MQ=255)
                    cagaAATGACTCGGCAAAGCGCCTGACAAACTCGGCCATCTGTGGGTATTAGGCCACCAGTGCTGGCTCGAACTCCTGGCGAATCAGGCCCTCGGCCAGTTAACTCACACGCTCATC‑GCCCATGCCCTGGCCGACCCAGGCATGTCCca                                                                                                                               >  6:70546/2‑149 (MQ=255)
                           gCCTCGGCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCt                                                                                                                        >  2:176559/1‑149 (MQ=255)
                                ggCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGc                                                                                                                     <  8:90032/147‑1 (MQ=255)
                                ggCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGc                                                                                                                     <  4:275273/147‑1 (MQ=255)
                                ggCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGc                                                                                                                     <  3:90699/147‑1 (MQ=255)
                                 gCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTa                                                                                                                  >  6:100640/1‑149 (MQ=255)
                                             accaccTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAAt                                                                                                      >  1:317174/1‑149 (MQ=255)
                                               caccTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATgg                                                                                                    >  3:160279/1‑149 (MQ=255)
                                                       ccATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCggg                                                                                            <  5:43769/149‑1 (MQ=255)
                                                         aTCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGcc                                                                                          >  4:80499/1‑149 (MQ=255)
                                                          tCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCg                                                                                         >  6:164236/1‑149 (MQ=255)
                                                              tgGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTcc                                                                                     <  1:176559/149‑1 (MQ=255)
                                                               gggTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCca                                                                                    >  1:135515/1‑149 (MQ=255)
                                                                 gTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAAt                                                                                                      <  4:38996/129‑1 (MQ=255)
                                                                 gTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAAt                                                                                                      >  3:38996/1‑129 (MQ=255)
                                                                      cGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCg                                                                              <  2:317174/148‑1 (MQ=255)
                                                                      cGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCg                                                                              <  4:160279/148‑1 (MQ=255)
                                                                       ggCCACCAGTGCGGGCGCGACCTGCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCgg                                                                             <  6:209017/148‑1 (MQ=255)
                                                                       ggCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGt                                                                            >  5:318698/1‑149 (MQ=255)
                                                                       ggCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGt                                                                            >  5:158134/1‑149 (MQ=255)
                                                                          caccaGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCt                                                                         <  5:164236/149‑1 (MQ=255)
                                                                           accaGTGCGGTCGCGACCTCCTGGAGAATCAGCCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCTTGCCCTGGCCGACGCAGGCATGGCCCATCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCTTTCCACAGGCGGTTCtt                                                                        <  7:66355/149‑1 (MQ=255)
                                                                           accaGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCtt                                                                        >  7:76494/1‑149 (MQ=255)
                                                                             caGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTtga                                                                      >  7:127530/1‑149 (MQ=255)
                                                                                  cGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATgatc                                                                 <  5:57601/149‑1 (MQ=255)
                                                                                   gggCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATCGGCCCATGCCCTGGCCGCCCCAGGCATGGCCCAGCAGGATGCCGTACCAGGCACCAATGGCGTGCGGGCCGTTCCACAGGCGGTTGTtgatgagc                                                                 >  1:89687/1‑147 (MQ=255)
                                                                                   gggCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATgatct                                                                <  5:85828/149‑1 (MQ=255)
                                                                                       gcgACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATgg                                                                                                    <  5:70632/109‑1 (MQ=255)
                                                                                       gcgACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATgg                                                                                                    >  6:70632/1‑109 (MQ=255)
                                                                                          aCCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATgatctggatctg                                                         <  2:135515/149‑1 (MQ=255)
                                                                                                  gCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGt                                                 >  5:106249/1‑149 (MQ=255)
                                                                                                  gCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGt                                                 >  6:16208/1‑149 (MQ=255)
                                                                                                    gAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGt                                                 >  2:90484/1‑147 (MQ=255)
                                                                                                    gAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCt                                               <  6:106249/149‑1 (MQ=255)
                                                                                                      aTCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTgg                                             <  3:80499/149‑1 (MQ=255)
                                                                                                          ggCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACt                                         >  5:68139/1‑149 (MQ=255)
                                                                                                             gcTCTGCGATTTCACTCACTCGCGCATC‑GCCCATGCCCTGGCCGACGCAGGCATGGCCCAGCATGCTGGCGTACCAGGCCACAATGTAGTGCGGGACGTTCCACAGGCGGTTCTTGATGATATGGATCTGGGCAATGTCTGGCACTGTc                                      <  5:70546/149‑1 (MQ=255)
                                                                                                                      gTTCGCTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGGGTACCAGGCCACAATGGCGTGCGGGCCGTGCCAGAGGCGGTGCTTGAGGGGCTGGCTCGGGGGAATGTCTGGCACTGGCTTCACCTgg                             >  2:119339/1‑149 (MQ=255)
                                                                                                                          actcacCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGgcg                           >  4:136136/1‑147 (MQ=255)
                                                                                                                          actcacCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTACACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGtcg                           >  8:146136/1‑147 (MQ=255)
                                                                                                                            tcacCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAgg                       >  2:10859/1‑149 (MQ=255)
                                                                                                                              acCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGTCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGc                      <  1:90484/148‑1 (MQ=255)
                                                                                                                                ccGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATgatctggatct                                                          >  8:52482/1‑110 (MQ=255)
                                                                                                                                ccGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATgatctggatct                                                          <  7:52482/110‑1 (MQ=255)
                                                                                                                                ccGCGCATC‑GCCCATGCACTGGCCGACCCAGGCATGGCACAGCAGGCTGGCGTAAAAGGCCACAATGGCGTGCGGGCCGTTCCACAGTCGGTTCTTGATGATCTGGATCTGGGCAATGTGTAGCACTTTCTTAACCTGGCGCAGGCGtt                   >  6:327610/1‑149 (MQ=255)
                                                                                                                                  gcgcATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCAc                                 >  2:58783/1‑133 (MQ=255)
                                                                                                                                  gcgcATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCAc                                 <  1:58783/133‑1 (MQ=255)
                                                                                                                                  gcgcATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCg                 <  2:89687/149‑1 (MQ=255)
                                                                                                                                      aTC‑GACCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAg               >  7:210826/1‑147 (MQ=255)
                                                                                                                                        c‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAg               >  5:322332/1‑145 (MQ=255)
                                                                                                                                         gggccATGCCTGGGTCGGCCAGGGCATGGGCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTgg                                             <  7:257483/112‑1 (MQ=37)
                                                                                                                                         gggccATGCCTGGGTCGGCCAGGGCATGGGCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTgg                                             >  8:257483/4‑115 (MQ=37)
                                                                                                                                          gCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGAt          <  3:136136/149‑1 (MQ=255)
                                                                                                                                              aTGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCt      >  7:137418/1‑149 (MQ=255)
                                                                                                                                               tGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTg     <  5:100640/149‑1 (MQ=255)
                                                                                                                                               tGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTg     <  6:28732/149‑1 (MQ=255)
                                                                                                                                                 cccTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGcc   <  6:68139/149‑1 (MQ=255)
                                                                                                                                                  ccTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGGTCCACAGGCGGTTATTGAAGCTCTGGATCTGGGCATTGTCTGGCAATCTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCcg  >  4:245613/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AACGAGGTTTTGCAGCGCTCGAGAAATGCCTCGGCAAAGCGGCTGACCACCTCGGCCATCTGTGGGTATTCGGCCACCAGTGCGGGCGCGACCTCCTGGCGAATCAGGCGCTCGGCCAGTTCACTCACCCGCGCATC‑GCCCATGCCCTGGCCGACCCAGGCATGGCCCAGCAGGCTGGCGTACCAGGCCACAATGGCGTGCGGGCCGTTCCACAGGCGGTTCTTGATGATCTGGATCTGGGCAATGTCTGGCACTGTCTTCACCTGGCGCAGGCGTTCGAGCAGATCGCTGCCG  >  NC_002947/2334697‑2334990

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