Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,736,097 T→A 5.1% intergenic (+138/+24) sodB → / ← ydhP superoxide dismutase (Fe)/putative transporter YdhP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,736,0970TA 5.1% 302.3 / 5.1 98intergenic (+138/+24)sodB/ydhPsuperoxide dismutase (Fe)/putative transporter YdhP
Reads supporting (aligned to +/- strand):  ref base T (42/51);  new base A (2/3);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

CTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCGGAAAT  >  NC_000913/1735952‑1736246
                                                                                                                                                 |                                                                                                                                                     
cTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATc                                                                                                                                                   >  2:1343046/1‑150 (MQ=255)
 tGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATcc                                                                                                                                                  <  1:501249/150‑1 (MQ=255)
    ataataACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCTGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCTGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAATGTTATCCTTc                                                                                                                                               <  2:1521946/150‑1 (MQ=255)
       ataaCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCATGCCTAAGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGc                                                                                                                                            >  1:1666212/1‑150 (MQ=255)
        taaCTGATTGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCg                                                                                                                                           <  2:602592/150‑1 (MQ=255)
         aaCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGt                                                                                                                                          <  1:446642/150‑1 (MQ=255)
         aaCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGt                                                                                                                                          >  1:1271567/1‑150 (MQ=255)
          aCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCGtt                                                                                                                                         <  1:236414/150‑1 (MQ=255)
          aCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCGtt                                                                                                                                         <  1:876830/150‑1 (MQ=255)
            tGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGtttt                                                                                                                                       <  1:1263162/150‑1 (MQ=255)
            tGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGtttt                                                                                                                                       <  1:261279/150‑1 (MQ=255)
                 gCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGttt                                                                                                                                  >  1:361908/1‑150 (MQ=255)
                 gCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGttt                                                                                                                                  >  1:125667/1‑150 (MQ=255)
                 gCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGttt                                                                                                                                  >  1:598835/1‑150 (MQ=255)
                 gCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGttt                                                                                                                                  <  2:893144/150‑1 (MQ=255)
                   aaaTGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTaa                                                                                                                                >  1:707650/1‑150 (MQ=255)
                    aaTGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAAt                                                                                                                               <  2:125667/150‑1 (MQ=255)
                      ttcatcaTTGCCTGAAGGGCTACGCTTATCAGGGATACGCGGATCATAGATGTAGGTCGGTTAAGGCACTCGCCGCATTCGTTAAGATTAATGGGACGTTGTCAGCATCTGTAACGCAGAAGGTTATCCTTTTGCGTTTTTGATTATTTa                                                                                                                             <  2:595799/145‑1 (MQ=255)
                       tcagcaTTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAg                                                                                                                            >  1:609725/2‑150 (MQ=255)
                        cagcaTTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGc                                                                                                                           >  2:1433762/1‑150 (MQ=255)
                          gcaTTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTg                                                                                                                         >  1:1009491/1‑150 (MQ=255)
                          gcaTTGCCTGAAGCGATACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTg                                                                                                                         <  2:1587421/150‑1 (MQ=255)
                             ttGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCTGATAAGGCACTCTCCGCATCCGGCAAGATAAATCTCACGTTGTCAGCAACTGTAACTCAGAATGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTa                                                                                                                      <  2:575924/150‑1 (MQ=255)
                             ttGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTa                                                                                                                      <  2:570348/150‑1 (MQ=255)
                                   gAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAg                                                                                                                     >  1:627901/1‑145 (MQ=255)
                                   gAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAg                                                                                                                     <  2:627901/145‑1 (MQ=255)
                                     aGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCa                                                                                                              <  1:1513268/150‑1 (MQ=255)
                                        gcTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAAc                                                                                                           <  2:1344477/150‑1 (MQ=255)
                                         cTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACt                                                                                                          >  1:23096/1‑150 (MQ=255)
                                         cTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACt                                                                                                          >  2:1157583/1‑150 (MQ=255)
                                            cGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGtt                                                                                                       >  2:1363818/1‑150 (MQ=255)
                                                tATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAg                                                                                                   <  1:1490333/150‑1 (MQ=255)
                                                 aTCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAgg                                                                                                  <  2:6263/150‑1 (MQ=255)
                                                  tCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgt                                                                                                 >  2:1161032/1‑150 (MQ=255)
                                                  tCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgt                                                                                                 <  1:444089/150‑1 (MQ=255)
                                                    aGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg                                                                                               <  1:579386/150‑1 (MQ=255)
                                                             cgGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGc                                                                                      <  2:1514483/150‑1 (MQ=255)
                                                              ggATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCt                                                                                     <  1:1089295/150‑1 (MQ=255)
                                                               gATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTg                                                                                    <  2:1271567/150‑1 (MQ=255)
                                                               gATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTg                                                                                    <  1:1271854/150‑1 (MQ=255)
                                                                atcatcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGGTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGa                                                                                   <  1:272455/150‑1 (MQ=255)
                                                                 tcatcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGAc                                                                                  <  1:296944/150‑1 (MQ=255)
                                                                  catcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACa                                                                                 <  2:609725/150‑1 (MQ=255)
                                                                    tcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATa                                                                               <  1:1343046/150‑1 (MQ=255)
                                                                     cGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATaa                                                                              <  2:1262537/150‑1 (MQ=255)
                                                                      gATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATaaa                                                                             >  1:305182/1‑150 (MQ=255)
                                                                      gATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTCGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATaaa                                                                             <  1:1105669/150‑1 (MQ=255)
                                                                        tGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAaca                                                                           <  1:1116145/150‑1 (MQ=255)
                                                                           aGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGc                                                                                      >  2:455528/1‑136 (MQ=255)
                                                                           aGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGc                                                                                      <  1:455528/136‑1 (MQ=255)
                                                                            ggTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAg                                                                       >  1:1638363/1‑150 (MQ=255)
                                                                            ggTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAg                                                                       >  1:1303425/1‑150 (MQ=255)
                                                                            ggTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAg                                                                       <  1:38769/150‑1 (MQ=255)
                                                                                 gATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaata                                                                  >  1:1484295/1‑150 (MQ=255)
                                                                                 gATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaata                                                                  >  1:795980/1‑150 (MQ=255)
                                                                                 gATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaata                                                                  >  1:878401/1‑150 (MQ=255)
                                                                                    aaGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATg                                                               <  2:503109/150‑1 (MQ=255)
                                                                                       gCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCa                                                            >  2:674384/1‑150 (MQ=255)
                                                                                       gCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCa                                                            >  1:1092571/1‑150 (MQ=255)
                                                                                        cACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAg                                                           <  2:23096/150‑1 (MQ=255)
                                                                                        cACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAg                                                           <  1:66451/150‑1 (MQ=255)
                                                                                             gccgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGGGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCcg                                                      >  2:797635/1‑150 (MQ=255)
                                                                                             gccgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCcg                                                      >  2:688150/1‑150 (MQ=255)
                                                                                               cgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCcgcg                                                    <  2:361908/150‑1 (MQ=255)
                                                                                                 cATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGAc                                                  <  1:1433762/150‑1 (MQ=255)
                                                                                                   tCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACaa                                                <  2:1303425/150‑1 (MQ=255)
                                                                                                     cGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATc                                              <  2:1009491/150‑1 (MQ=255)
                                                                                                     cGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATc                                              >  1:620695/1‑150 (MQ=255)
                                                                                                        cAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCt                                           <  1:674384/150‑1 (MQ=255)
                                                                                                         aaGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTc                                          <  2:707650/150‑1 (MQ=255)
                                                                                                          aGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTcc                                         <  1:1143057/150‑1 (MQ=255)
                                                                                                            aTAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCa                                       >  1:3301/1‑150 (MQ=255)
                                                                                                            aTAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCa                                       <  1:1201736/150‑1 (MQ=255)
                                                                                                            aTAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCa                                       <  2:305182/150‑1 (MQ=255)
                                                                                                                aTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCAc                                   >  1:255639/1‑150 (MQ=255)
                                                                                                                aTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCAc                                   >  1:1109303/1‑150 (MQ=255)
                                                                                                                   gCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCgg                                <  1:688150/150‑1 (MQ=255)
                                                                                                                     aCGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGca                              <  2:527684/150‑1 (MQ=255)
                                                                                                                      cGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGcac                             >  1:1509373/1‑150 (MQ=255)
                                                                                                                         tGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATATTGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACaaa                          >  2:450271/1‑150 (MQ=255)
                                                                                                                           tCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGc                        >  2:1531916/1‑150 (MQ=255)
                                                                                                                             aGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTg                      <  2:374043/150‑1 (MQ=255)
                                                                                                                              gCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGt                     >  2:1321276/1‑150 (MQ=255)
                                                                                                                               cAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTa                    <  1:1161032/150‑1 (MQ=255)
                                                                                                                                  cTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcc                 <  2:482495/150‑1 (MQ=255)
                                                                                                                                    gTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATccca               >  1:15186/1‑150 (MQ=255)
                                                                                                                                     tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTg                      <  1:1021929/142‑1 (MQ=255)
                                                                                                                                     tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTg                      >  2:1021929/1‑142 (MQ=255)
                                                                                                                                     tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcccag              <  2:255639/150‑1 (MQ=255)
                                                                                                                                     tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcccag              >  1:123352/1‑150 (MQ=255)
                                                                                                                                       aCGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcccagcc            >  2:645853/1‑150 (MQ=255)
                                                                                                                                       aCGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcccagcc            >  2:817639/1‑150 (MQ=255)
                                                                                                                                       aCGCAGAAGGATAACCTTCTGCGTTTTTGGGTAATTAGCTGTTAGCAACGCCAACTGTTTAAGGTTGTTTTCTGGCTGAAATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCGCAGCGCCGGCGCACAGCTGTATcccggcc            >  2:951954/1‑150 (MQ=255)
                                                                                                                                        cGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcccagccc           >  1:1660015/1‑150 (MQ=255)
                                                                                                                                         gCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcccagccca          >  1:1137277/1‑150 (MQ=255)
                                                                                                                                            gAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGAATCCCAGCCCAGCg       >  1:382800/1‑150 (MQ=255)
                                                                                                                                             aaGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCgg      <  1:975517/150‑1 (MQ=255)
                                                                                                                                               ggTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCGGaa    >  1:975794/1‑150 (MQ=255)
                                                                                                                                                 ttATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCGGAAAt  >  1:1171565/1‑150 (MQ=255)
                                                                                                                                                 ttATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTACCGCGACATTCTCTCCCATCAGCGGCACAAAGCTGTATCCTAGCCCAGCTGatat  <  2:598835/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCGGAAAT  >  NC_000913/1735952‑1736246

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