Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,897,911 C→G 6.6% intergenic (‑41/+53) ygcS ← / ← ygcU putative transporter YgcS/putative FAD‑containing dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,897,9110CG 6.6% 363.7 / 12.5 106intergenic (‑41/+53)ygcS/ygcUputative transporter YgcS/putative FAD‑containing dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (50/49);  new base G (2/5);  total (52/54)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

Only 100 of 108 total aligned reads displayed.

CGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCGT  >  NC_000913/2897762‑2898058
                                                                                                                                                     |                                                                                                                                                   
cGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGc                                                                                                                                                     <  2:122783/150‑1 (MQ=255)
 gCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGcc                                                                                                                                                    <  1:942872/150‑1 (MQ=255)
   cGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTg                                                                                                                                                  >  1:30312/1‑150 (MQ=255)
   cGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTg                                                                                                                                                  >  1:959895/1‑150 (MQ=255)
    gagaACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTga                                                                                                                                                 <  1:1425430/150‑1 (MQ=255)
    gagaACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTga                                                                                                                                                 <  1:1372333/150‑1 (MQ=255)
    gagaACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTga                                                                                                                                                 <  1:1330378/150‑1 (MQ=255)
       aaCATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGaa                                                                                                                                              <  1:851050/150‑1 (MQ=255)
         cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGc                                                                                                                                            >  2:523710/1‑150 (MQ=255)
         cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGc                                                                                                                                            <  2:207635/150‑1 (MQ=255)
           tAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCg                                                                                                                                          <  2:57346/150‑1 (MQ=255)
           tAACCGTCGGGCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCg                                                                                                                                          <  1:779019/150‑1 (MQ=255)
            aaCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgt                                                                                                                                         <  2:9443/150‑2 (MQ=255)
            aaCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgt                                                                                                                                         <  1:401987/150‑2 (MQ=255)
             aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgtt                                                                                                                                        <  2:650768/150‑3 (MQ=255)
               cGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgttct                                                                                                                                      <  1:904034/150‑5 (MQ=255)
               cGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgttct                                                                                                                                      <  1:911676/150‑5 (MQ=255)
                gTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAg                                                                                                                                     <  1:974175/150‑1 (MQ=255)
                   ggTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCa                                                                                                                                  >  1:510265/1‑150 (MQ=255)
                    gTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCGGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCaa                                                                                                                                 >  2:1042377/1‑150 (MQ=255)
                    gTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCac                                                                                                                                 >  2:619657/1‑149 (MQ=255)
                     tCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAAc                                                                                                                                <  2:599842/150‑1 (MQ=255)
                      cAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACg                                                                                                                               >  1:340415/1‑150 (MQ=255)
                      cAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACg                                                                                                                               <  1:689627/150‑1 (MQ=255)
                           gtgCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCg                                                                                                                          >  2:528847/1‑150 (MQ=255)
                           gtgCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACACCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCg                                                                                                                          <  2:554460/150‑1 (MQ=255)
                             gcgcCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTc                                                                                                                        >  1:468127/1‑150 (MQ=255)
                              cgcCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCa                                                                                                                       >  1:610500/1‑150 (MQ=255)
                                ccGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATc                                                                                                                     >  2:375647/1‑150 (MQ=255)
                                    aaGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCgg                                                                                                                 <  2:1546767/150‑1 (MQ=255)
                                      gTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGaa                                                                                                               >  2:426235/1‑150 (MQ=255)
                                       tGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAAc                                                                                                              <  2:1233845/150‑1 (MQ=255)
                                        gAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACg                                                                                                             >  2:799578/1‑150 (MQ=255)
                                            gcagcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGc                                                                                                         >  2:378260/1‑150 (MQ=255)
                                             cagcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGcc                                                                                                        <  1:1006599/150‑1 (MQ=255)
                                             cagcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGcc                                                                                                        >  1:1315422/1‑150 (MQ=255)
                                               gcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTg                                                                                                      >  2:572585/1‑150 (MQ=255)
                                               gcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTg                                                                                                      >  1:1620128/1‑150 (MQ=255)
                                                caATGCCGCAGTGAAAACGGATAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGGTCATCAAGACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGa                                                                                                     >  2:1123898/1‑150 (MQ=255)
                                                        cAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACattat                                                                                             >  1:1094590/1‑150 (MQ=255)
                                                         aGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACattatt                                                                                            >  2:441669/1‑150 (MQ=255)
                                                            gAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTAttttt                                                                                         >  2:1606692/1‑150 (MQ=255)
                                                             aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                        <  1:528847/150‑1 (MQ=255)
                                                             aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                        >  2:214183/1‑150 (MQ=255)
                                                             aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                        >  2:259044/1‑150 (MQ=255)
                                                             aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                        <  2:497600/150‑1 (MQ=255)
                                                              aaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCa                                                                                       <  2:30312/150‑1 (MQ=255)
                                                               aaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCaa                                                                                      <  2:1521148/150‑1 (MQ=255)
                                                                  ggTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCg                                                                                   <  1:209547/150‑1 (MQ=255)
                                                                     tAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgat                                                                                <  2:468127/150‑1 (MQ=255)
                                                                      aaGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgata                                                                               >  1:122236/1‑150 (MQ=255)
                                                                      aaGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgata                                                                               <  1:1138495/150‑1 (MQ=255)
                                                                       aGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatag                                                                              >  2:830463/1‑150 (MQ=255)
                                                                       aGCGGTAAATCATCCATTCGCACAGGTGAAGTGTTCATCTTTACACCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCAGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatag                                                                              >  2:1206748/1‑150 (MQ=255)
                                                                         cGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatagat                                                                            <  2:525615/150‑1 (MQ=255)
                                                                          ggTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatagata                                                                           >  1:630704/1‑150 (MQ=255)
                                                                            tAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGt                                                                         >  2:124613/1‑150 (MQ=255)
                                                                             aaaTCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTa                                                                        >  1:1209603/1‑150 (MQ=255)
                                                                                    ccATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgat                                                                                <  2:174162/135‑1 (MQ=255)
                                                                                    ccATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgat                                                                                >  1:174162/1‑135 (MQ=255)
                                                                                     cATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTAtt                                                                >  1:1663166/1‑150 (MQ=255)
                                                                                      aTTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTc                                                               >  1:52575/1‑150 (MQ=255)
                                                                                      aTTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTc                                                               <  1:284017/150‑1 (MQ=255)
                                                                                        tCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCAt                                                             <  2:1094590/150‑1 (MQ=255)
                                                                                         cGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATa                                                            >  1:702168/1‑150 (MQ=255)
                                                                                           cACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAAt                                                          <  1:1606692/150‑1 (MQ=255)
                                                                                              cGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGcc                                                       >  1:1311832/1‑150 (MQ=255)
                                                                                              cGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGcc                                                       >  1:297038/1‑150 (MQ=255)
                                                                                                gTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCAt                                                     >  2:69746/1‑150 (MQ=255)
                                                                                                 tGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCAtt                                                    <  1:214183/150‑1 (MQ=255)
                                                                                                   aaGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAg                                                  <  2:110923/150‑1 (MQ=255)
                                                                                                       gtTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATc                                              <  2:340415/150‑1 (MQ=255)
                                                                                                       gtTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATc                                              <  1:572585/150‑1 (MQ=255)
                                                                                                        ttCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCg                                             <  2:941576/150‑1 (MQ=255)
                                                                                                         tCATCTTTACACCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGa                                            >  2:357887/1‑150 (MQ=255)
                                                                                                           aTCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAAc                                          <  1:441669/150‑1 (MQ=255)
                                                                                                                  ccccTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCtttt                                   <  2:773144/150‑1 (MQ=255)
                                                                                                                    ccTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTc                                 <  2:959895/150‑1 (MQ=255)
                                                                                                                     cTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCa                                <  2:1105094/150‑1 (MQ=255)
                                                                                                                      ttCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAg                               <  1:357887/150‑1 (MQ=255)
                                                                                                                             tCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTcc                        >  1:881769/1‑150 (MQ=255)
                                                                                                                               cGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAg                      <  2:630704/150‑1 (MQ=255)
                                                                                                                                 gTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCa                    <  1:1147309/150‑1 (MQ=255)
                                                                                                                                   aataaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAAc                  >  1:1185639/1‑150 (MQ=255)
                                                                                                                                    ataaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACg                 <  2:702168/150‑1 (MQ=255)
                                                                                                                                     taaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGc                >  1:818085/1‑150 (MQ=255)
                                                                                                                                      aaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGcc               >  2:1434230/1‑150 (MQ=255)
                                                                                                                                      aaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGcc               <  1:378260/150‑1 (MQ=255)
                                                                                                                                      aaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGcc               >  2:158736/1‑150 (MQ=255)
                                                                                                                                      aaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGcc               >  1:1486660/1‑150 (MQ=255)
                                                                                                                                       ttttcGGCTTCTCAGGCGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGccc              >  1:218067/5‑150 (MQ=255)
                                                                                                                                          tcGGCTTCTCAGGCGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcg           >  1:1078101/2‑150 (MQ=255)
                                                                                                                                          aCGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcg           <  2:1194066/150‑1 (MQ=255)
                                                                                                                                              cTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCt       >  1:1504752/1‑150 (MQ=255)
                                                                                                                                                tctcGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGc     >  1:243245/1‑150 (MQ=255)
                                                                                                                                                tctcGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGc     >  1:243239/1‑150 (MQ=255)
                                                                                                                                                tctcGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGc     <  1:334968/150‑1 (MQ=255)
                                                                                                                                                  tcGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCg   >  1:565254/1‑150 (MQ=255)
                                                                                                                                                  tcGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCAATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCg   >  1:196233/1‑150 (MQ=255)
                                                                                                                                                   cGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCgt  >  2:1504582/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCGT  >  NC_000913/2897762‑2898058

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: