Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,077,282 G→T 7.8% intergenic (+29/+79) PP_4470 → / ← mgtE Arc domain‑containing transcriptional regulator/magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,077,2820GT 7.8% 290.2 / 8.7 90intergenic (+29/+79)PP_4470/mgtEArc domain‑containing transcriptional regulator/magnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (39/44);  new base T (3/4);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TGTCCCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTG  >  NC_002947/5077138‑5077427
                                                                                                                                                |                                                                                                                                                 
tGTCCCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGc                                                                                                                                              >  1:448062/1‑150 (MQ=255)
    ccTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGc                                                                                                                                          >  2:616169/1‑150 (MQ=255)
     cTGCATGAGCGCGAGTTGCGGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCt                                                                                                                                         <  1:745007/150‑1 (MQ=255)
      tGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCTTACGCTGGCtt                                                                                                                                        <  2:285342/150‑1 (MQ=255)
         aTGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCttttt                                                                                                                                     <  1:939144/150‑1 (MQ=255)
          tGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTACGCTGGCtttttc                                                                                                                                    <  2:1891947/150‑2 (MQ=255)
            aGCGCGGGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTtgt                                                                                                                                  >  2:1532216/1‑150 (MQ=255)
             gcgcGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGCATAAGCCAGCGTAAGCTGGCTTTTTtgtg                                                                                                                                 <  1:52697/150‑1 (MQ=255)
             gcgcGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTtgtg                                                                                                                                 >  2:482263/1‑150 (MQ=255)
              cgcgAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTgg                                                                                                                                <  1:1775128/150‑1 (MQ=255)
                   gTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTc                                                                                                                           >  2:657540/1‑150 (MQ=255)
                    ttGCTGCAGCGTTTCCGCCAACTGCCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCt                                                                                                                          <  1:605991/150‑1 (MQ=255)
                      gctgcAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGa                                                                                                                        >  1:1064574/1‑150 (MQ=255)
                         gcAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTg                                                                                                                     >  2:1278827/1‑150 (MQ=255)
                          cAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGCCTGCTGGGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAAACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                    <  1:100317/150‑1 (MQ=255)
                          cAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                    >  2:654123/1‑150 (MQ=255)
                          cAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                    >  2:288862/1‑150 (MQ=255)
                           aGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGgc                                                                                                                   <  1:288862/150‑1 (MQ=255)
                                 tCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAg                                                                                                             <  1:1532216/150‑1 (MQ=255)
                                   cGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGc                                                                                                           >  1:1480223/1‑150 (MQ=255)
                                    gCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGcc                                                                                                          <  2:1480223/150‑1 (MQ=255)
                                    gCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGcc                                                                                                          >  1:1133653/1‑150 (MQ=255)
                                     ccAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGCGCTAGTCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTCTGGGCGTATGACTGGCAGCAAGGCCa                                                                                                         >  1:1089275/1‑150 (MQ=255)
                                     ccAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGTGTCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGACGCCGTCTAACTGTCAGCAAGGGCa                                                                                                         >  1:120931/1‑150 (MQ=255)
                                     ccAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCa                                                                                                         >  2:1911784/1‑150 (MQ=255)
                                           gTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTccc                                                                                                   >  1:920550/1‑150 (MQ=255)
                                              ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGc                                                                                                <  1:654123/150‑1 (MQ=255)
                                               cACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCg                                                                                               <  1:692722/150‑1 (MQ=255)
                                               cACCGACAGCAGAACGCTTTGGTCGCACTCATCGCTCATGATGTGGAAATGGCCGCAGACGCCTACTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGAGTGGCAGCAAGGCCATTTCCCTGCg                                                                                               <  2:394201/150‑1 (MQ=255)
                                                      agcagAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGt                                                                                        >  2:1909879/1‑150 (MQ=255)
                                                      agcagAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTACGCTGGCtttttc                                                                                                                                    <  1:1232847/106‑2 (MQ=38)
                                                      agcagAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTACGCTGGCtttttc                                                                                                                                    >  2:1232847/1‑105 (MQ=38)
                                                         agAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTcc                                                                                     <  1:482263/150‑1 (MQ=255)
                                                          gAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGa                                                                                                                        <  1:1863687/114‑1 (MQ=255)
                                                          gAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGa                                                                                                                        >  2:1863687/1‑114 (MQ=255)
                                                           aaCGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCca                                                                                   >  2:1064325/1‑150 (MQ=255)
                                                           aaCGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCca                                                                                   >  2:1541967/1‑150 (MQ=255)
                                                            aCGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCcac                                                                                  <  2:1133653/150‑1 (MQ=255)
                                                             cGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCcaca                                                                                 <  1:657540/150‑1 (MQ=255)
                                                              gCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTgg                                                                                                                                >  2:790115/1‑102 (MQ=255)
                                                              gCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTgg                                                                                                                                <  1:790115/102‑1 (MQ=255)
                                                                 ttGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGgga                                                                             <  1:1541967/150‑1 (MQ=255)
                                                                       gCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACt                                                                                                                      <  2:1308484/103‑1 (MQ=255)
                                                                       gCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACt                                                                                                                      >  1:1308484/1‑103 (MQ=255)
                                                                         ccTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGAtt                                                                        >  1:1430575/1‑147 (MQ=255)
                                                                         ccTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGAtt                                                                        <  2:1430575/147‑1 (MQ=255)
                                                                         ccTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGc                                                                     >  1:459229/1‑150 (MQ=255)
                                                                          cTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCCGCGGGATTGTTCCCACAGGGAGGATTGGCt                                                                    >  1:793509/1‑150 (MQ=255)
                                                                            cATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCATCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTa                                                                  <  1:1005183/150‑1 (MQ=255)
                                                                            cATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTa                                                                  <  1:732984/150‑1 (MQ=255)
                                                                             aTCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTac                                                                 <  1:616169/150‑1 (MQ=255)
                                                                                 cTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCCGCGGGATTGTTCCCACAGGGAGGATTGGCTTacaaca                                                             <  2:793509/150‑1 (MQ=255)
                                                                                  tCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGTATTGTTCCCACAGGGAGGATTGGCTTACAACAg                                                            >  1:1583179/1‑150 (MQ=255)
                                                                                  tCACGATGTGGAAATGGCCGCAGACGCCTACTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAg                                                            >  1:586806/1‑150 (MQ=255)
                                                                                        tgtgGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAg                                                      >  1:1663116/1‑150 (MQ=255)
                                                                                          tgGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                    <  2:1312473/86‑1 (MQ=255)
                                                                                          tgGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                    >  1:1312473/1‑86 (MQ=255)
                                                                                            gAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTg                                                  >  1:1464160/1‑150 (MQ=255)
                                                                                            gAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTg                                                  <  2:920550/150‑1 (MQ=255)
                                                                                              aaTGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGc                                                <  2:1064574/150‑1 (MQ=255)
                                                                                                tGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCa                                              <  2:1464160/150‑1 (MQ=255)
                                                                                                   ccGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTACGCTgg                                                                                                                                           <  2:1458174/54‑1 (MQ=21)
                                                                                                   ccGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTACGCTgg                                                                                                                                           >  1:1458174/1‑54 (MQ=21)
                                                                                                      cAGCCGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGGGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGCAAGAACGTGGCCAGCCTaa                                        <  2:120931/150‑1 (MQ=255)
                                                                                                      cAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCTTGTCCAGCCCaa                                        <  2:1089275/150‑1 (MQ=255)
                                                                                                       aGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGt                                                                                        >  2:552772/1‑101 (MQ=255)
                                                                                                       aGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGt                                                                                        <  1:552772/101‑1 (MQ=255)
                                                                                                        gACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAgg                                      <  2:459229/150‑1 (MQ=255)
                                                                                                         aCGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCGATGTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGa                                     >  1:1393743/1‑150 (MQ=255)
                                                                                                          cGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGaa                                    <  1:1064325/150‑1 (MQ=255)
                                                                                                             ctcctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGAt                                 >  1:152311/1‑150 (MQ=255)
                                                                                                             ctcctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGAt                                 >  1:1371499/1‑150 (MQ=255)
                                                                                                             ctcctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGAt                                 <  2:1583179/150‑1 (MQ=255)
                                                                                                              tcctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATg                                >  1:268521/1‑150 (MQ=255)
                                                                                                                ctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATgaa                              <  2:268521/150‑1 (MQ=255)
                                                                                                                ctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATgaa                              <  2:1371499/150‑1 (MQ=255)
                                                                                                                  gAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATgaaga                            <  2:1663116/150‑1 (MQ=255)
                                                                                                                  gAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATgaaga                            >  1:343811/1‑150 (MQ=255)
                                                                                                                   aGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGc                                                     >  1:1423552/1‑124 (MQ=255)
                                                                                                                   aGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGc                                                     <  2:1423552/124‑1 (MQ=255)
                                                                                                                             cGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGc                                                                     <  2:965203/98‑1 (MQ=255)
                                                                                                                             cGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGc                                                                     >  1:965203/1‑98 (MQ=255)
                                                                                                                               aaCACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAg                                       <  2:233819/126‑1 (MQ=255)
                                                                                                                               aaCACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAg                                       >  1:233819/1‑126 (MQ=255)
                                                                                                                               aaCACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTc               <  1:1911784/150‑1 (MQ=255)
                                                                                                                                ccacaaaaaaGCCAGCTTACGCTGGCTTTTTCGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCg              >  2:1242793/6‑150 (MQ=255)
                                                                                                                                  acGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTg                                                  <  1:8948/112‑1 (MQ=255)
                                                                                                                                  acGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTg                                                  >  2:8948/1‑112 (MQ=255)
                                                                                                                                    gAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTca          <  1:1909879/150‑1 (MQ=255)
                                                                                                                                            cAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTg  <  2:586806/150‑1 (MQ=255)
                                                                                                                                                gTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTacaaca                                                             <  2:1535228/87‑1 (MQ=255)
                                                                                                                                                gTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTacaaca                                                             >  1:1535228/1‑87 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TGTCCCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTG  >  NC_002947/5077138‑5077427

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