Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,858,761 C→A 5.2% P237H (CCC→CAC)  PP_5133 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,858,7610CA 5.2% 1059.0 / 30.4 306P237H (CCC→CAC) PP_5133membrane protein
Reads supporting (aligned to +/- strand):  ref base C (168/122);  new base A (8/8);  total (176/130)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 309 total aligned reads displayed.

ATGCTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGTGGCGCTGGC  >  NC_002947/5858613‑5858910
                                                                                                                                                    |                                                                                                                                                     
aTGCTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGccc                                                                                                                                                      <  1:2673824/150‑1 (MQ=255)
   cTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCcgac                                                                                                                                                   >  1:3170480/1‑150 (MQ=255)
    tCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCcgacg                                                                                                                                                  >  2:227489/1‑150 (MQ=255)
    tCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCcgacg                                                                                                                                                  >  1:3236604/1‑150 (MQ=255)
      aCCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGTCGGGCACGGCGaa                                                                                                                                                <  2:2135644/150‑1 (MQ=255)
       ccGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGTCGGGCACGGCGAAg                                                                                                                                               <  1:2492175/150‑1 (MQ=255)
       ccGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGTCGGGCACGGCGAAg                                                                                                                                               <  1:2043645/150‑1 (MQ=255)
          cgcTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGGTTTGGCTTTGTCTGGGTAACCTCCTTGCTCGCCGCCGACCCGTTCATTCAGCTGACCCAAGCACTGGGCGCCCTCCCTTCGTTGCTGGGCTTCGCCGTGCCCGAGGAAGCCa                                                                                                                                            <  2:4206135/150‑1 (MQ=255)
          cgcTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCa                                                                                                                                            <  1:2963996/150‑1 (MQ=255)
               ggCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATc                                                                                                                                       >  2:2569778/1‑150 (MQ=255)
                  aTGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCcgc                                                                                                                                    >  2:1405062/1‑150 (MQ=255)
                      tACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCa                                                                                                                                >  2:3946881/1‑150 (MQ=255)
                      tACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCa                                                                                                                                >  1:3698636/1‑150 (MQ=255)
                          ggTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcgg                                                                                                                            <  2:3877511/150‑1 (MQ=255)
                          ggTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcgg                                                                                                                            <  1:2784471/150‑1 (MQ=255)
                           gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                           >  2:2055908/1‑150 (MQ=255)
                           gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                           >  1:2303970/1‑150 (MQ=255)
                           gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                           >  1:3920350/1‑150 (MQ=255)
                            tCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggcg                                                                                                                          >  1:3049559/1‑150 (MQ=255)
                              ttGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCCGCGGCGAc                                                                                                                        >  2:914168/1‑150 (MQ=255)
                              ttGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGAc                                                                                                                        >  2:3238713/1‑150 (MQ=255)
                                     cGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCAACAGCCTGc                                                                                                                 <  2:3149121/150‑1 (MQ=255)
                                       aCCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCa                                                                                                               >  2:3767734/1‑150 (MQ=255)
                                        cccGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAg                                                                                                              >  2:3474046/1‑150 (MQ=255)
                                         ccGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGc                                                                                                             >  1:2228403/1‑150 (MQ=255)
                                         ccGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGc                                                                                                             >  1:681354/1‑150 (MQ=255)
                                         ccGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCGGCCAGc                                                                                                             >  1:4155555/1‑150 (MQ=255)
                                                 aTGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACg                                                                                                     >  2:1870558/1‑150 (MQ=255)
                                                 aTGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACg                                                                                                     >  1:1494067/1‑150 (MQ=255)
                                                 aTGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCG‑CCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGt                                                                                                    <  1:896290/150‑1 (MQ=255)
                                                  tGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGt                                                                                                    >  1:3612159/1‑150 (MQ=255)
                                                         gTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGc                                                                                             <  1:67760/150‑1 (MQ=255)
                                                           cTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCa                                                                                           <  2:287886/150‑1 (MQ=255)
                                                           cTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCa                                                                                           <  2:148376/150‑1 (MQ=255)
                                                            tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAg                                                                                          >  1:261752/1‑150 (MQ=255)
                                                            tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAg                                                                                          >  2:4476264/1‑150 (MQ=255)
                                                            tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAg                                                                                          >  1:1770802/1‑150 (MQ=255)
                                                             gggAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAgg                                                                                         >  1:1803602/1‑150 (MQ=255)
                                                             gggAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAgg                                                                                         >  2:1431564/1‑150 (MQ=255)
                                                              ggAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGc                                                                                        >  2:2615090/1‑150 (MQ=255)
                                                               gAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCCTGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGcc                                                                                       >  1:646483/1‑150 (MQ=255)
                                                               gAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGcc                                                                                       >  2:3278721/1‑150 (MQ=255)
                                                               gAAACCACCTTGCTCGCCGCCGACCCGAACATTCACCTGACCCAAGCACAGGGCGCCCTCACGTCGTAGCTGGGCTACGCCGTGCCCGACGAAGCCAAGATACCCGACACCCCCGACAGACTGCCAGCACAGGACGTGGCCCGCCAGGcc                                                                                       >  1:2909230/1‑150 (MQ=255)
                                                                    caccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGc                                                                                  <  2:283021/150‑1 (MQ=255)
                                                                      ccgtGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACCGCCGGCCCGCCCGGGACGTGGCCCGCCCGGCCGGGGCCa                                                                                >  1:542808/4‑150 (MQ=255)
                                                                         tGCTCGCCGCCGACCCGTTCATGCAGCTGACCCCATCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGGCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgct                                                                             <  2:1764429/150‑1 (MQ=255)
                                                                             cgccgcCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgctggct                                                                         >  2:2462677/1‑150 (MQ=255)
                                                                                 gccgACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTc                                                                     >  2:338932/1‑150 (MQ=255)
                                                                                    gACCCGTTCATTCACCTGACCCAAGCACGGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGc                                                                  <  2:3650706/150‑1 (MQ=255)
                                                                                       ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTTGGCTTAGCCGTGACCGACGAAGCCATGATCCGCGCCAGCGGCGACACCCTGCCACCCCAGGACGTGGCCCCCCAGGCCAGGGCCAGCTGGCTGCTCGGCgtg                                                               >  1:410602/1‑150 (MQ=255)
                                                                                       ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                               >  1:1415178/1‑150 (MQ=255)
                                                                                       ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                               >  1:3500437/1‑150 (MQ=255)
                                                                                       ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                               >  1:3825167/1‑150 (MQ=255)
                                                                                       ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                               >  2:3926568/1‑150 (MQ=255)
                                                                                         gTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATTATACGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtggt                                                             >  1:800727/1‑150 (MQ=255)
                                                                                           tcattcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGc                                                           <  1:4415406/150‑1 (MQ=255)
                                                                                            cattcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCt                                                          <  1:147453/150‑1 (MQ=255)
                                                                                            cattcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCt                                                          <  2:2020020/150‑1 (MQ=255)
                                                                                            cattcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCt                                                          <  1:1234196/150‑1 (MQ=255)
                                                                                             attcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTg                                                         >  2:2694586/1‑150 (MQ=255)
                                                                                             attcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTg                                                         >  1:1363489/1‑150 (MQ=255)
                                                                                               tcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGt                                                       <  1:1641002/150‑1 (MQ=255)
                                                                                               tcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGt                                                       <  1:1069803/150‑1 (MQ=255)
                                                                                               tcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGt                                                       <  2:3368792/150‑1 (MQ=255)
                                                                                                    tGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCGGGTATACg                                                  >  1:2120492/1‑150 (MQ=255)
                                                                                                        ccAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCct                                              >  2:4394976/1‑150 (MQ=255)
                                                                                                        ccAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCct                                              >  2:119673/1‑150 (MQ=255)
                                                                                                         cAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCctg                                             >  2:3519913/1‑150 (MQ=255)
                                                                                                         cAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCAGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCctg                                             >  2:4316403/1‑150 (MQ=255)
                                                                                                          aaGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCctgc                                            >  1:215043/1‑150 (MQ=255)
                                                                                                                tGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGcccc                                      <  2:3109425/150‑1 (MQ=255)
                                                                                                                   gcgcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGcc                                   >  2:2831781/1‑150 (MQ=255)
                                                                                                                     gcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCCACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctg                                 >  1:1309553/1‑150 (MQ=255)
                                                                                                                       ccTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgct                               <  2:3049559/150‑1 (MQ=255)
                                                                                                                        cTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgctg                              >  2:2238052/1‑150 (MQ=255)
                                                                                                                        cTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgctg                              <  1:4236107/150‑1 (MQ=255)
                                                                                                                         tCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgctgc                             >  2:739888/1‑150 (MQ=255)
                                                                                                                           ccGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgctgctg                           >  2:1948631/1‑150 (MQ=255)
                                                                                                                               cgtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCg                       <  2:584088/150‑1 (MQ=255)
                                                                                                                                gtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg                      >  1:2715766/1‑150 (MQ=255)
                                                                                                                                gtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg                      <  2:283766/150‑1 (MQ=255)
                                                                                                                                gtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg                      <  1:4363912/150‑1 (MQ=255)
                                                                                                                                    cTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTa                  >  1:2199955/1‑150 (MQ=255)
                                                                                                                                     tgtgCTTCGCCGTGCCCGACGAAGCCATTATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGACATGCCTGTGCCAGCTGGCTGCTCGGCGTGTTGCTTGTATACTGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTAt                 <  1:3356539/147‑1 (MQ=255)
                                                                                                                                     tGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTAt                 >  1:1726864/1‑150 (MQ=255)
                                                                                                                                      tggCTTCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATg                >  2:1352633/2‑150 (MQ=255)
                                                                                                                                      gggCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATg                >  1:953512/1‑150 (MQ=255)
                                                                                                                                      gggCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATg                >  2:1728138/1‑150 (MQ=255)
                                                                                                                                      gggCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATg                >  1:464830/1‑150 (MQ=255)
                                                                                                                                      gggCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATg                >  1:3663968/1‑150 (MQ=255)
                                                                                                                                       ggCTTCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGc               >  2:2636821/1‑150 (MQ=255)
                                                                                                                                        gCTTCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGcc              >  1:3864224/1‑150 (MQ=255)
                                                                                                                                         cTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCt             >  2:410785/1‑150 (MQ=255)
                                                                                                                                          ttCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGCCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTCTGCCtg            >  2:4518186/1‑150 (MQ=255)
                                                                                                                                          ttCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGccct            <  2:1136312/150‑3 (MQ=255)
                                                                                                                                          ttCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCtg            <  2:1539022/150‑1 (MQ=255)
                                                                                                                                          ttCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCtg            >  1:3736385/1‑150 (MQ=255)
                                                                                                                                              ccGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGtggc        <  1:3133742/150‑1 (MQ=255)
                                                                                                                                                gTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGtggcgc      >  2:3613275/1‑150 (MQ=255)
                                                                                                                                                    ccGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGtggcgctggc  <  1:1157303/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGCTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGTGGCGCTGGC  >  NC_002947/5858613‑5858910

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑