Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,813,291:1 +A 5.3% coding (325/852 nt) PP_3369 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,813,2911.A 5.3% 472.3 / 18.4 114coding (325/852 nt)PP_3369LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (50/58);  new base A (2/4);  total (52/62)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 119 total aligned reads displayed.

ACGAAGGCCAACGTGTTGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCAACTCTGGTGGTCAGCCTGCTG  >  NC_002947/3813145‑3813440
                                                                                                                                                   |                                                                                                                                                      
aCGAAGGCCAACGTGTTGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGgc                                                                                                                                                     >  2:487579/1‑150 (MQ=255)
aCGAAGGCCAACGTGTTGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGgc                                                                                                                                                     >  1:885098/1‑150 (MQ=255)
 cGAAGGCCAACGTGTTGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGgcg                                                                                                                                                    >  1:2764182/1‑150 (MQ=255)
     ggCCAACGTGTTGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTg                                                                                                                                                <  2:1552418/150‑1 (MQ=255)
      gCCAACGTGTTGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTgg                                                                                                                                               >  1:2879696/1‑150 (MQ=255)
      gCCAACGTGTTGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTgg                                                                                                                                               <  1:3035201/150‑1 (MQ=255)
        cAACGTGTTGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGc                                                                                                                                             >  2:2829772/1‑150 (MQ=255)
          aCGTGTTGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCAGGCGCCGGGCGa                                                                                                                                           <  2:1352294/150‑1 (MQ=255)
           cGTGTTGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAg                                                                                                                                          <  2:1172737/150‑1 (MQ=255)
            ttgtTGCGGAACTCGCCGAGCAGATGGCCATGTACTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGc                                                                                                                                         <  1:97092/149‑1 (MQ=255)
             tcttGCCGAACTCGCCGAGCAGATTGCCATGTCTTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTTGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCAGGCGCCGGGCGatga                                                                                                                                        <  1:27657/148‑4 (MQ=255)
               ttGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGt                                                                                                                                      <  1:487579/150‑1 (MQ=255)
                tGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTa                                                                                                                                     >  2:2920403/1‑150 (MQ=255)
                   ggAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCTAAGTGGTGCGTTCGGCAGCGCGGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCa                                                                                                                                  <  1:1637312/150‑1 (MQ=255)
                   ggAACTCACCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCa                                                                                                                                  <  1:3048855/150‑1 (MQ=255)
                    gAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAg                                                                                                                                 >  1:2896926/1‑150 (MQ=255)
                     aaCTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAgc                                                                                                                                <  2:2450293/150‑1 (MQ=255)
                        tCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCt                                                                                                                             <  1:1132060/150‑1 (MQ=255)
                          gCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTAc                                                                                                                           >  1:1027218/1‑150 (MQ=255)
                          gCCGAGCAGATGGCCATGTCCTCTTTTGACCTGGAGCATTTCGCCCGAGCCGGCCAACAGAAAGTGGAAGGCGAAGTGGAGCTCTCTTCACCACCGTCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTAc                                                                                                                           >  2:347851/1‑150 (MQ=255)
                           ccGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTAcc                                                                                                                          <  1:1512727/150‑1 (MQ=255)
                                 aGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGc                                                                                                                    >  2:1162611/1‑150 (MQ=255)
                                    tGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTgc                                                                                                                 <  1:1063442/150‑1 (MQ=255)
                                    tGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTgc                                                                                                                 <  2:2171022/150‑1 (MQ=255)
                                         aTGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCAGGCGCCGGGCGatga                                                                                                                                        >  2:2815753/1‑119 (MQ=255)
                                           gTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGaa                                                                                                          <  2:2896926/150‑1 (MQ=255)
                                            tCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAg                                                                                                         >  2:14774/1‑150 (MQ=255)
                                            tCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAg                                                                                                         >  1:312390/1‑150 (MQ=255)
                                              ctctTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCt                                                                                                       >  2:2643672/1‑150 (MQ=255)
                                               tctTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTg                                                                                                      <  1:1417023/150‑1 (MQ=255)
                                                  tttGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCt                                                                                                                                        <  2:311158/113‑1 (MQ=255)
                                                  tttGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCt                                                                                                                                        >  1:311158/1‑113 (MQ=255)
                                                  tttGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTa                                                                                                   <  1:1673436/150‑1 (MQ=255)
                                                     gCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGgc                                                                                                <  1:1834202/150‑1 (MQ=255)
                                                     acccTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGgc                                                                                                >  2:2109155/2‑150 (MQ=255)
                                                      cccTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGgca                                                                                               >  2:1502267/1‑150 (MQ=255)
                                                       ccTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCt                                                                                                                                        <  2:2151437/108‑1 (MQ=255)
                                                       ccTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCt                                                                                                                                        >  1:2151437/1‑108 (MQ=255)
                                                       ccTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTg                                                                                                      <  1:2210111/142‑1 (MQ=255)
                                                       ccTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTg                                                                                                      >  2:2210111/1‑142 (MQ=255)
                                                       ccTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGgcag                                                                                              >  1:1350755/1‑150 (MQ=255)
                                                            aGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTACTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAAt                                                                                         <  2:2262742/150‑1 (MQ=255)
                                                              cATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATcc                                                                                      >  2:578267/1‑150 (MQ=255)
                                                                tttCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGTGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCaa                                                                                    <  2:387782/150‑1 (MQ=255)
                                                                 ttCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAg                                                                                   <  1:541544/150‑1 (MQ=255)
                                                                 ttCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAg                                                                                   <  2:2504523/150‑1 (MQ=255)
                                                                        gcgcCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCt                                                                            <  1:578267/150‑1 (MQ=255)
                                                                           ccGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTgg                                                                         <  2:145121/150‑1 (MQ=255)
                                                                           ccGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTgg                                                                         <  2:2556507/150‑1 (MQ=255)
                                                                            cGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGt                                                                        >  1:1931508/1‑150 (MQ=255)
                                                                                  acagacagTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTccgccg                                                                  <  2:2879696/150‑1 (MQ=255)
                                                                                    agacagTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTa                                                                                                   >  2:3057331/1‑116 (MQ=255)
                                                                                    agacagTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTa                                                                                                   <  1:3057331/116‑1 (MQ=255)
                                                                                    agacagTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCcgcg                                                                >  1:564329/1‑150 (MQ=255)
                                                                                        agTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGc                                                            <  1:14774/150‑1 (MQ=255)
                                                                                         gTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCt                                                           <  1:2829772/150‑1 (MQ=255)
                                                                                              aGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACAt                                                      >  1:1659269/1‑150 (MQ=255)
                                                                                              aGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTACGCATCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACAt                                                      <  1:1597998/150‑1 (MQ=255)
                                                                                               ggCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATc                                                     <  1:2142144/150‑1 (MQ=255)
                                                                                                gCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGTAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCa                                                    <  1:2920403/150‑1 (MQ=255)
                                                                                                 cGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAg                                                   >  2:536348/1‑150 (MQ=255)
                                                                                                  gAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCAGGCGCCGGGCGa                                                                                                                                           <  1:1447615/62‑1 (MQ=38)
                                                                                                  gAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCAGGCGCCGGGCGa                                                                                                                                           >  2:1447615/1‑62 (MQ=38)
                                                                                                  gAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGc                                                  <  2:1027692/150‑1 (MQ=255)
                                                                                                  gAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGc                                                  >  1:2159567/1‑150 (MQ=255)
                                                                                                   aaGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCa                                                 >  2:1926363/1‑150 (MQ=255)
                                                                                                        gtgCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCaa                                                                                    >  2:2954870/1‑110 (MQ=255)
                                                                                                          gCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCt                                          >  1:2599968/1‑150 (MQ=255)
                                                                                                           ctctCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAg                                                                                                                     >  1:1213559/1‑75 (MQ=255)
                                                                                                           ctctCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCt                                          <  2:87848/149‑1 (MQ=255)
                                                                                                           ctctCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCt                                          >  1:87848/1‑149 (MQ=255)
                                                                                                           ctctCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTg                                         <  1:2109155/150‑1 (MQ=255)
                                                                                                            tctcGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTgg                                        >  1:2800027/1‑150 (MQ=255)
                                                                                                             ctcGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCaa                                                                                    <  1:2848715/105‑1 (MQ=255)
                                                                                                             ctcGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCaa                                                                                    >  2:2848715/1‑105 (MQ=255)
                                                                                                             ctcGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGc                                       <  1:2193736/150‑1 (MQ=255)
                                                                                                              tcGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGcc                                      >  1:2724026/1‑150 (MQ=255)
                                                                                                              tcGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGAGCAAGCTGTACAAGCGCTACCACCGGATGCGCCTGTAGCTGGTAGGCGCTAAAGCCAAGGCAGCCTTTGTCCTCCGCGAAGCTGACATCAGCATCGCCCTGGcc                                      >  2:345686/1‑150 (MQ=255)
                                                                                                                 gCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAgg                                   >  2:11539/1‑150 (MQ=255)
                                                                                                                  caccacCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGc                                  <  1:1742223/150‑1 (MQ=255)
                                                                                                                  caccacCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGc                                  <  2:2764182/150‑1 (MQ=255)
                                                                                                                  caccacCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCAAAGCTGCGAATGACGACGGTAGGCAGCAAATCAAAGGCATGCTTGGTCCGCAGCGAAGATCACATCAGCATCGCCCTGGCCAGGc                                  >  2:187376/1‑150 (MQ=255)
                                                                                                                   accaccGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGcc                                 <  1:1162611/150‑1 (MQ=255)
                                                                                                                      accGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAACTGTAACATCGCTACCCACAGCTGCGCCTGAAGCTGGTATGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGAc                              >  1:1598257/1‑150 (MQ=255)
                                                                                                                       ccGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACa                             <  2:183343/150‑1 (MQ=255)
                                                                                                                         ggCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACaga                           >  2:2966608/1‑150 (MQ=255)
                                                                                                                             ctgctgGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCa                       >  1:557037/1‑150 (MQ=255)
                                                                                                                              tgctgGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCaa                      <  2:1931508/150‑1 (MQ=255)
                                                                                                                                   ggCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGCATCGCCTGCCGAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCAACTCtg                 <  2:885098/150‑1 (MQ=255)
                                                                                                                                   ggCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCAACTGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCAACTCtg                 <  1:187376/150‑1 (MQ=255)
                                                                                                                                    gcagcaTCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCTAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCAACTCtgg                >  1:657145/1‑150 (MQ=255)
                                                                                                                                    gcagcaTCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCAACTCtgg                <  2:607632/150‑1 (MQ=255)
                                                                                                                                       gcaTCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTccgcc                                                                   <  1:2555423/96‑1 (MQ=255)
                                                                                                                                       gcaTCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTccgcc                                                                   >  2:2555423/1‑96 (MQ=255)
                                                                                                                                            atcGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACaga                           <  2:1189933/131‑1 (MQ=255)
                                                                                                                                            atcGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACaga                           >  1:1189933/1‑131 (MQ=255)
                                                                                                                                                cccGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCAACTCtggtg              <  1:866602/140‑1 (MQ=255)
                                                                                                                                                cccGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCAACTCtggtg              >  2:866602/1‑140 (MQ=255)
                                                                                                                                                 ccGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCAACTCTGGTGGTCAGCctgct   <  2:1714270/150‑1 (MQ=255)
                                                                                                                                                  cGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCAACTCTGGTGGTCAGCctgctg  <  1:347851/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ACGAAGGCCAACGTGTTGCGGAACTCGCCGAGCAGATGGCCATGTCCTCTTTTGCCCTGGAGCATTTCGCCAGCGCCGGCCAACAGACAGTGGAAGGCGAAGTGGTGCTCTCGGCACCACCGGCGCTGCTGGGCAGCATCATCGCCCGGCGCCTGGGCGAGCTGTACCAGCGCTACCCACAGCTGCGCCTGAAGCTGGTAGGCAGCAAATCCAAGGCATCCTTGGTCCGCCGCGAAGCTGACATCAGCATCGCCCTGGCCAGGCCGACAGAGCCAACTCTGGTGGTCAGCCTGCTG  >  NC_002947/3813145‑3813440

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