Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 827,358 G→T 7.2% G182G (GGG→GGT)  ycaC‑I → putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947827,3580GT 7.2% 454.3 / 18.6 139G182G (GGG→GGT) ycaC‑Iputative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (62/67);  new base T (4/6);  total (66/73)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 142 total aligned reads displayed.

GTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCAT  >  NC_002947/827209‑827507
                                                                                                                                                     |                                                                                                                                                     
gTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGGCTGGCGCAACGACATCGAAggg                                                                                                                                                       >  1:3819093/1‑150 (MQ=255)
   tttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTg                                                                                                                                                    >  1:4616580/4‑150 (MQ=2)
   ttCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTg                                                                                                                                                    <  2:3025898/150‑1 (MQ=255)
    ttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTgg                                                                                                                                                   <  2:177167/148‑1 (MQ=2)
    ttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTgg                                                                                                                                                   <  2:4762051/148‑1 (MQ=2)
    tcgtcgTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTgg                                                                                                                                                   <  2:2622799/150‑1 (MQ=255)
      gtcgtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGt                                                                                                                                                 >  2:1575548/1‑150 (MQ=255)
       tcgtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTa                                                                                                                                                >  2:1945089/1‑150 (MQ=255)
       tcgtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTa                                                                                                                                                >  1:2302036/1‑150 (MQ=255)
       tcgtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTa                                                                                                                                                >  1:396995/1‑150 (MQ=255)
         gtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTAcc                                                                                                                                              <  1:4426018/150‑1 (MQ=255)
          tcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTAccc                                                                                                                                             <  1:3227054/150‑1 (MQ=255)
          tcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGAGACTGGCGCAACGACATCGAAGGGCTGGGTAccc                                                                                                                                             >  2:4849537/1‑150 (MQ=255)
            aCCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTg                                                                                                                                           <  1:4939395/150‑1 (MQ=255)
              cGACGCCTCTGGCACGTTCAACGAGCTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGtt                                                                                                                                         <  2:293613/150‑1 (MQ=255)
              cGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGtt                                                                                                                                         <  2:772851/150‑1 (MQ=255)
              cGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGtt                                                                                                                                         <  2:3747650/150‑1 (MQ=255)
              cGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGtt                                                                                                                                         <  2:772831/150‑1 (MQ=255)
                aCGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTtct                                                                                                                                       <  2:4020464/150‑1 (MQ=255)
                  gCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTcc                                                                                                                                     <  2:4271239/150‑1 (MQ=255)
                    ctctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCaa                                                                                                                                   <  1:1886555/150‑1 (MQ=255)
                     tctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAc                                                                                                                                  <  1:122094/150‑1 (MQ=255)
                     tctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAc                                                                                                                                  <  2:3611939/150‑1 (MQ=255)
                     tctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACAACATCGAAGGGCTGGGTACCCTGTTCTCCAAc                                                                                                                                  >  1:5053191/1‑150 (MQ=255)
                      ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc                                                                                                                                 <  1:2719344/150‑1 (MQ=255)
                      ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc                                                                                                                                 <  1:2719414/150‑1 (MQ=255)
                      ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc                                                                                                                                 <  1:4852602/150‑1 (MQ=255)
                      ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc                                                                                                                                 <  1:2719304/150‑1 (MQ=255)
                               tCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCgg                                                                                                                        >  1:3792496/1‑150 (MQ=2)
                                cAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCGGa                                                                                                                       <  1:1764462/150‑1 (MQ=2)
                                cAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGa                                                                                                                       >  1:869204/1‑150 (MQ=255)
                                 aatgAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCGGat                                                                                                                      <  1:2996529/147‑2 (MQ=2)
                                 aaCGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGAc                                                                                                                      >  1:1341281/1‑150 (MQ=255)
                                 aaCGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGAc                                                                                                                      >  1:531433/1‑150 (MQ=255)
                                    gAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTAc                                                                                                                   >  2:3385532/1‑150 (MQ=255)
                                         gACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTACACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCACaa                                                                                                              >  1:545994/1‑150 (MQ=255)
                                            ccGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCt                                                                                                           >  1:1350917/1‑150 (MQ=255)
                                            ccGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCt                                                                                                           >  2:4145747/1‑150 (MQ=255)
                                              gTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtga                                                                                                         <  2:577744/150‑1 (MQ=255)
                                              gTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtga                                                                                                         <  2:3827102/150‑1 (MQ=255)
                                               tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                                                        <  2:5013918/150‑1 (MQ=255)
                                               tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                                                        >  2:4294491/1‑150 (MQ=255)
                                               tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                                                        >  2:1840083/1‑150 (MQ=255)
                                               tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                                                        <  2:1473991/150‑1 (MQ=255)
                                               tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                                                        <  2:1474045/150‑1 (MQ=255)
                                               tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                                                        >  2:1375630/1‑150 (MQ=255)
                                               tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                                                        <  1:1839208/150‑1 (MQ=255)
                                                   tCGGCCTGGCGACGCATGGAGGCTGCCGGGGCACAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGAcc                                                                                                    >  2:4269834/1‑150 (MQ=255)
                                                    cGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCa                                                                                                   <  1:4034069/150‑1 (MQ=255)
                                                     ggCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAg                                                                                                  <  2:2603396/150‑1 (MQ=255)
                                                        cTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTa                                                                                               <  2:2374112/150‑1 (MQ=255)
                                                          ggCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTaca                                                                                             >  1:1842447/1‑150 (MQ=255)
                                                            cGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTacaac                                                                                           >  1:5040482/1‑150 (MQ=255)
                                                             gACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTacaaca                                                                                          >  2:4266242/1‑150 (MQ=255)
                                                                gCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGc                                                                                       >  2:2480328/1‑150 (MQ=255)
                                                                gCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATAGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGc                                                                                       >  2:4304505/1‑150 (MQ=255)
                                                                 cATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCt                                                                                      <  2:4805599/150‑1 (MQ=255)
                                                                  aTGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTg                                                                                     <  2:1889806/150‑1 (MQ=255)
                                                                       ggCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCaa                                                                                <  1:1988901/150‑1 (MQ=255)
                                                                       ggCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCCCAACCTGATGACCACCTACAACAAGCTGGCCaa                                                                                <  2:4122377/150‑1 (MQ=255)
                                                                        gCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCaaa                                                                               <  1:2356112/150‑1 (MQ=255)
                                                                          tGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATa                                                                             <  1:4299568/150‑1 (MQ=255)
                                                                          tGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATa                                                                             >  1:4935944/1‑150 (MQ=255)
                                                                          tGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATa                                                                             <  2:1274536/150‑1 (MQ=255)
                                                                            gcGGGGCTCAACTGATGACCTGGTTCTGCTTGGACTGCGAACTGCACCGTGACTGGCGAAACGACATCGAAGGGATGTGTACCCTGTTCTCAAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGt                                                                           <  1:1390189/149‑1 (MQ=21)
                                                                                    cAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTg                                                                   >  2:4202326/1‑150 (MQ=255)
                                                                                     aaCTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTgg                                                                  >  1:2032024/1‑150 (MQ=255)
                                                                                     aaCTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTgg                                                                  >  2:1663223/1‑150 (MQ=255)
                                                                                      aCTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTggg                                                                 <  1:1139431/150‑1 (MQ=255)
                                                                                      aCTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTggg                                                                 <  2:2122097/150‑1 (MQ=255)
                                                                                       cTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTgggg                                                                <  1:4515825/150‑1 (MQ=255)
                                                                                        tgatgaCCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTggggg                                                               >  1:701835/1‑150 (MQ=255)
                                                                                         gatgaCCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGc                                                              <  2:2059347/150‑1 (MQ=255)
                                                                                              ccTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCAc                                                         <  2:2243626/150‑1 (MQ=255)
                                                                                               cTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACg                                                        >  2:1249645/1‑150 (MQ=255)
                                                                                                  gTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACa                                                     <  1:1685665/150‑1 (MQ=255)
                                                                                                  gTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACa                                                     >  2:432210/1‑150 (MQ=255)
                                                                                                       gCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGcccc                                                >  2:4505902/1‑150 (MQ=255)
                                                                                                        cGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGccccc                                               >  2:1013886/1‑150 (MQ=255)
                                                                                                          tGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCcac                                             >  2:4929935/1‑150 (MQ=255)
                                                                                                           ggCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCcaga                                            <  1:3144714/150‑3 (MQ=255)
                                                                                                              cTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCa                                         >  1:3746518/1‑150 (MQ=255)
                                                                                                                  gAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCt                                     >  2:2079086/1‑150 (MQ=255)
                                                                                                                       gCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCaag                                >  2:3875309/1‑150 (MQ=255)
                                                                                                                         aCCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCaagcg                              >  1:4476104/1‑150 (MQ=255)
                                                                                                                         aCCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCaagcg                              >  1:1073534/1‑150 (MQ=255)
                                                                                                                         aCCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCaagcg                              >  1:2975963/1‑150 (MQ=255)
                                                                                                                         aCCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCaagcg                              >  2:3024100/1‑150 (MQ=255)
                                                                                                                         aCCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTACTGGGGGCTGCACGACAGCCCCCACATCAATCTGCaagcg                              >  2:3593845/1‑150 (MQ=255)
                                                                                                                              gACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCaagcgaagcg                         <  2:5105484/150‑1 (MQ=255)
                                                                                                                                   gcgcAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGTCAGCCCCCACATCAATCTGCAAGCGAAGCGGCaca                    >  2:3122609/1‑148 (MQ=255)
                                                                                                                                   gcgcAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCaaa                    <  2:3926793/150‑1 (MQ=255)
                                                                                                                                     gcAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCACAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACt                  <  2:545994/150‑1 (MQ=255)
                                                                                                                                       aaCGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTgg                <  2:1433119/150‑1 (MQ=255)
                                                                                                                                        aCGACATCCAACCGCTCGGAACCCTGTTCACCAACCATATCCCGCACAACCCCAACCACAAGACCACCAACAACAACCTGGCCAAATACTCCCTACTGGCGGCAGCACGACAGCCCCCACATCAATCTCCAAGCGAACCCCCAAACTGGc               <  2:4836230/150‑1 (MQ=255)
                                                                                                                                         cGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCa              <  1:4363798/150‑1 (MQ=255)
                                                                                                                                              tCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAAt         >  2:2110607/1‑150 (MQ=255)
                                                                                                                                                gAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCg       >  1:1815783/1‑150 (MQ=255)
                                                                                                                                                 aaGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGa      >  1:4407527/1‑150 (MQ=255)
                                                                                                                                                   gggCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGAcc    <  1:208469/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCAT  >  NC_002947/827209‑827507

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: ‑