Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,077,282 G→T 6.2% 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 6.2% 312.8 / 5.7 96intergenic (+29/+79)PP_4470/mgtEArc domain‑containing transcriptional regulator/magnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (45/45);  new base T (3/3);  total (48/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.66e-01

CTGTCCCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTG  >  NC_002947/5077137‑5077427
                                                                                                                                                 |                                                                                                                                                 
cTGTCCCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTacg                                                                                                                                               >  1:1285832/1‑148 (MQ=255)
  gTCCCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCt                                                                                                                                             >  2:1678044/1‑150 (MQ=255)
  gTCCCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCt                                                                                                                                             <  1:1472565/150‑1 (MQ=255)
     ccTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGc                                                                                                                                          >  2:1212114/1‑150 (MQ=255)
      cTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCt                                                                                                                                         >  1:690874/1‑150 (MQ=255)
      cTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCt                                                                                                                                         >  1:773199/1‑150 (MQ=255)
      cTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTACGCTGGCt                                                                                                                                         <  1:63929/150‑1 (MQ=255)
        gCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCTTACGCTGGCttt                                                                                                                                       <  2:1285832/150‑1 (MQ=255)
              gcgcGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTtgtg                                                                                                                                 <  1:1284119/150‑1 (MQ=255)
              gcgcGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTtgtg                                                                                                                                 <  1:1309347/150‑1 (MQ=255)
              gcgcGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTtgtg                                                                                                                                 <  1:1358238/150‑1 (MQ=255)
              gcgcGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTtgtg                                                                                                                                 <  1:1130722/150‑1 (MQ=255)
                gcgAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTggg                                                                                                                               >  1:777645/1‑150 (MQ=255)
                 cgAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGc                                                                                                                              <  2:1112481/150‑1 (MQ=255)
                 cgAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGc                                                                                                                              >  1:1028291/1‑150 (MQ=255)
                  gAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCg                                                                                                                             >  1:103567/1‑150 (MQ=255)
                       gctgcAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGa                                                                                                                        >  1:639208/1‑150 (MQ=255)
                          gcAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTg                                                                                                                     >  1:957438/1‑150 (MQ=255)
                           cAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                    >  1:1891377/1‑150 (MQ=255)
                                tttCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCGCGCATCGCTCATGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGgcagca                                                                                                               >  1:800721/1‑150 (MQ=255)
                                tttCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGgcagca                                                                                                               >  1:1313208/1‑150 (MQ=255)
                                  tCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAg                                                                                                             <  2:103567/150‑1 (MQ=255)
                                  tCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAg                                                                                                             <  1:978714/150‑1 (MQ=255)
                                   ccGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGCACATGTAGGAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAgg                                                                                                            <  2:265717/150‑1 (MQ=255)
                                    cGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGc                                                                                                           <  2:773199/150‑1 (MQ=255)
                                    cGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCATCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGTGTACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTTGGCGTCTGACTGGCAGCAAGGc                                                                                                           <  2:670440/150‑1 (MQ=255)
                                     gCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGcc                                                                                                          >  1:429663/1‑150 (MQ=255)
                                     gCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGcc                                                                                                          <  1:541186/150‑1 (MQ=255)
                                     gCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGcc                                                                                                          >  2:778350/1‑150 (MQ=255)
                                        aaCTGTCCCACCGCCAGCGGAACGGCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCAtt                                                                                                       <  1:6963/150‑1 (MQ=255)
                                         aCTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCAttt                                                                                                      <  1:294758/150‑1 (MQ=255)
                                               ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAGAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGc                                                                                                >  1:976792/1‑150 (MQ=255)
                                               ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGc                                                                                                >  2:619237/1‑150 (MQ=255)
                                               ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGc                                                                                                <  1:1804646/150‑1 (MQ=255)
                                                cACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCg                                                                                               >  1:1048982/1‑150 (MQ=255)
                                                cACCGCCAGCAGAACGCATTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCTTACGCTGGCtttt                                                                                                                                      >  2:1150232/1‑111 (MQ=255)
                                                cACCGCCAGCAGAACGCATTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCTTACGCTGGCtttt                                                                                                                                      <  1:1150232/111‑1 (MQ=255)
                                                cACCGCCAGCAGAACGACTTGGTCGCCCTCATCGCTCACGATGGGGAAAGGGCCGCATACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCg                                                                                               >  2:498618/1‑150 (MQ=255)
                                                 aCCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCgg                                                                                              <  1:1678044/150‑1 (MQ=255)
                                                     ccAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAtt                                                                                          >  1:1004592/1‑150 (MQ=255)
                                                     ccAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAtt                                                                                          <  2:1459493/150‑1 (MQ=255)
                                                       agcagAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGt                                                                                        >  1:241213/1‑150 (MQ=255)
                                                        gcagAACGCCTTGGTCGCCCTCATCGCTCCCCCTGTCCTAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGtt                                                                                       <  2:1095875/150‑1 (MQ=255)
                                                           gAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTccc                                                                                    >  1:789972/1‑150 (MQ=255)
                                                            aaCGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAtt                                                                                          <  2:589900/143‑1 (MQ=255)
                                                            aaCGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAtt                                                                                          >  1:589900/1‑143 (MQ=255)
                                                            aaCGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCca                                                                                   >  2:1661342/1‑150 (MQ=255)
                                                               gCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGGCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAg                                                                                <  1:1153073/150‑1 (MQ=255)
                                                               gCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCAGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAg                                                                                >  1:373270/1‑150 (MQ=255)
                                                                ccTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGgca                                                                                                                  >  2:1608749/1‑115 (MQ=255)
                                                                ccTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGgca                                                                                                                  <  1:1608749/115‑1 (MQ=255)
                                                                ccTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAtt                                                                                          >  2:717368/1‑139 (MQ=255)
                                                                ccTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAtt                                                                                          <  1:717368/139‑1 (MQ=255)
                                                                  ttGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGgga                                                                             <  2:1905654/150‑1 (MQ=255)
                                                                   tGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGggag                                                                            <  2:429663/150‑1 (MQ=255)
                                                                         cccTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTgg                                                                      <  2:227279/150‑1 (MQ=255)
                                                                          ccTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGc                                                                     >  2:1258074/1‑150 (MQ=255)
                                                                            tCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCtt                                                                   >  1:182978/1‑150 (MQ=255)
                                                                             cATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGCTTGTTCCCACAGGGAGGGTTGGCTTa                                                                  >  2:155958/1‑150 (MQ=255)
                                                                                  cTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTacaaca                                                             <  2:789972/150‑1 (MQ=255)
                                                                                    cACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAgg                                                           <  2:690874/150‑1 (MQ=255)
                                                                                         tgtgGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAg                                                      >  2:861833/1‑150 (MQ=255)
                                                                                             gAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTg                                                  >  2:650800/1‑150 (MQ=255)
                                                                                                  ggCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAg                                             >  1:1512789/1‑150 (MQ=255)
                                                                                                   gCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGc                                            <  1:861833/150‑1 (MQ=255)
                                                                                                       cAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAACGCCATTTCCCTGCGGGATTCTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCaa                                        <  2:193915/150‑1 (MQ=255)
                                                                                                         gACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAgg                                      <  2:1512789/150‑1 (MQ=255)
                                                                                                           cGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTcc                                                                                     <  1:191470/101‑1 (MQ=255)
                                                                                                           cGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTcc                                                                                     >  2:191470/1‑101 (MQ=255)
                                                                                                            gCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAg                                   >  1:542860/1‑150 (MQ=255)
                                                                                                              ctcctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGAt                                 >  2:1593576/1‑150 (MQ=255)
                                                                                                              ctcctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGAt                                 >  1:1200748/1‑150 (MQ=255)
                                                                                                                 ctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATgaa                              <  1:1346927/150‑1 (MQ=255)
                                                                                                                 ctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATgaa                              <  1:650800/150‑1 (MQ=255)
                                                                                                                    aGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATgaa                              >  1:601388/1‑147 (MQ=255)
                                                                                                                    aGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATgaa                              <  2:601388/147‑1 (MQ=255)
                                                                                                                          gCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAaccacc                     <  1:1548730/150‑1 (MQ=255)
                                                                                                                             gCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCt                  >  2:183048/1‑150 (MQ=255)
                                                                                                                               gAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGt                <  1:1593576/150‑1 (MQ=255)
                                                                                                                                aaCACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAt                                                                                           >  2:814290/1‑74 (MQ=255)
                                                                                                                                aaCACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAt                                                                                           <  1:814290/74‑1 (MQ=255)
                                                                                                                                aaCACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTc               <  1:619237/150‑1 (MQ=255)
                                                                                                                                 acacGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTNTCg              <  2:1200748/150‑1 (MQ=255)
                                                                                                                                  cacaAAAAAGCCAGCTTACGCTGGCTTTTTCGTGGGTGTATGACTGGGAATGCGCGGAGTTGCCTTCGGGAGTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCgg             >  1:515289/1‑150 (MQ=255)
                                                                                                                                     gAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTg                                                                       >  2:622178/1‑89 (MQ=255)
                                                                                                                                     gAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTg                                                                       <  1:622178/89‑1 (MQ=255)
                                                                                                                                     gAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTca          >  2:571873/1‑150 (MQ=255)
                                                                                                                                     gAAAAAGCCAGCGTAAGCTGGATGTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTca          >  1:824895/1‑150 (MQ=255)
                                                                                                                                       aaaaGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTcaca        >  1:1694453/1‑150 (MQ=255)
                                                                                                                                       aaaaGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTcaca        <  1:832280/150‑1 (MQ=255)
                                                                                                                                       aaaaGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTcaca        <  1:1132741/150‑1 (MQ=255)
                                                                                                                                          aGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCAcagca     <  2:139179/150‑1 (MQ=255)
                                                                                                                                           gCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGCCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCAcagcag    >  2:1083084/1‑150 (MQ=255)
                                                                                                                                             cAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGgga                                                                             <  1:101139/75‑1 (MQ=255)
                                                                                                                                             cAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGgga                                                                             >  2:101139/1‑75 (MQ=255)
                                                                                                                                             cAGCGTAAGCTGGCTTTTTTGTCGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTg  <  2:957438/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTGTCCCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTG  >  NC_002947/5077137‑5077427

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: