Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 827,349 C→T 6.8% D179D (GAC→GAT)  ycaC‑I → putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947827,3490CT 6.8% 479.0 / 15.5 147D179D (GAC→GAT) ycaC‑Iputative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (63/74);  new base T (4/6);  total (67/80)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 151 total aligned reads displayed.

GGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCA  >  NC_002947/827200‑827495
                                                                                                                                                     |                                                                                                                                                  
ggCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGAc                                                                                                                                                    <  1:1267305/150‑1 (MQ=255)
  tttCGATGTCTTTGTCGTCACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGatat                                                                                                                                                  <  1:58776/149‑1 (MQ=2)
  tttCGATGTCTTTGTCGTCACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGatat                                                                                                                                                  <  2:3616879/149‑1 (MQ=2)
  cTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACAt                                                                                                                                                  >  2:3062638/1‑150 (MQ=255)
   ttCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAGCTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATc                                                                                                                                                 >  1:293613/1‑150 (MQ=255)
   ttCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATc                                                                                                                                                 <  2:1245524/150‑1 (MQ=255)
      gAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGaa                                                                                                                                              <  1:4059641/150‑1 (MQ=255)
      gAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGaa                                                                                                                                              <  1:2220625/150‑1 (MQ=255)
       aGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAg                                                                                                                                             <  1:785362/150‑1 (MQ=255)
         gTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGGCTGGCGCAACGACATCGAAggg                                                                                                                                           >  1:3819093/1‑150 (MQ=255)
            ttCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTg                                                                                                                                        <  2:3025898/150‑1 (MQ=255)
             ttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTgg                                                                                                                                       <  2:177167/148‑1 (MQ=2)
             ttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTgg                                                                                                                                       <  2:4916040/148‑1 (MQ=2)
             ttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTgg                                                                                                                                       <  2:4762051/148‑1 (MQ=2)
             tcgtcgTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTgg                                                                                                                                       <  1:3961710/150‑1 (MQ=255)
               gtcgtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGt                                                                                                                                     >  2:1575548/1‑150 (MQ=255)
                  gtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTAcc                                                                                                                                  <  1:4426018/150‑1 (MQ=255)
                   tcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTAccc                                                                                                                                 >  2:391437/1‑150 (MQ=255)
                     aCCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTg                                                                                                                               <  1:4939395/150‑1 (MQ=255)
                      ccGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGt                                                                                                                              <  1:498594/150‑1 (MQ=255)
                       cGACGCCTCTGGCACGTTCAACGAGCTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGtt                                                                                                                             <  2:293613/150‑1 (MQ=255)
                       cGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGtt                                                                                                                             <  2:772851/150‑1 (MQ=255)
                       cGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGtt                                                                                                                             <  2:2576872/150‑1 (MQ=255)
                       cGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGtt                                                                                                                             <  2:3747650/150‑1 (MQ=255)
                         aCGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTtct                                                                                                                           <  2:4020464/150‑1 (MQ=255)
                           gCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTcc                                                                                                                         <  2:4271239/150‑1 (MQ=255)
                           gCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTcc                                                                                                                         <  1:2233488/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:2719414/150‑1 (MQ=255)
                               ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc                                                                                                                     <  1:2719344/150‑1 (MQ=255)
                               ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc                                                                                                                     <  1:4852602/150‑1 (MQ=255)
                                  gCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCata                                                                                                                  >  1:4846493/1‑150 (MQ=255)
                                     ccttCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCacattc                                                                                                               >  2:1993294/3‑148 (MQ=2)
                                      gTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAg                                                                                                                                             <  2:4444431/119‑1 (MQ=255)
                                      gTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAg                                                                                                                                             >  1:4444431/1‑119 (MQ=255)
                                         cAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCGGa                                                                                                           <  2:4170473/150‑1 (MQ=2)
                                         cAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGa                                                                                                           >  1:869204/1‑150 (MQ=255)
                                          aaCGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGAc                                                                                                          >  1:5105800/1‑150 (MQ=255)
                                          aaCGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGAc                                                                                                          >  1:1341281/1‑150 (MQ=255)
                                           aCGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACt                                                                                                         >  2:208469/1‑150 (MQ=255)
                                             gAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTAc                                                                                                       >  2:3385532/1‑150 (MQ=255)
                                              aaCTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTAcc                                                                                                      <  1:4501505/150‑1 (MQ=255)
                                               aCTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAg                                                                                                                                             >  1:1117567/1‑110 (MQ=255)
                                                gtGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCGGATTATCGc                                                                                                    >  2:5117779/2‑150 (MQ=2)
                                                  gACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTACACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCACaa                                                                                                  >  1:545994/1‑150 (MQ=255)
                                                   aCCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAAc                                                                                                 >  2:406948/1‑150 (MQ=255)
                                                    cccGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAAcc                                                                                                >  2:4536301/1‑150 (MQ=255)
                                                     ccGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCt                                                                                               >  2:4145747/1‑150 (MQ=255)
                                                     ccGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCt                                                                                               >  1:1350917/1‑150 (MQ=255)
                                                       gTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtga                                                                                             <  2:577744/150‑1 (MQ=255)
                                                        tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                                            >  2:4294491/1‑150 (MQ=255)
                                                        tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                                            <  1:1765973/150‑1 (MQ=255)
                                                        tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                                            <  2:1474045/150‑1 (MQ=255)
                                                        tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                                            >  2:1840083/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)
                                                                      gACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTacaaca                                                                              >  2:4266242/1‑150 (MQ=255)
                                                                           aTGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTg                                                                         <  2:1889806/150‑1 (MQ=255)
                                                                                ggCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCaa                                                                    <  1:1988901/150‑1 (MQ=255)
                                                                                 gCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCaaa                                                                   <  1:2356112/150‑1 (MQ=255)
                                                                                 gCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCaaa                                                                   <  1:2356584/150‑1 (MQ=255)
                                                                                   tGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATa                                                                 <  1:4299568/150‑1 (MQ=255)
                                                                                     gcGGGGCTCAACTGATGACCTGGTTCTGCTTGGACTGCGAACTGCACCGTGACTGGCGAAACGACATCGAAGGGATGTGTACCCTGTTCTCAAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGt                                                               <  1:1390189/149‑1 (MQ=21)
                                                                                      cGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTg                                                              <  2:137385/150‑1 (MQ=255)
                                                                                              aaCTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTgg                                                      >  1:2032024/1‑150 (MQ=255)
                                                                                              aaCTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTgg                                                      >  2:1663223/1‑150 (MQ=255)
                                                                                               aCTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTggg                                                     <  2:2122097/150‑1 (MQ=255)
                                                                                                cTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTgggg                                                    <  1:4515825/150‑1 (MQ=255)
                                                                                                 tgatgaCCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTggggg                                                   >  1:701835/1‑150 (MQ=255)
                                                                                                 tgatgaCCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTggggg                                                   >  1:3596790/1‑150 (MQ=255)
                                                                                                  gatgaCCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGc                                                  <  2:2059347/150‑1 (MQ=255)
                                                                                                   atgaCCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCt                                                 <  1:3028781/150‑1 (MQ=255)
                                                                                                       ccTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCAc                                             <  2:2243626/150‑1 (MQ=255)
                                                                                                        cTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACg                                            >  2:1249645/1‑150 (MQ=255)
                                                                                                           gTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACa                                         >  2:432210/1‑150 (MQ=255)
                                                                                                           gTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACa                                         <  1:1685665/150‑1 (MQ=255)
                                                                                                           gTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACa                                         <  2:481266/150‑1 (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)
                                                                                                                                gCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCaag                    >  2:3875309/1‑150 (MQ=255)
                                                                                                                                  aCCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCaagcg                  >  1:2975963/1‑150 (MQ=255)
                                                                                                                                  aCCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCaagcg                  >  1:1073534/1‑150 (MQ=255)
                                                                                                                                       gACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCaagcgaagcg             <  2:5105484/150‑1 (MQ=255)
                                                                                                                                          tGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCa          >  2:3734898/1‑150 (MQ=255)
                                                                                                                                          tGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCa          >  2:2800543/1‑150 (MQ=255)
                                                                                                                                            gcgcAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGTCAGCCCCCACATCAATCTGCAAGCGAAGCGGCaca        >  2:3122609/1‑148 (MQ=255)
                                                                                                                                            gcgcAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCaaa        <  2:617496/150‑1 (MQ=255)
                                                                                                                                                aaCGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTgg    <  2:1433210/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)
                                                                                                                                                     |                                                                                                                                                  
GGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCA  >  NC_002947/827200‑827495

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