Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 152,158 A→T 5.2% intergenic (‑52/+42) PP_0143 ← / ← PP_0144 membrane protein/insulinase family metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947152,1580AT 5.2% 713.0 / 19.6 212intergenic (‑52/+42)PP_0143/PP_0144membrane protein/insulinase family metalloprotease
Reads supporting (aligned to +/- strand):  ref base A (99/99);  new base T (4/7);  total (104/108)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 215 total aligned reads displayed.

CTGCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTAGAG  >  NC_002947/152009‑152307
                                                                                                                                                     |                                                                                                                                                     
cTGCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCaa                                                                                                                                                       <  2:825293/150‑1 (MQ=255)
  gCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGc                                                                                                                                                     >  2:2414316/1‑150 (MQ=255)
     gATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCcg                                                                                                                                                  >  2:2043467/1‑150 (MQ=255)
      aTACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTGGAGTAGTTCGATGAGGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCcgc                                                                                                                                                 <  2:1714478/150‑1 (MQ=255)
        aCCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                               <  1:1673426/150‑1 (MQ=255)
         ccGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              <  1:2335885/150‑1 (MQ=255)
         ccGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                              <  1:2196026/150‑1 (MQ=255)
          cGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGaaa                                                                                                                                             <  2:4114041/150‑1 (MQ=255)
           ggAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGaaga                                                                                                                                            <  1:1565034/150‑3 (MQ=255)
           ggAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGaaga                                                                                                                                            <  1:2421299/150‑3 (MQ=255)
           ggAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAg                                                                                                                                            >  1:2869423/1‑150 (MQ=255)
            gagaGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAgg                                                                                                                                           >  1:2735231/1‑150 (MQ=255)
              gagGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGt                                                                                                                                         <  2:1277143/150‑1 (MQ=255)
              gagGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGt                                                                                                                                         <  1:3703516/150‑1 (MQ=255)
              gagGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGt                                                                                                                                         <  1:680302/150‑1 (MQ=255)
               agGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTc                                                                                                                                        >  2:3518865/1‑150 (MQ=255)
                ggAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTTGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGaagaggctg                                                                                                                                       <  1:3988601/150‑8 (MQ=255)
                ggAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGaagaggctg                                                                                                                                       <  2:2154156/150‑8 (MQ=255)
                ggAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCg                                                                                                                                       >  2:1633468/1‑150 (MQ=255)
                ggAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCg                                                                                                                                       <  1:1225345/150‑1 (MQ=255)
                 gAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCgg                                                                                                                                      >  2:3686443/1‑150 (MQ=255)
                 gAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCgg                                                                                                                                      <  1:3711787/150‑1 (MQ=255)
                 gAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCgg                                                                                                                                      >  2:2710604/1‑150 (MQ=255)
                   aaCCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGtc                                                                                                                                    >  2:4187095/1‑149 (MQ=255)
                     ccgccgAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACt                                                                                                                                  <  1:1179030/150‑1 (MQ=255)
                        ccgAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGc                                                                                                                               >  1:3256577/1‑150 (MQ=255)
                          gAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCaa                                                                                                                             >  1:1193390/1‑150 (MQ=255)
                            gcgcAATGATGCCCGGGTAGCGATTGAGTAGTTCGATGATGCGCTGGCCGATCTGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAATCCCGCGAAAGGTTCGGTACTGGCaaaa                                                                                                                           <  1:368218/150‑1 (MQ=255)
                              gcAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAgg                                                                                                                         >  1:1385966/1‑150 (MQ=255)
                               cAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGa                                                                                                                        <  2:1416764/150‑1 (MQ=255)
                                aaTGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGAt                                                                                                                       >  1:3215976/1‑150 (MQ=255)
                                aaTGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGAt                                                                                                                       >  1:2040231/1‑150 (MQ=255)
                                aaTGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGAt                                                                                                                       >  2:3981925/1‑150 (MQ=255)
                                aaTGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGAt                                                                                                                       >  1:1006507/1‑150 (MQ=255)
                                aaTGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGAt                                                                                                                       <  1:4260840/150‑1 (MQ=255)
                                  tgatgCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATAc                                                                                                                     <  2:679060/150‑1 (MQ=255)
                                      gCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGTGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCCCGGGCAAGACCGCGAAACGGTCGGTACTGGCAAAATGATaacacc                                                                                                                 >  1:3768156/1‑145 (MQ=255)
                                      gCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATAccgcc                                                                                                                 <  2:1306839/150‑1 (MQ=255)
                                       cccGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATAccgccg                                                                                                                <  2:12294/150‑1 (MQ=255)
                                                 gcTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCt                                                                                                      >  2:1828446/1‑150 (MQ=255)
                                                   ttGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGc                                                                                                    <  2:2788494/150‑1 (MQ=255)
                                                    tGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGcc                                                                                                   <  2:4260514/150‑1 (MQ=255)
                                                       gtagtTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCgtg                                                                                                >  2:4427206/1‑150 (MQ=255)
                                                        tagtTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCgtgt                                                                                               <  2:4256607/150‑1 (MQ=255)
                                                         agtTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTa                                                                                              <  1:4249582/150‑1 (MQ=255)
                                                         agtTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTa                                                                                              <  1:2772067/150‑1 (MQ=255)
                                                           ttCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGc                                                                                            >  1:2534675/1‑150 (MQ=255)
                                                           ttCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGc                                                                                            >  2:543747/1‑150 (MQ=255)
                                                           ttCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATAACGCCGATTTACCGCTGCCGTGTAGc                                                                                            >  1:127897/1‑150 (MQ=255)
                                                              gatgatGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGcacc                                                                                         <  1:601291/150‑1 (MQ=255)
                                                              gatgatGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGcacc                                                                                         >  2:29666/1‑150 (MQ=255)
                                                              gatgatGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGcacc                                                                                         >  1:3464602/1‑150 (MQ=255)
                                                              gatgatGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGcacc                                                                                         >  2:1993312/1‑150 (MQ=255)
                                                                tgatgCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGcaccac                                                                                       >  1:3881532/1‑150 (MQ=255)
                                                                 gatgCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACt                                                                                      <  1:517427/150‑1 (MQ=255)
                                                                  atgCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTg                                                                                     <  2:1845823/150‑1 (MQ=255)
                                                                    gcgcTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTggg                                                                                   <  2:3082050/150‑1 (MQ=255)
                                                                          gACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTc                                                                             >  2:587604/1‑150 (MQ=255)
                                                                             gATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAgg                                                                          >  2:2528109/1‑150 (MQ=255)
                                                                              atatGGGGGAGTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGt                                                                         <  2:2486029/150‑1 (MQ=255)
                                                                               tatGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTc                                                                        >  1:3204561/1‑150 (MQ=255)
                                                                                   ggggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCCGGCAAACCCCCGAAAGGGTCGGTACTGGCAAAACGCTACCGCCGATTTACCCCAGCCGGGTCGCACCCCTGGGCCGTGAAGGTCCaag                                                                    >  1:1255284/1‑148 (MQ=255)
                                                                                    gggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAGGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATgc                                                                   >  2:3776360/1‑150 (MQ=255)
                                                                                           ggggCATTGGGGGCTACGAATGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGTAAGGGTCGGTACTGTCAACAGGATTCGGCCGATTTACCGCTGCCGTGTAGCAACAGTGGGCCGTTCAGGTCCATgccgccgcc                                                            <  2:2103363/150‑1 (MQ=255)
                                                                                               cATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAgcg                                                        <  1:3795998/150‑1 (MQ=255)
                                                                                                   gggggCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAAc                                                    <  1:2414316/150‑1 (MQ=255)
                                                                                                   gggggCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAAc                                                    <  2:2356343/150‑1 (MQ=255)
                                                                                                    ggggCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACa                                                   >  1:3343634/1‑150 (MQ=255)
                                                                                                    ggggCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACa                                                   >  1:3939881/1‑150 (MQ=255)
                                                                                                         tACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGcccc                                              >  2:4081835/1‑150 (MQ=255)
                                                                                                            gACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCccgcc                                           >  2:2930363/1‑150 (MQ=255)
                                                                                                            gACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCccgcc                                           >  2:2930395/1‑150 (MQ=255)
                                                                                                               aGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCgccgcc                                        >  1:1843799/1‑150 (MQ=255)
                                                                                                                   aaGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGc                                    <  1:3399560/150‑1 (MQ=255)
                                                                                                                     gTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGccc                                  <  2:2735231/150‑1 (MQ=255)
                                                                                                                      tGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCa                                 <  2:3595543/150‑1 (MQ=255)
                                                                                                                      tGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGGTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCa                                 >  1:2731980/1‑150 (MQ=255)
                                                                                                                       gAGGGGGCTGTGCCAGCCTCTTCGCGGGCACGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAg                                <  1:2855352/150‑1 (MQ=255)
                                                                                                                        aGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGc                               <  1:1901817/150‑1 (MQ=255)
                                                                                                                         gctGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGcg                              >  2:3539241/3‑149 (MQ=255)
                                                                                                                          gTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGccc                             >  1:1796439/1‑150 (MQ=255)
                                                                                                                          gTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGccc                             >  2:4005904/1‑150 (MQ=255)
                                                                                                                             gCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCAc                          >  2:3044207/1‑150 (MQ=255)
                                                                                                                               tgtgCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGa                        <  2:3275007/150‑1 (MQ=255)
                                                                                                                               tgtgCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGa                        <  2:1553848/150‑1 (MQ=255)
                                                                                                                                   ccAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCaa                    >  1:3125398/1‑150 (MQ=255)
                                                                                                                                   ccAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCaa                    <  2:1127263/150‑1 (MQ=255)
                                                                                                                                    cAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAg                   <  2:3423434/150‑1 (MQ=255)
                                                                                                                                    cAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAg                   <  2:485734/150‑1 (MQ=255)
                                                                                                                                     aGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAgg                  <  1:2710604/150‑1 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGc                 >  2:1918429/1‑150 (MQ=255)
                                                                                                                                        cctttCGCGGGCTTGCCCGCGAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCAt              >  1:4109236/4‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCAt              <  2:109160/150‑1 (MQ=255)
                                                                                                                                          tttCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATc             >  1:98891/2‑150 (MQ=255)
                                                                                                                                               cgGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAAccc        <  2:617579/150‑1 (MQ=255)
                                                                                                                                                gggCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGTCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCg       >  1:2316374/1‑150 (MQ=255)
                                                                                                                                                  gCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTa     >  1:4442445/1‑150 (MQ=255)
                                                                                                                                                  gCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTa     <  1:2153668/150‑1 (MQ=255)
                                                                                                                                                     aGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTagag  >  1:1094123/1‑150 (MQ=255)
                                                                                                                                                     aGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGACGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGACCAAGGCCATCAAccctaacag  >  1:4498282/1‑144 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTAGAG  >  NC_002947/152009‑152307

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: