Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,028,664 C→G 6.3% A295G (GCC→GGC)  PP_4432 → Xaa‑Pro aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,028,6640CG 6.3% 284.6 / 4.8 93A295G (GCC→GGC) PP_4432Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (42/44);  new base G (2/4);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGA  >  NC_002947/5028515‑5028813
                                                                                                                                                     |                                                                                                                                                      
cGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGgccgc                                                                                                                                                        <  1:1681105/150‑1 (MQ=255)
  gcgcCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGGCa                                                                                                                                                      <  1:2062520/150‑1 (MQ=255)
  gcgcCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGGCa                                                                                                                                                      <  2:2267333/150‑1 (MQ=255)
   cgcCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGCCCTCCGGGCCGCCaa                                                                                                                                                     <  2:1422300/150‑1 (MQ=255)
       ggCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCt                                                                                                                                                 >  2:196415/1‑150 (MQ=255)
       ggCACCTTCTTCGAAATCACCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCt                                                                                                                                                 <  1:241396/150‑1 (MQ=255)
             ttcttcGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCCGACCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAac                                                                                                                                           >  1:1507479/1‑150 (MQ=255)
               cttcGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAacac                                                                                                                                         <  1:1314462/150‑1 (MQ=255)
                ttcGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAAcacc                                                                                                                                        >  1:2104168/1‑150 (MQ=255)
                   gAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACt                                                                                                                                     >  1:854286/1‑150 (MQ=255)
                     aaTCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTgg                                                                                                                                   <  2:1575474/150‑1 (MQ=255)
                     aaTCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTgg                                                                                                                                   <  2:1183443/150‑1 (MQ=255)
                      aTCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGt                                                                                                                                  <  2:403276/150‑1 (MQ=255)
                      aTCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGt                                                                                                                                  <  2:1923312/150‑1 (MQ=255)
                      aTCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGACAACACCACTGGt                                                                                                                                  <  2:1264089/150‑1 (MQ=255)
                       tCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTg                                                                                                                                 >  1:1518834/1‑150 (MQ=255)
                       tCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTg                                                                                                                                 <  2:2362174/150‑1 (MQ=255)
                                    cAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGc                                                                                                                    <  1:196415/150‑1 (MQ=255)
                                       gcgcTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGc                                                                                                                                              >  2:1241045/1‑121 (MQ=255)
                                       gcgcTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGc                                                                                                                                              <  1:1241045/121‑1 (MQ=255)
                                          cTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCcttctt                                                                                                              >  2:891751/1‑150 (MQ=255)
                                          cTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCcttctt                                                                                                              >  1:1513780/1‑150 (MQ=255)
                                           tATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGAACGTCCcctctgc                                                                                                             >  2:1046503/1‑148 (MQ=255)
                                             tCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCgg                                                                                                           >  2:1587909/1‑150 (MQ=255)
                                                   cccGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCt                                                                                                     >  1:1519329/1‑150 (MQ=255)
                                                      gCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGc                                                                                                                                              <  2:52971/106‑1 (MQ=255)
                                                      gCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGc                                                                                                                                              >  1:52971/1‑106 (MQ=255)
                                                      gCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTga                                                                                                  >  2:2045987/1‑150 (MQ=255)
                                                      gCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTga                                                                                                  <  2:1757020/150‑1 (MQ=255)
                                                       cTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTgag                                                                                                 >  1:431385/1‑150 (MQ=255)
                                                         gTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGaa                                                                                               <  2:854286/150‑1 (MQ=255)
                                                         gTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGaa                                                                                               <  2:2047073/150‑1 (MQ=255)
                                                         gTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGaa                                                                                               <  1:1179763/150‑1 (MQ=255)
                                                         gTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGaa                                                                                               <  1:2045987/150‑1 (MQ=255)
                                                          tCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAg                                                                                              >  2:122375/1‑150 (MQ=255)
                                                           cccGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGt                                                                                             >  2:1982142/1‑150 (MQ=255)
                                                            ccGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGtt                                                                                            <  2:1838887/150‑1 (MQ=255)
                                                            ccGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGtt                                                                                            >  2:1827806/1‑150 (MQ=255)
                                                             cGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTc                                                                                           >  1:2180917/1‑150 (MQ=255)
                                                              gAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCg                                                                                          >  2:784632/1‑150 (MQ=255)
                                                               aaCTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGa                                                                                         <  2:2391298/150‑1 (MQ=255)
                                                               aaCTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGa                                                                                         <  1:278910/150‑1 (MQ=255)
                                                                   gTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTg                                                                                                                                                <  1:1316348/91‑1 (MQ=255)
                                                                   gTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTg                                                                                                                                                >  2:1316348/1‑91 (MQ=255)
                                                                   gTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATc                                                                                     <  2:1518834/150‑1 (MQ=255)
                                                                      tACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCcac                                                                                  <  1:1827806/150‑1 (MQ=255)
                                                                        ccTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACaa                                                                                >  1:1824771/1‑150 (MQ=255)
                                                                             gCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACa                                                                           <  2:2104168/150‑1 (MQ=255)
                                                                              cAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGGGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAg                                                                          >  1:870931/1‑150 (MQ=255)
                                                                                 gCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTTCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCg                                                                       >  1:1266375/1‑150 (MQ=255)
                                                                                  ccAACGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGt                                                                      >  2:761839/1‑150 (MQ=255)
                                                                                    aCCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTg                                                                    <  1:630699/150‑1 (MQ=255)
                                                                                               tCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCaa                                                        >  2:2335291/1‑150 (MQ=255)
                                                                                                 aTCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATc                                                      >  1:281338/1‑150 (MQ=255)
                                                                                                 aTCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATc                                                      <  2:281338/150‑1 (MQ=255)
                                                                                                  tCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CCGTTACCCAATCg                                                     >  2:1022625/1‑150 (MQ=255)
                                                                                                  tCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGAAGTCGCACGTGCCTTCGTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCg                                                     >  1:1495304/1‑150 (MQ=255)
                                                                                                          gagaAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAAt                                                       >  1:1668406/1‑140 (MQ=255)
                                                                                                          gagaAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAAt                                                       <  2:1668406/140‑1 (MQ=255)
                                                                                                             aaaGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTAt                                          <  2:1824771/150‑1 (MQ=255)
                                                                                                                gCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATcca                                       >  1:1677664/1‑150 (MQ=255)
                                                                                                                   gtcgtcGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATccacca                                    <  2:2180917/150‑1 (MQ=255)
                                                                                                                      gtcgAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGAt                                 >  1:1419379/1‑150 (MQ=255)
                                                                                                                        cgAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTCGCAAAATCCCTGGTGAACCCTTACCGCGCGGCGTGGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTg                               <  1:891751/150‑1 (MQ=255)
                                                                                                                         gAAGGTATCGCTGCCGGCCTCGAGGCCGTCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTgg                              >  2:1375202/1‑150 (MQ=255)
                                                                                                                         gAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCcac                                                                                  <  2:2212507/99‑1 (MQ=255)
                                                                                                                         gAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCcac                                                                                  >  1:2212507/1‑99 (MQ=255)
                                                                                                                         gAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTAt                                          >  2:2229074/1‑138 (MQ=255)
                                                                                                                         gAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTAt                                          <  1:2229074/138‑1 (MQ=255)
                                                                                                                         gAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTgg                              <  2:1677664/150‑1 (MQ=255)
                                                                                                                          aaGGTATCGCTGCCGGCCTCGAGGCGGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGATTGAGAAGTTCGATAT‑TTCTGGAACATCCTCTGTCCGTTACCCAATCGGTATCAGCTATCCACCAGATTggg                             <  1:1022625/150‑1 (MQ=255)
                                                                                                                              tATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGa                         <  1:2335291/150‑1 (MQ=255)
                                                                                                                                tCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTggg                             >  2:1930480/1‑144 (MQ=255)
                                                                                                                                tCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTggg                             <  1:1930480/144‑1 (MQ=255)
                                                                                                                                tCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAgc                       >  1:94474/1‑150 (MQ=255)
                                                                                                                                 gcttggcGGCCTCGAGGCCGGCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGt                                                   <  1:1902059/115‑1 (MQ=255)
                                                                                                                                 gcttggcGGCCTCGAGGCCGGCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGt                                                   >  2:1902059/7‑121 (MQ=255)
                                                                                                                                   cTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCa                    <  2:2389618/150‑1 (MQ=255)
                                                                                                                                    tGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCAc                   <  1:1982142/150‑1 (MQ=255)
                                                                                                                                     ggcGGCCTCGAGGCCGGCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGa                         >  2:798729/3‑143 (MQ=255)
                                                                                                                                     ggcGGCCTCGAGGCCGGCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGa                         <  1:798729/141‑1 (MQ=255)
                                                                                                                                       cGGCCTCGAGGCCGCCAAGCCTGGCATCAAAGGTGGTGATATCGCACGTGCCTTCTTCGGCTTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCAt                <  2:1507479/150‑1 (MQ=255)
                                                                                                                                       cGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTgg                              >  1:410920/1‑136 (MQ=255)
                                                                                                                                       cGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTgg                              <  2:410920/136‑1 (MQ=255)
                                                                                                                                         gCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGa              <  1:1587909/150‑1 (MQ=255)
                                                                                                                                           cTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑Cg                                                                 <  2:1852795/97‑1 (MQ=255)
                                                                                                                                           cTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCg                                                     >  1:1522413/1‑109 (MQ=255)
                                                                                                                                           cTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCg                                                     <  2:1522413/109‑1 (MQ=255)
                                                                                                                                           cTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCAAAAGGACAGCCGTTG‑Cg                                                                 >  1:1852795/1‑97 (MQ=255)
                                                                                                                                            tCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGcc           >  1:2045185/1‑150 (MQ=255)
                                                                                                                                            tCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTAGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGcc           >  1:49845/1‑150 (MQ=255)
                                                                                                                                              gAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTc                                                                                           >  2:282630/1‑69 (MQ=255)
                                                                                                                                              gAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTc                                                                                           <  1:282630/69‑1 (MQ=255)
                                                                                                                                              gAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTg         >  2:1927168/1‑150 (MQ=255)
                                                                                                                                                   cgccAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTAGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGccc    <  2:49845/150‑1 (MQ=255)
                                                                                                                                                    gccAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTg         >  1:1311958/1‑144 (MQ=255)
                                                                                                                                                    gccAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTg         <  2:1311958/144‑1 (MQ=255)
                                                                                                                                                     ccAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGa  >  2:1210172/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTG‑CGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGA  >  NC_002947/5028515‑5028813

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