Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 836,219 G→T 5.1% intergenic (‑138/+39) PP_0718 ← / ← ychF sulfate transporter/ribosome‑associated potassium‑dependent ATP/GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947836,2190GT 5.1% 311.1 / 2.1 99intergenic (‑138/+39)PP_0718/ychFsulfate transporter/ribosome‑associated potassium‑dependent ATP/GTPase
Reads supporting (aligned to +/- strand):  ref base G (49/42);  new base T (2/3);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 103 total aligned reads displayed.

GGGTTGGGGCATGGTGTCGGCTCTGGGCGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCTTCACCCTTGAACTGGAT  >  NC_002947/836070‑836368
                                                                                                                                                       |                                                                                                                                                     
gggTTGGGGCATGGTGTCGGCTCTGGGCGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCTGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAgcg                                                                                                                                                       >  2:194706/1‑150 (MQ=255)
   ttGGGGCATGGTGTCGGCTCTGGGCGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCaa                                                                                                                                                    <  2:2090404/150‑1 (MQ=255)
     ggggCATGGTGTCGGCTCTGGGCGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGc                                                                                                                                                  >  1:546404/1‑150 (MQ=255)
        gCATGGTGTCGGCTCTGGGCGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTc                                                                                                                                               >  2:2129900/1‑150 (MQ=255)
         ttcaagggtgaagccggtgccaaggaagccggtaagtggcgtctggaaggcaaagactacatcgtcaaagacggcgacgtgatgcacttccgctttaacgtctaagcggtaacctgcagcacacacA‑‑AACCCCTTGAAGCGCAAGCCTCa                                                                                                                                              <  1:233613‑M1/24‑1 (MQ=255)
           tGGTGTCGGCTCTGGGCGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAg                                                                                                                                            >  2:700617/1‑150 (MQ=255)
            aagggtgaagccggtgccaaggaagccggtaagtggcgtctggaaggcaaagactacatcgtcaaagacggcgacgtgatgcacttccgctttaacgtctaagcggtaacctgcagcacacacA‑‑AACCCCTTGAAGCGCAAGCCTCAAgg                                                                                                                                           >  2:664374‑M1/124‑150 (MQ=255)
             agggtgaagccggtgccaaggaagccggtaagtggcgtctggaaggcaaagactacatcgtcaaagacggcgacgtgatgcacttccgctttaacgtctaagcggtaacctgcagcacacacA‑‑AACCCCTTGAAGCGCAAGCCTCAAggg                                                                                                                                          >  2:233613‑M1/123‑150 (MQ=255)
               gtCGGCTCTGGGCGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGt                                                                                                                                        <  2:229869/150‑1 (MQ=255)
                 cGGCTCTGGGCGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGttt                                                                                                                                      <  2:2070907/150‑1 (MQ=255)
                  ggCTCTGGGCGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTtg                                                                                                                                     >  2:2058612/1‑150 (MQ=255)
                      ctGGGCGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTtgtgtg                                                                                                                                 <  1:2058612/150‑1 (MQ=255)
                           cGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGtgctg                                                                                                                            <  1:2334900/150‑1 (MQ=255)
                            gagaAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGtgctgc                                                                                                                           <  1:700617/150‑1 (MQ=255)
                              gaAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAg                                                                                                                         >  1:1473603/1‑150 (MQ=255)
                                aGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGt                                                                                                                       <  1:843508/150‑1 (MQ=255)
                                       tGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCt                                                                                                                <  2:1182112/150‑1 (MQ=255)
                                         taagtggcgtctggaaggcaaagactacatcgtcaaagacggcgacgtgatgcacttccgctttaacgtctaagcggtaacctgcagcacacacA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTtgt                                                                                                                                    <  1:277405‑M1/34‑1 (MQ=255)
                                         taagtggcgtctggaaggcaaagactacatcgtcaaagacggcgacgtgatgcacttccgctttaacgtctaagcggtaacctgcagcacacacA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTtgt                                                                                                                                    >  2:277405‑M1/95‑128 (MQ=255)
                                         aaCCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTa                                                                                                              <  2:1210072/150‑1 (MQ=255)
                                           ccTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCa                                                                                                                          >  2:1889550/1‑136 (MQ=255)
                                           ccTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCa                                                                                                                          <  1:1889550/136‑1 (MQ=255)
                                            cTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGAc                                                                                                           <  2:1474861/150‑1 (MQ=255)
                                             tCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCaa                                                                                                                                             >  1:261974/1‑115 (MQ=255)
                                             tCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCaa                                                                                                                                             <  2:972991/115‑1 (MQ=255)
                                             tCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCaa                                                                                                                                             >  1:972991/1‑115 (MQ=255)
                                             tCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCaa                                                                                                                                             <  2:261974/115‑1 (MQ=255)
                                             tCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACg                                                                                                          >  1:1836659/1‑150 (MQ=255)
                                              cAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGt                                                                                                         <  2:1473603/150‑1 (MQ=255)
                                                      aTGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCgg                                                                                                 <  2:1758656/150‑1 (MQ=255)
                                                        gCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTGAAGGCGGaa                                                                                               >  1:1196972/1‑150 (MQ=255)
                                                          tGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAAGGGGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGt                                                                                             >  2:1958627/1‑150 (MQ=255)
                                                            aaGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGc                                                                                           <  1:1251236/150‑1 (MQ=255)
                                                            aaGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGAGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGc                                                                                           >  1:176897/1‑150 (MQ=255)
                                                              ggCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGAGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCAt                                                                                         <  2:176897/150‑1 (MQ=255)
                                                               gCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATc                                                                                        <  1:1642174/150‑1 (MQ=255)
                                                                   gggCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCACAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGt                                                                                    <  1:1262097/150‑1 (MQ=255)
                                                                         cTTTGCGGCCTGGCACGACGCAAGGCCGCACCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGt                                                                              >  1:777893/1‑150 (MQ=255)
                                                                          tttGCGGCCTGGTACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTc                                                                             >  2:36003/1‑150 (MQ=255)
                                                                               ggCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGa                                                                        >  1:172988/1‑150 (MQ=255)
                                                                               ggCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGa                                                                        >  1:161084/1‑150 (MQ=255)
                                                                                gCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGAc                                                                       >  2:2325204/1‑150 (MQ=255)
                                                                                 ccTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACg                                                                      >  1:1431585/1‑150 (MQ=255)
                                                                                  cTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGa                                                                     >  1:1968825/1‑150 (MQ=255)
                                                                                   tGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCACAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGAt                                                                    >  2:1306962/1‑150 (MQ=255)
                                                                                     gCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGt                                                                  >  2:213955/1‑150 (MQ=255)
                                                                                      cACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTa                                                                 <  1:213955/150‑1 (MQ=255)
                                                                                         gacgCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTc                                                              >  1:991603/1‑150 (MQ=255)
                                                                                           cgCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAACACTACCATGGCCCT‑‑ACCTTTATCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCtt                                                            <  2:1228630/150‑1 (MQ=255)
                                                                                            gCAAGGCCGCCCCTGCCAAGGGCGGGTTTGTTGGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCttt                                                           >  1:1228630/1‑150 (MQ=255)
                                                                                              aaGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCg                                                                                                  <  2:1117004/109‑1 (MQ=255)
                                                                                              aaGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCg                                                                                                  >  1:1117004/1‑109 (MQ=255)
                                                                                              aaGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGc                                                         <  2:875421/150‑1 (MQ=255)
                                                                                              aaGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGc                                                         >  2:1416373/1‑150 (MQ=255)
                                                                                                    gccCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCa                                                   >  2:2113163/1‑150 (MQ=255)
                                                                                                     ccccTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAg                                                  >  2:599574/1‑150 (MQ=255)
                                                                                                      cccTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTtgtgtgt                                                                                                                                <  2:424904/71‑1 (MQ=255)
                                                                                                      cccTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTtgtgtgt                                                                                                                                >  1:424904/1‑71 (MQ=255)
                                                                                                       ccTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGAc                                                <  1:2129900/150‑1 (MQ=255)
                                                                                                       ccTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGAc                                                >  1:2187777/1‑150 (MQ=255)
                                                                                                         tGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGc                                              <  2:546404/150‑1 (MQ=255)
                                                                                                              aGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGa                                                 >  2:2119925/1‑142 (MQ=255)
                                                                                                              aGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGa                                                 <  1:2119925/142‑1 (MQ=255)
                                                                                                              aGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACtt                                         >  2:1280518/1‑150 (MQ=255)
                                                                                                                   gcgcATGTCTGGGTAAAACAA‑‑AACCCCTTGAATCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCgg                                    <  2:455548/150‑1 (MQ=255)
                                                                                                                    cgcATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGc                                   <  2:2154820/150‑1 (MQ=255)
                                                                                                                     gcATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGc                                                                                                   <  2:658332/85‑1 (MQ=255)
                                                                                                                     gcATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGc                                                                                                   >  1:658332/1‑85 (MQ=255)
                                                                                                                       aTGTCTGGGCAACACAA‑‑ACCCCCTGGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTc                                >  2:573582/1‑150 (MQ=255)
                                                                                                                       aTGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTc                                <  1:194706/150‑1 (MQ=255)
                                                                                                                        tGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTcc                               <  1:1620964/150‑1 (MQ=255)
                                                                                                                         gTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCt                              >  1:1252178/1‑150 (MQ=255)
                                                                                                                          tCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCtt                             <  2:172988/150‑1 (MQ=255)
                                                                                                                          tCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCtt                             <  1:84106/150‑1 (MQ=255)
                                                                                                                          ctgcagcacacacA‑‑AACCCCTTGAGGCTTGCGCTTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGaa                                                                                               <  2:2087942‑M1/71‑1 (MQ=255)
                                                                                                                          ctgcagcacacacA‑‑AACCCCTTGAGGCTTGCGCTTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGaa                                                                                               >  1:2087942‑M1/14‑84 (MQ=255)
                                                                                                                          ctgcagcacacacA‑‑AACCCCTTGAGGCTTGCGCTTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCtt                             >  2:247359‑M1/14‑150 (MQ=255)
                                                                                                                          ctgcagcacacacA‑‑AACCCCTTGAGGCTTGCGCTTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCtt                             <  1:247359‑M1/137‑1 (MQ=255)
                                                                                                                              ggCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCa                         >  2:704439/1‑150 (MQ=255)
                                                                                                                               gCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCAt                                                                                         >  1:864029/1‑85 (MQ=255)
                                                                                                                               gCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCAt                                                                                         <  2:864029/85‑1 (MQ=255)
                                                                                                                               gCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTc                                                              >  2:1259250/1‑112 (MQ=255)
                                                                                                                               gCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTc                                                              <  1:1259250/112‑1 (MQ=255)
                                                                                                                                caaaacaa‑‑aaCCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCAcc                       >  1:66484/1‑150 (MQ=255)
                                                                                                                                 aaaacaa‑‑aacCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCg                      >  1:1603814/1‑150 (MQ=255)
                                                                                                                                 aaaacaa‑‑aacCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCg                      >  2:990619/1‑150 (MQ=255)
                                                                                                                                   aacaa‑‑aacCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGc                                                         >  1:1727738/1‑113 (MQ=255)
                                                                                                                                   aacaa‑‑aacCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGc                                                         <  2:1727738/113‑1 (MQ=255)
                                                                                                                                    acaa‑‑aacCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCt                   <  2:1603814/150‑1 (MQ=255)
                                                                                                                                     caa‑‑aacCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAgg                                                                                                                        <  2:2154768/48‑1 (MQ=255)
                                                                                                                                     caa‑‑aacCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAgg                                                                                                                        >  1:2154768/1‑48 (MQ=255)
                                                                                                                                          AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAATTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCTTCAc               >  2:777893‑M1/1‑150 (MQ=255)
                                                                                                                                          AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCTTCAc               >  2:1252178‑M1/1‑150 (MQ=255)
                                                                                                                                             CCCTTGAAGCCCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCTTCACCCt            <  1:497351‑M1/150‑1 (MQ=255)
                                                                                                                                              CCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCTTCACCCtt           >  1:990619‑M1/1‑150 (MQ=255)
                                                                                                                                               CTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCTTCACCCTTg          <  1:92022‑M1/150‑1 (MQ=255)
                                                                                                                                                 TGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCTTCACCCTTGaa        >  2:2187777‑M1/1‑150 (MQ=255)
                                                                                                                                                   AAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCTTCACCCTTGAACt      <  2:362432‑M1/150‑1 (MQ=255)
                                                                                                                                                       gcAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCTTCACCCTTGAACTGGAt  <  1:995938/150‑1 (MQ=255)
                                                                                                                                                       gcAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCTTCACCCTTGAACTGGAt  <  1:36003/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GGGTTGGGGCATGGTGTCGGCTCTGGGCGAGAAGGGCGCTGAACCTCAGCTTAGATGCTGAAGGCGGGGGCCGCTTTGCGGCCTGGCACGACGCAAGGCCGCCCCTGCCAAGGGCGCGCATGTCTGGGCAAAACAA‑‑AACCCCTTGAAGCGCAAGCCTCAAGGGGTTTGTGTGTGCTGCAGGTTACCGCTTAGACGTTAAAGCGGAAGTGCATCACGTCGCCGTCTTTGACGATGTAGTCTTTGCCTTCCAGACGCCACTTACCGGCTTCCTTGGCACCGGCTTCACCCTTGAACTGGAT  >  NC_002947/836070‑836368

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: