Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 833,814 C→G 5.7% P163A (CCG→GCG)  PP_0717 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947833,8140CG 5.7% 650.4 / 21.1 194P163A (CCG→GCG) PP_0717membrane protein
Reads supporting (aligned to +/- strand):  ref base C (101/82);  new base G (8/3);  total (109/85)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 201 total aligned reads displayed.

ACCATGATCGCTGCTGCCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAG  >  NC_002947/833670‑833961
                                                                                                                                                |                                                                                                                                                   
aCCATGATCGCTGCTGCCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTg                                                                                                                                                >  2:2019364/1‑150 (MQ=255)
 ccATGATCGCTGCTGCCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTgc                                                                                                                                               <  2:3289044/150‑1 (MQ=255)
  cATGATCGCTGCTGCCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCAGCGCCGGTgag                                                                                                                                              >  2:4847246/1‑148 (MQ=255)
  cATGATCGCTGCTGCCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTgcc                                                                                                                                              >  1:386132/1‑150 (MQ=255)
  cATGATCGCTGCTGCCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTgcc                                                                                                                                              >  1:1996386/1‑150 (MQ=255)
   aTGATCGCTGCTGCCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTgccg                                                                                                                                             >  1:2744425/1‑150 (MQ=255)
    tGATCGCTGCTGCCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATCCGTTGTTCCCCGAACCTGCCCATGCGCCAGCACCGCCGGTgccgc                                                                                                                                            >  2:2398340/1‑150 (MQ=255)
     gATCGCTGCTGCCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCGGCACCGGCGGTgctggc                                                                                                                                           <  1:4019965/150‑5 (MQ=255)
     gATCGCTGCTGCCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCGGCACCGGCGGTgctggc                                                                                                                                           <  1:1440160/150‑5 (MQ=255)
             ctgcCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCCTACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTgg                                                                                                                                   >  2:250716/1‑150 (MQ=255)
              tgcCGGCACCCGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTggg                                                                                                                                  <  1:3103705/150‑1 (MQ=255)
                ccGGCACCAGGGACTTCACCCTGGCATTGTCGGGGGTCGAAGCGCTGGCCAACGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTa                                                                                                                                <  2:3869641/150‑1 (MQ=255)
                  ggCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGc                                                                                                                              >  2:3520783/1‑150 (MQ=255)
                    cACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGa                                                                                                                            >  1:1668390/1‑150 (MQ=255)
                    cACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGa                                                                                                                            <  2:2422492/150‑1 (MQ=255)
                       cAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTCGCGAGCa                                                                                                                         >  2:5029071/1‑150 (MQ=255)
                         gTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATg                                                                                                                       >  2:346700/1‑150 (MQ=255)
                             cTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAg                                                                                                                   >  2:4859481/1‑150 (MQ=255)
                                cACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTg                                                                                                                >  2:1969999/1‑150 (MQ=255)
                                   ccTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGt                                                                                                             <  1:642464/150‑1 (MQ=255)
                                     tGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTgg                                                                                                           >  1:964350/1‑150 (MQ=255)
                                         aTTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCt                                                                                                       <  1:3831631/150‑1 (MQ=255)
                                                gtggtCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGcc                                                                                                >  2:2830808/1‑150 (MQ=255)
                                                 tggtCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCACTGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCa                                                                                               <  2:3631208/150‑1 (MQ=255)
                                                     cGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGtt                                                                                           <  2:2820003/150‑1 (MQ=255)
                                                         gcgcTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATc                                                                                       <  2:1031712/150‑1 (MQ=255)
                                                         gcgcTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATc                                                                                       <  1:1013571/150‑1 (MQ=255)
                                                           gcTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGt                                                                                     >  2:359347/1‑150 (MQ=255)
                                                           gcTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACATTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGt                                                                                     <  2:3838839/150‑1 (MQ=255)
                                                            cTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTg                                                                                    >  2:4148386/1‑150 (MQ=255)
                                                             tGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGa                                                                                   >  1:2877933/1‑150 (MQ=255)
                                                             tGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGa                                                                                   >  1:2668898/1‑150 (MQ=255)
                                                                ccAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGc                                                                                >  1:1034868/1‑150 (MQ=255)
                                                                 cAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGcc                                                                               >  1:294150/1‑150 (MQ=255)
                                                                 cAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGcc                                                                               <  2:476228/150‑1 (MQ=255)
                                                                 cAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGcc                                                                               >  1:43951/1‑150 (MQ=255)
                                                                 cAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGcc                                                                               >  1:3032418/1‑150 (MQ=255)
                                                                 cAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCACTGCGTGCCACGTTGATCGTGATGcc                                                                               >  1:4719952/1‑150 (MQ=255)
                                                                  aaaGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCg                                                                              >  1:393936/1‑150 (MQ=255)
                                                                    aGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGc                                                                            <  2:4345581/150‑1 (MQ=255)
                                                                       cATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCtt                                                                         <  1:374399/150‑1 (MQ=255)
                                                                        aTGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCttt                                                                        >  2:1378720/1‑150 (MQ=255)
                                                                               tgGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTgg                                                                 <  1:2067095/150‑1 (MQ=255)
                                                                               tgGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTgg                                                                 <  2:249798/150‑1 (MQ=255)
                                                                               tgGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTgg                                                                 >  1:3449264/1‑150 (MQ=255)
                                                                               tgGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTgg                                                                 >  2:4084086/1‑150 (MQ=255)
                                                                               tgGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTgg                                                                 <  2:10870/150‑1 (MQ=255)
                                                                               tgGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTgg                                                                 <  1:1052499/150‑1 (MQ=255)
                                                                               tgGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTgg                                                                 <  2:543276/150‑1 (MQ=255)
                                                                                ggCAGTGCTGAGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGc                                                                >  1:4779570/1‑150 (MQ=255)
                                                                                  cAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCAt                                                              <  2:1892916/150‑1 (MQ=255)
                                                                                     tGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGa                                                           >  1:1207061/1‑150 (MQ=255)
                                                                                       cTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCAGTGATc                                                         >  2:668696/1‑150 (MQ=255)
                                                                                         gggCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGc                                                       >  1:2915292/1‑150 (MQ=255)
                                                                                           gCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGccc                                                     >  2:2878395/1‑150 (MQ=255)
                                                                                           gCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGccc                                                     <  1:4593956/150‑1 (MQ=255)
                                                                                           gCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGccc                                                     <  1:4350110/150‑1 (MQ=255)
                                                                                            cACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGcccc                                                    <  1:1964722/150‑1 (MQ=255)
                                                                                            cACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGcccc                                                    <  1:3891964/150‑1 (MQ=255)
                                                                                               gTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGa                                                 <  2:3583409/150‑1 (MQ=255)
                                                                                                   gtggtgCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAg                                             >  2:389626/1‑150 (MQ=255)
                                                                                                      gtgCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTAc                                          <  2:2038697/150‑1 (MQ=255)
                                                                                                           tACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGcc                                     >  1:118244/1‑150 (MQ=255)
                                                                                                           tACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGcc                                     >  2:4002861/1‑150 (MQ=255)
                                                                                                           tACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGcc                                     >  1:1154012/1‑150 (MQ=255)
                                                                                                              gttgttCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGtt                                  <  1:2242360/150‑1 (MQ=255)
                                                                                                               ttgttCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTg                                 >  1:3963141/1‑150 (MQ=255)
                                                                                                                     cccGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATg                           >  2:821319/1‑150 (MQ=255)
                                                                                                                     cccGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATg                           >  2:4556527/1‑150 (MQ=255)
                                                                                                                     cccGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATg                           >  1:4165092/1‑150 (MQ=255)
                                                                                                                     cccGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATg                           <  2:3486941/150‑1 (MQ=255)
                                                                                                                       cGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTTCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGaa                         <  1:1120531/150‑1 (MQ=255)
                                                                                                                       cGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGaa                         <  2:277524/150‑1 (MQ=255)
                                                                                                                       cGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGaa                         <  1:174604/150‑1 (MQ=255)
                                                                                                                        gAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAg                        <  1:250464/150‑1 (MQ=255)
                                                                                                                        gAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAg                        <  2:43951/150‑1 (MQ=255)
                                                                                                                             tGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTgg                                                                 >  1:721869/1‑104 (MQ=255)
                                                                                                                             tGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTgg                                                                 <  2:721869/104‑1 (MQ=255)
                                                                                                                             tGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGt                   <  2:2208413/150‑1 (MQ=255)
                                                                                                                             tGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGt                   <  1:4748220/150‑1 (MQ=255)
                                                                                                                               ccAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGa                 <  1:5062928/150‑1 (MQ=255)
                                                                                                                               ccAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGa                 >  2:489263/1‑150 (MQ=255)
                                                                                                                                 aacgcgcCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAgt               >  1:3547065/4‑150 (MQ=255)
                                                                                                                                 aaTGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAgt               >  1:1279883/1‑150 (MQ=255)
                                                                                                                                 aaTGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATGGATCGCCCCTGAGCAGTACATGCCGTTGATCCTGAAGGCCGGGAgg               >  2:2290598/1‑149 (MQ=255)
                                                                                                                                  aTGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAgtg              >  1:719842/1‑150 (MQ=255)
                                                                                                                                    gcgcCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGtgtg            <  1:1485172/150‑1 (MQ=255)
                                                                                                                                      ggcggCACCGGCGGTGCTGGCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTggg          >  1:3408659/5‑150 (MQ=255)
                                                                                                                                      gcCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTggg          >  1:4254545/1‑150 (MQ=255)
                                                                                                                                      gcCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTggg          >  2:1747930/1‑150 (MQ=255)
                                                                                                                                      gcCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTggg          <  1:781052/150‑1 (MQ=255)
                                                                                                                                      gcCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTggg          <  1:4105181/150‑1 (MQ=255)
                                                                                                                                      gcCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTggg          <  1:3211195/150‑1 (MQ=255)
                                                                                                                                       ccAGCAGCGGCGGCACCGGCGGTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCAGTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGGGTgggg         >  1:1302673/1‑150 (MQ=255)
                                                                                                                                        cagcagCGGCGGCACCGGCGGTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGc        >  1:61326/1‑150 (MQ=255)
                                                                                                                                         agcagcGGCGGCACCGGCGGTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGca       >  1:3639002/1‑150 (MQ=255)
                                                                                                                                          gCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGcaa      >  1:2175074/1‑150 (MQ=255)
                                                                                                                                           cACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGcaac     <  2:1628126/150‑1 (MQ=255)
                                                                                                                                             ccgccgGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTctg                                                                     <  1:4962927/84‑1 (MQ=255)
                                                                                                                                             ccgccgGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTctg                                                                     >  2:4962927/1‑84 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ACCATGATCGCTGCTGCCGGCACCAGTGACTTCACCCTGGCATTGTCGGTGGTCGAAGCGCTGGCCAAAGGCATGCTGCTGGCAGTGCTGGGCACGTCGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAG  >  NC_002947/833670‑833961

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