Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,651,399 C→T 5.6% P80S (CCG→TCG)  PP_3217 → aliphatic sulfonate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,651,3990CT 5.6% 268.3 / 3.4 89P80S (CCG→TCG) PP_3217aliphatic sulfonate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (46/37);  new base T (3/2);  total (50/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AGAGTTACAACGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCGGCCATCTGGAAA  >  NC_002947/3651250‑3651548
                                                                                                                                                     |                                                                                                                                                     
agagTTACAACGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACgcgc                                                                                                                                                       >  2:1640371/1‑150 (MQ=255)
    ttACAACGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGcgccgc                                                                                                                                                    >  2:472803/1‑149 (MQ=255)
    ttACAACGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGcgccgc                                                                                                                                                    <  1:472803/149‑1 (MQ=255)
     tACAACGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAATAACTGCCTTACAGCCTGGAGTGGAATCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGccc                                                                                                                                                  <  2:455068/150‑1 (MQ=255)
      aCAACGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGcccc                                                                                                                                                 >  2:1493237/1‑150 (MQ=255)
        aaCGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTg                                                                                                                                               >  1:717887/1‑150 (MQ=255)
        aaCGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTg                                                                                                                                               >  2:1819061/1‑150 (MQ=255)
         aCGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTCGCCGAACCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGt                                                                                                                                              <  1:331276/150‑1 (MQ=255)
          cGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGtt                                                                                                                                             >  1:1242233/1‑150 (MQ=255)
          cGCTCGTGCCGTGATGGAAGCCGCCGGCGGCCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGGCGCGTCGcccagcag                                                                                                                                             <  1:1211034/150‑6 (MQ=255)
               gTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCg                                                                                                                                        >  1:470030/1‑150 (MQ=255)
                 gCCATGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCACCCCTGTTCCTCGGt                                                                                                                                      >  2:1677180/1‑150 (MQ=255)
                    gTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGcc                                                                                                                                   >  2:1593551/1‑150 (MQ=255)
                    gTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGcc                                                                                                                                   >  2:2138307/1‑150 (MQ=255)
                       aTGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTg                                                                                                                                <  1:996208/150‑1 (MQ=255)
                         ggAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTggg                                                                                                                              >  2:254300/1‑150 (MQ=255)
                           aaGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTggggg                                                                                                                            <  2:2357669/150‑1 (MQ=255)
                                     cGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCaa                                                                                                                  <  1:1362079/150‑1 (MQ=255)
                                      gTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGTTCacg                                                                                                                 <  2:1242233/150‑3 (MQ=255)
                                      gTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAg                                                                                                                 <  1:2054728/150‑1 (MQ=255)
                                       tGCTAGAACAACTACCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAgg                                                                                                                <  2:2244095/150‑1 (MQ=255)
                                          tAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGa                                                                                                                      >  1:2302156/1‑141 (MQ=255)
                                          tAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGa                                                                                                                      <  2:2302156/141‑1 (MQ=255)
                                           aGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGAt                                                                                                            >  1:1220540/1‑150 (MQ=255)
                                            gAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCCCCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCt                                                                                                                                          <  1:532313/119‑1 (MQ=255)
                                            gAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCt                                                                                                                                          >  2:532313/1‑119 (MQ=255)
                                               caacTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCGCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGCGCGCCGATCAAGGTGATTGCg                                                                                                        <  2:1026945/150‑1 (MQ=255)
                                                  cTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTc                                                                                                     >  2:2144102/1‑150 (MQ=255)
                                                   tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTgggg                                                                                                                             >  2:1059292/1‑125 (MQ=255)
                                                   tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTgggg                                                                                                                             <  1:1059292/125‑1 (MQ=255)
                                                   tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTg                                                                                                                                <  1:1027242/122‑1 (MQ=255)
                                                   tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTg                                                                                                                                >  2:1027242/1‑122 (MQ=255)
                                                   tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAg                                                                                                                 >  2:980370/1‑137 (MQ=255)
                                                   tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAg                                                                                                                 <  1:980370/137‑1 (MQ=255)
                                                   tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCa                                                                                                    <  1:1493237/150‑1 (MQ=255)
                                                     ccTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGcgccg                                                                                                                                                     <  1:179329/99‑1 (MQ=255)
                                                     ccTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGcgccg                                                                                                                                                     >  2:179329/1‑99 (MQ=255)
                                                      cTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGcc                                                                                                 >  1:398559/1‑150 (MQ=255)
                                                           acCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAg                                                                                            >  1:1454512/1‑150 (MQ=255)
                                                             ccTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGaa                                                                                          >  2:1178390/1‑150 (MQ=255)
                                                                ggAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCt                                                                                       <  1:1634839/150‑1 (MQ=255)
                                                                   gTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGa                                                                                    <  2:1224935/150‑1 (MQ=255)
                                                                    tGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGaa                                                                                   >  1:110231/1‑150 (MQ=255)
                                                                     ggAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAg                                                                                  <  2:1454512/150‑1 (MQ=255)
                                                                      gAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAgg                                                                                 >  1:1633488/1‑150 (MQ=255)
                                                                      gAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAgg                                                                                 >  2:1300856/1‑150 (MQ=255)
                                                                       aaGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAg                                                                                                                 >  2:776719/1‑117 (MQ=255)
                                                                       aaGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAg                                                                                                                 <  1:776719/117‑1 (MQ=255)
                                                                         gCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGt                                                                              <  1:1300856/150‑1 (MQ=255)
                                                                          cAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTg                                                                             <  1:229924/150‑1 (MQ=255)
                                                                           aGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTgg                                                                            <  1:1178390/150‑1 (MQ=255)
                                                                           aGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTgg                                                                            >  2:947608/1‑150 (MQ=255)
                                                                           aGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTgg                                                                            <  2:1448060/150‑1 (MQ=255)
                                                                             ttCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGcc                                                                          >  2:460819/1‑150 (MQ=255)
                                                                             ttCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTCGcc                                                                          >  1:1214406/1‑150 (MQ=255)
                                                                               cACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCAt                                                                        >  1:1448158/1‑150 (MQ=255)
                                                                                aCCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAAACTCGAAGGCGTGGCCATc                                                                       >  1:1767896/1‑150 (MQ=255)
                                                                                aCCGCAGGCTCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATcc                                                                      <  1:1733746/150‑1 (MQ=255)
                                                                                 ccGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCg                                                                                                                                        >  1:1489402/1‑84 (MQ=255)
                                                                                 ccGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCg                                                                                                                                        <  2:1489402/84‑1 (MQ=255)
                                                                                           cccGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCaa                                                            <  1:1819061/150‑1 (MQ=255)
                                                                                            ccGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAg                                                           <  1:252874/150‑1 (MQ=255)
                                                                                                tGGCCGAAGCGCTCAACGTCGGCAGATTCGATATAGGCTTGCTGGGCGTTCCATCTCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGGGTGGCCATCATGGCGCGCAAGGAtt                                                       >  1:909693/1‑150 (MQ=255)
                                                                                                     gAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCg                                                  >  2:2153111/1‑150 (MQ=255)
                                                                                                     gAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGAGACGAAAATCGTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCg                                                  >  2:593505/1‑150 (MQ=255)
                                                                                                      aaGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCg                                                                                                        >  1:1677267/1‑95 (MQ=255)
                                                                                                      aaGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCg                                                                                                        <  2:1677267/95‑1 (MQ=255)
                                                                                                      aaGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGa                                                 <  1:2138307/150‑1 (MQ=255)
                                                                                                        gcgcTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATc                                               <  2:909693/150‑1 (MQ=255)
                                                                                                         cgcTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGaa                                                                                   <  2:422390/113‑1 (MQ=255)
                                                                                                         cgcTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGaa                                                                                   >  1:422390/1‑113 (MQ=255)
                                                                                                           cTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCcac                                            >  2:2334951/1‑150 (MQ=255)
                                                                                                             cAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCa                                                                                             >  2:1711539/1‑99 (MQ=255)
                                                                                                             cAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCa                                                                                             <  1:1711539/99‑1 (MQ=255)
                                                                                                                  tCGGCAGTTTAGATATTGTCCCGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCCCCGATGAAGGTGCTTGCTGTTAGCCTTGAGAACCTCGAAGGTGTGGCCATCCTCGCGCGCAAGGATTCTCCGATCCACAGCCTGc                                     <  2:470030/150‑1 (MQ=255)
                                                                                                                    ggCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTg                                                                             >  2:753000/1‑108 (MQ=255)
                                                                                                                    ggCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTg                                                                             <  1:753000/108‑1 (MQ=255)
                                                                                                                     gCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGACTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAgg                                  >  1:139609/1‑150 (MQ=255)
                                                                                                                           tcGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTgc                            >  1:1072385/1‑150 (MQ=255)
                                                                                                                             gATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTgcgc                          >  2:1814607/1‑150 (MQ=255)
                                                                                                                                    gCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGaa                                                                                   <  2:1388475/86‑1 (MQ=255)
                                                                                                                                    gCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGaa                                                                                   >  1:1388475/1‑86 (MQ=255)
                                                                                                                                    gCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAAc                   >  2:1547273/1‑150 (MQ=255)
                                                                                                                                      tcaaggggCGGCGCGTCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAAcgc                 >  2:590374/6‑150 (MQ=255)
                                                                                                                                      ctgctgGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCa                                                                                                    <  2:738420/67‑1 (MQ=255)
                                                                                                                                      ctgctgGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCa                                                                                                    >  1:738420/1‑67 (MQ=255)
                                                                                                                                       tgctgGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAAcgcg                >  2:1645124/1‑150 (MQ=255)
                                                                                                                                       tgctgGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGCTCAGCCGGCAGAACCTCGAAGGCGTGGCCGTCCTTGCGCGCAGTGCTGCTCCGAGCCACAGCCTGCAGTACCTTCGCGGCAAAcgcg                >  2:508360/1‑150 (MQ=255)
                                                                                                                                       acaggggCGGCGCGTCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGc                       >  1:1383421/5‑143 (MQ=255)
                                                                                                                                       acaggggCGGCGCGTCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGc                       <  2:1383421/139‑1 (MQ=255)
                                                                                                                                          tgGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCgg             >  2:1176501/1‑150 (MQ=255)
                                                                                                                                          tgGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCgg             >  2:1651383/1‑150 (MQ=255)
                                                                                                                                           gggCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCGGc            <  2:1220540/150‑1 (MQ=255)
                                                                                                                                              cGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCGGCCAt         >  1:2029533/1‑150 (MQ=255)
                                                                                                                                                     ccgccCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCGGCCATCTGGaaa  <  2:110231/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGAGTTACAACGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAGGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCGGCCATCTGGAAA  >  NC_002947/3651250‑3651548

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