Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 683,826 T→C 7.1% S10G (AGC→GGC)  PP_0585 ← MerR family metal‑sensitive transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947683,8260TC 7.1% 243.1 / 6.6 84S10G (AGC→GGC) PP_0585MerR family metal‑sensitive transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (35/43);  new base C (4/2);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

AAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCA  >  NC_002947/683679‑683971
                                                                                                                                                   |                                                                                                                                                 
aaaaGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGGTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTgc                                                                                                                                                 >  2:819110/1‑150 (MQ=255)
    gCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCCgcagcg                                                                                                                                             >  1:515713/1‑150 (MQ=255)
     cccAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGcggcgg                                                                                                                                            <  1:490359/150‑1 (MQ=255)
       cAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCggcgggc                                                                                                                                          >  1:817519/1‑150 (MQ=255)
        aGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCggcgggcg                                                                                                                                         >  2:62247/1‑150 (MQ=255)
         ggTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCCGCagcgggctc                                                                                                                                        <  2:515713/150‑3 (MQ=255)
          gTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGggcggc                                                                                                                                       >  2:1425486/1‑150 (MQ=255)
                     gcTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGt                                                                                                                            <  1:1133102/150‑1 (MQ=255)
                     gcTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGt                                                                                                                            >  2:285523/1‑150 (MQ=255)
                         gatgaAGGCCAGGTTGTGCAAGCCTTGCGCCGGGTACAGGCGGTAGCCGCTGGCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCAt                                                                                                                        <  1:200464/150‑1 (MQ=255)
                           tgaAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATgg                                                                                                                      >  1:1389926/1‑150 (MQ=255)
                           tgaAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATgg                                                                                                                      <  1:653605/150‑1 (MQ=255)
                                 ccAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGcctc                                                                                                                >  1:288694/1‑150 (MQ=255)
                                  cAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGcctcc                                                                                                               >  2:1325231/1‑150 (MQ=255)
                                    ggTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTg                                                                                                             <  1:41807/150‑1 (MQ=255)
                                     gTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTgcgg                                                                                                                                               <  1:1389782/115‑1 (MQ=255)
                                     gTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTgcgg                                                                                                                                               >  2:1389782/1‑115 (MQ=255)
                                     gTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTgg                                                                                                            <  1:644579/150‑1 (MQ=255)
                                     gTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCAgg                                                                                                            <  1:1425486/150‑1 (MQ=255)
                                            aaGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGtgctgc                                                                                                     >  1:502800/1‑150 (MQ=255)
                                                 ttGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCtggtg                                                                                                >  1:1166355/1‑150 (MQ=255)
                                                 ttGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCtggtg                                                                                                <  2:160953/150‑1 (MQ=255)
                                                 ttGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCtggtg                                                                                                >  1:317201/1‑150 (MQ=255)
                                                  tGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCtggtgg                                                                                               <  2:502800/150‑1 (MQ=255)
                                                  tGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCtggtgg                                                                                               <  2:821548/150‑1 (MQ=255)
                                                  tGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCtggtgg                                                                                               <  1:1400402/150‑1 (MQ=255)
                                                      ccTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGcggcgg                                                                                                                                            <  2:478252/101‑1 (MQ=255)
                                                      ccTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGcggcgg                                                                                                                                            >  1:478252/1‑101 (MQ=255)
                                                      ccTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTggg                                                                                              >  2:262435/1‑147 (MQ=255)
                                                      ccTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTggg                                                                                              <  1:262435/147‑1 (MQ=255)
                                                      ccTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGttt                                                                                           <  1:825979/150‑1 (MQ=255)
                                                       cTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTGTGAGCAGGCCGATGGATTCGGAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGGCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTc                                                                                          >  2:53830/1‑150 (MQ=255)
                                                          gTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGGGCCCCCTGGGTGTTGCTGGTAGGCTTCCAc                                                                                       >  2:605648/1‑150 (MQ=255)
                                                             cAGGCGGTAGCCGCTGTCGCTGCGAGGGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTTCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGtt                                                                                    <  2:169669/150‑1 (MQ=255)
                                                                  ggTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAg                                                                               <  2:1389926/150‑1 (MQ=255)
                                                                   gTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGc                                                                              <  2:1166355/150‑1 (MQ=255)
                                                                       ccGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCcaca                                                                          <  2:817519/150‑1 (MQ=255)
                                                                        cGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAg                                                                         <  1:1407611/150‑1 (MQ=255)
                                                                          cTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTg                                                                       >  1:869310/1‑150 (MQ=255)
                                                                           tGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGc                                                                      <  1:1401372/150‑1 (MQ=255)
                                                                                 tGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCt                                                                >  2:419083/1‑150 (MQ=255)
                                                                                   cGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTga                                                              <  2:869310/150‑1 (MQ=255)
                                                                                    gAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTgac                                                             >  1:1378396/1‑150 (MQ=255)
                                                                                         gCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGc                                                        <  1:125542/150‑1 (MQ=255)
                                                                                          cGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCt                                                       <  1:62247/150‑1 (MQ=255)
                                                                                              tttGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCAt                                                                                                                        <  2:442655/81‑1 (MQ=255)
                                                                                              tttGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCAt                                                                                                                        >  1:442655/1‑81 (MQ=255)
                                                                                                 gAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTg                                                <  1:192802/150‑1 (MQ=255)
                                                                                                       gCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGCCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGCGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTAACGCCGCTGACACTGGACAgg                                          >  1:1321325/1‑150 (MQ=255)
                                                                                                        ccGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAggg                                         <  1:605648/150‑1 (MQ=255)
                                                                                                        ccGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAggg                                         <  2:41071/150‑1 (MQ=255)
                                                                                                         cGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGc                                        >  1:858432/1‑150 (MQ=255)
                                                                                                             ggATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGt                                                                                             <  2:985133/93‑1 (MQ=255)
                                                                                                             ggATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGt                                                                                             >  1:985133/1‑93 (MQ=255)
                                                                                                              gATTCGTAGTAACGGATCATCTTGGTGCTGAGTCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATgg                                   <  1:419083/150‑1 (MQ=255)
                                                                                                              gATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAg                                           >  1:1339544/1‑142 (MQ=255)
                                                                                                              gATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAg                                           <  2:1339544/142‑1 (MQ=255)
                                                                                                                 tCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATggcgg                                <  1:181041/150‑1 (MQ=255)
                                                                                                                   gtagtaACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCt                              >  1:31846/1‑150 (MQ=255)
                                                                                                                   gtagtaACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCt                              >  2:378398/1‑150 (MQ=255)
                                                                                                                      gtaACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCt                           <  2:1378396/150‑1 (MQ=255)
                                                                                                                        aaCGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGc                         >  2:862960/1‑150 (MQ=255)
                                                                                                                          cGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTgg                                                                                                            <  1:1393729/65‑1 (MQ=255)
                                                                                                                          cGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTgg                                                                                                            >  2:1393729/1‑65 (MQ=255)
                                                                                                                           ggATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCACTGTGGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAg                      >  1:154690/1‑150 (MQ=255)
                                                                                                                            gATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCa                                                                                        <  2:868707/83‑1 (MQ=255)
                                                                                                                            gATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCa                                                                                        >  1:868707/1‑83 (MQ=255)
                                                                                                                                atcTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCg                 <  1:285523/150‑1 (MQ=255)
                                                                                                                                  cTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCggg               >  2:854059/1‑150 (MQ=255)
                                                                                                                                   ttGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGt              >  2:639341/1‑150 (MQ=255)
                                                                                                                                    tGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGtt             >  2:188910/1‑150 (MQ=255)
                                                                                                                                    tGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGtt             <  1:378398/150‑1 (MQ=255)
                                                                                                                                      gTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGa           <  2:31846/150‑1 (MQ=255)
                                                                                                                                       tGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAg          <  2:858432/150‑1 (MQ=255)
                                                                                                                                        gCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAgg         <  2:317201/150‑1 (MQ=255)
                                                                                                                                           gAGCCCGCTGCGGCGGGCGGCCTGGCTGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTaa      <  1:650499/150‑1 (MQ=255)
                                                                                                                                            cgCCCGCCGCAGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGt                                                                                             <  2:1109759/61‑1 (MQ=255)
                                                                                                                                            cgCCCGCCGCAGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGt                                                                                             >  1:1109759/2‑62 (MQ=255)
                                                                                                                                            cgCCCGCCGCAGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAAc     >  1:1233586/2‑150 (MQ=255)
                                                                                                                                            cgCCCGCCGCAGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAAc     >  2:849165/2‑150 (MQ=255)
                                                                                                                                            aGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAAc     >  2:783584/1‑150 (MQ=255)
                                                                                                                                              cccGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCa                       >  1:81223/1‑130 (MQ=255)
                                                                                                                                              cccGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCa                       <  2:81223/130‑1 (MQ=255)
                                                                                                                                               ccGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCAAAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCa  <  2:288694/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGGTGCTGAGCCCGCTGCGGCGGGCGGCCTGGCCGATGTTCATGGCGCCTCCTGGGTGCTGCTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCA  >  NC_002947/683679‑683971

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: