Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,270,512 C→A 20.3% S376S (TCC→TCA PP_4643 → xanthine/uracil permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,270,5120CA20.3% 403.2 / 89.7 172S376S (TCC→TCAPP_4643xanthine/uracil permease family protein
Reads supporting (aligned to +/- strand):  ref base C (55/82);  new base A (21/14);  total (76/96)
Fisher's exact test for biased strand distribution p-value = 3.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

GTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGAT  >  NC_002947/5270363‑5270661
                                                                                                                                                     |                                                                                                                                                     
gtggtgGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTcc                                                                                                                                                       <  1:1498972/150‑1 (MQ=255)
  ggtggCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCaa                                                                                                                                                     <  2:2076130/150‑1 (MQ=255)
  ggtggCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCaa                                                                                                                                                     <  1:2525458/150‑1 (MQ=255)
  ggtggCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCaaa                                                                                                                                                     <  2:149430/150‑1 (MQ=255)
  ggtggCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCaaa                                                                                                                                                     <  1:519012/150‑1 (MQ=255)
  ggtggCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCaaa                                                                                                                                                     <  1:2939756/150‑1 (MQ=255)
  ggtggCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCaaa                                                                                                                                                     <  1:768574/150‑1 (MQ=255)
     ggCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGt                                                                                                                                                  >  1:2106666/1‑150 (MQ=255)
      gcaaCGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTc                                                                                                                                                 <  1:4811402/147‑1 (MQ=255)
       cTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCa                                                                                                                                                <  2:1746752/150‑1 (MQ=255)
        tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGTGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGACAg                                                                                                                                               >  2:1291679/1‑150 (MQ=255)
        aaCGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTcga                                                                                                                                               >  2:3738023/2‑148 (MQ=255)
         aCGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGc                                                                                                                                              <  1:4398780/150‑1 (MQ=255)
             gcggcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACa                                                                                                                                          <  1:4694033/150‑1 (MQ=255)
                gcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGa                                                                                                                                       <  2:5066834/150‑1 (MQ=255)
                 cgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGat                                                                                                                                      <  2:4819094/150‑2 (MQ=255)
                  gTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaac                                                                                                                                     <  1:1686651/150‑1 (MQ=255)
                  gTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaac                                                                                                                                     >  2:5072919/1‑150 (MQ=255)
                     aTTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaacaac                                                                                                                                  <  2:2195285/150‑1 (MQ=255)
                         tggtggTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCa                                                                                                                              >  2:4272855/1‑150 (MQ=255)
                           gtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAc                                                                                                                            <  2:4380798/150‑1 (MQ=255)
                           gtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAc                                                                                                                            <  1:2405660/150‑1 (MQ=255)
                               tACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGa                                                                                                                        <  1:2894281/150‑1 (MQ=255)
                               tACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGa                                                                                                                        <  2:1804900/150‑1 (MQ=255)
                                aCTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCTTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGat                                                                                                                       >  2:1159488/1‑150 (MQ=255)
                                 cTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatc                                                                                                                      >  1:1667635/1‑150 (MQ=255)
                                 cTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatc                                                                                                                      >  1:4558170/1‑150 (MQ=255)
                                 cTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatc                                                                                                                      <  1:4059281/150‑1 (MQ=255)
                                  tCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatca                                                                                                                     >  2:4947319/1‑150 (MQ=255)
                                  tCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatca                                                                                                                     >  2:713279/1‑150 (MQ=255)
                                    ggCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatcatc                                                                                                                   >  1:4856913/1‑150 (MQ=255)
                                     gCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCAtcg                                                                                                                  >  2:4661375/1‑150 (MQ=255)
                                       ttGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCAtcgtc                                                                                                                <  2:2532517/150‑1 (MQ=255)
                                        ttctgcCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCAttgtgg                                                                                                               >  1:1344580/3‑148 (MQ=255)
                                         tctgcCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTcgc                                                                                                              <  2:2875540/149‑1 (MQ=255)
                                          ctgcCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTcgcg                                                                                                             <  1:3838102/150‑1 (MQ=255)
                                             ccAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGcc                                                                                                          >  1:3890207/1‑150 (MQ=255)
                                             ccAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGcc                                                                                                          <  2:729641/150‑1 (MQ=255)
                                                 tCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTCTCCGCACCCTGACCAAGGTCAGCTACAAGAACAACGTCAACCTCATCAGCGACGCGGCCTCGc                                                                                                      >  2:3171008/1‑150 (MQ=255)
                                                 tCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGc                                                                                                      <  1:3172046/150‑1 (MQ=255)
                                                    tGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTgg                                                                                                   <  2:199322/150‑1 (MQ=255)
                                                    tGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTgg                                                                                                   >  2:3818725/1‑150 (MQ=255)
                                                     gggCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTggg                                                                                                  <  2:3848154/150‑1 (MQ=255)
                                                              gATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCAt                                                                                         <  2:871578/150‑1 (MQ=255)
                                                              gATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCAt                                                                                         <  1:5073071/150‑1 (MQ=255)
                                                               aTTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATg                                                                                        <  1:3703536/150‑1 (MQ=255)
                                                                   ccgccgTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATcc                                                                                    >  1:1339580/1‑150 (MQ=255)
                                                                       cgTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGAt                                                                                <  1:3762789/150‑1 (MQ=255)
                                                                       cgTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGAt                                                                                <  1:1169313/150‑1 (MQ=255)
                                                                           ccTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGcc                                                                            >  1:3110636/1‑150 (MQ=255)
                                                                           ccTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGcc                                                                            <  2:4558170/150‑1 (MQ=255)
                                                                            cTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCg                                                                           <  1:3284562/150‑1 (MQ=255)
                                                                                agcagtgcTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCAcc                                                                       >  2:2276481/8‑150 (MQ=255)
                                                                                 ccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCg                                                                      >  1:1181783/1‑150 (MQ=255)
                                                                                    gtgcTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCAccagaa                                                                   >  2:716111/4‑146 (MQ=255)
                                                                                    gTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGAcc                                                                   <  2:698733/150‑1 (MQ=255)
                                                                                      gcTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCAccagaatt                                                                 <  2:1070928/149‑7 (MQ=255)
                                                                                        tGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCt                                                               <  1:2651249/150‑1 (MQ=255)
                                                                                        tGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCt                                                               <  1:3194471/150‑1 (MQ=255)
                                                                                         gggCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTa                                                              <  1:614301/150‑1 (MQ=255)
                                                                                           gcggcgCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTacc                                                            <  2:3575522/150‑1 (MQ=255)
                                                                                             ggcgCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTaccac                                                          <  1:4833517/150‑1 (MQ=255)
                                                                                               tgcCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTAccatca                                                        <  2:848889/149‑1 (MQ=255)
                                                                                                    gTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTcc                                                   <  2:2663985/150‑1 (MQ=255)
                                                                                                        cGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGaa                                               <  1:2624323/150‑1 (MQ=255)
                                                                                                         gTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAAc                                              >  2:4840785/1‑150 (MQ=255)
                                                                                                           gCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCggcctg                                            <  2:3131713/150‑6 (MQ=255)
                                                                                                             tGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGt                                          >  2:3188399/1‑150 (MQ=255)
                                                                                                             tGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGt                                          <  2:1700030/150‑1 (MQ=255)
                                                                                                             tGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGt                                          >  2:4488339/1‑150 (MQ=255)
                                                                                                             tGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGt                                          <  2:2064861/150‑1 (MQ=255)
                                                                                                             tGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGt                                          >  2:2976275/1‑150 (MQ=255)
                                                                                                               ttCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTc                                        >  1:934921/1‑150 (MQ=255)
                                                                                                               ttCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCAGAACTGGTTc                                        >  1:1633209/1‑150 (MQ=255)
                                                                                                                   gCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCgaga                                    <  2:4872465/150‑1 (MQ=255)
                                                                                                                   gCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCgaga                                    >  1:2285501/1‑150 (MQ=255)
                                                                                                                   gCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAAAGTCAAACTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCgaga                                    <  1:2744090/150‑1 (MQ=255)
                                                                                                                            cggcCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTcc                           <  1:4948307/150‑1 (MQ=255)
                                                                                                                            cggcCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTcc                           <  1:4948350/150‑1 (MQ=255)
                                                                                                                             ggcCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCa                          >  1:840835/1‑150 (MQ=255)
                                                                                                                             cgcaaGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCa                          <  1:207339/146‑1 (MQ=255)
                                                                                                                                cAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTc                       >  1:3173480/1‑150 (MQ=255)
                                                                                                                                 aGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTcc                      <  1:1291679/150‑1 (MQ=255)
                                                                                                                                 aGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTcc                      >  2:2254966/1‑150 (MQ=255)
                                                                                                                                  gCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTcag                     <  1:4237579/150‑3 (MQ=255)
                                                                                                                                  gCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTcag                     >  1:4054603/1‑148 (MQ=255)
                                                                                                                                     gTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCa                  <  2:1668128/150‑1 (MQ=255)
                                                                                                                                       aTCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATc                >  2:2140631/1‑150 (MQ=255)
                                                                                                                                       aTCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATc                >  2:2140708/1‑150 (MQ=255)
                                                                                                                                        tCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCa               <  1:1731121/150‑1 (MQ=255)
                                                                                                                                        tCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCa               <  1:4272855/150‑1 (MQ=255)
                                                                                                                                        tCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCa               <  2:181953/150‑1 (MQ=255)
                                                                                                                                         ccGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAg              <  2:935302/150‑1 (MQ=255)
                                                                                                                                          cGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGc             >  1:4007818/1‑150 (MQ=255)
                                                                                                                                          cGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGc             <  2:1648491/150‑1 (MQ=255)
                                                                                                                                            cACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTc           <  1:552755/150‑1 (MQ=255)
                                                                                                                                               ccTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGACCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTcggc        <  1:713279/150‑1 (MQ=255)
                                                                                                                                                     cAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGag  >  2:3028709/1‑149 (MQ=255)
                                                                                                                                                     cAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGAt  >  1:1903880/1‑150 (MQ=255)
                                                                                                                                                     cAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGAt  >  2:3569220/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGAT  >  NC_002947/5270363‑5270661

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: