Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,897,913 T→C 6.6% intergenic (‑43/+51) 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,9130TC 6.6% 321.1 / 9.8 106intergenic (‑43/+51)ygcS/ygcUputative transporter YgcS/putative FAD‑containing dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (50/49);  new base C (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 = 9.53e-01

Only 100 of 106 total aligned reads displayed.

CGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCGT  >  NC_000913/2897765‑2898058
                                                                                                                                                    |                                                                                                                                                 
cGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTg                                                                                                                                                  >  1:30312/1‑150 (MQ=255)
cGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTg                                                                                                                                                  >  1:959895/1‑150 (MQ=255)
 gagaACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTga                                                                                                                                                 <  1:1330378/150‑1 (MQ=255)
 gagaACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTga                                                                                                                                                 <  1:1372333/150‑1 (MQ=255)
 gagaACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTga                                                                                                                                                 <  1:1425430/150‑1 (MQ=255)
    aaCATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGaa                                                                                                                                              <  1:851050/150‑1 (MQ=255)
      cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGc                                                                                                                                            <  2:207635/150‑1 (MQ=255)
      cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGc                                                                                                                                            >  2:523710/1‑150 (MQ=255)
        tAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCg                                                                                                                                          <  2:57346/150‑1 (MQ=255)
        tAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCg                                                                                                                                          <  2:584532/150‑1 (MQ=255)
        tAACCGTCGGGCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCg                                                                                                                                          <  1:779019/150‑1 (MQ=255)
         aaCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgt                                                                                                                                         <  1:401987/150‑2 (MQ=255)
         aaCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgt                                                                                                                                         <  2:9443/150‑2 (MQ=255)
          aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgtt                                                                                                                                        <  2:650768/150‑3 (MQ=255)
            cGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgttct                                                                                                                                      <  1:904034/150‑5 (MQ=255)
             gTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAg                                                                                                                                     <  1:974175/150‑1 (MQ=255)
             gTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAg                                                                                                                                     <  1:561688/150‑1 (MQ=255)
              tCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCACACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAgg                                                                                                                                    <  1:1472377/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:689627/150‑1 (MQ=255)
                   cAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACg                                                                                                                               >  1:340415/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                                                                                                      <  2:342901/150‑1 (MQ=255)
                                             caATGCCGCAGTGAAAACGGATAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGGTCATCAAGACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGa                                                                                                     >  2:1123898/1‑150 (MQ=255)
                                                     cAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACattat                                                                                             >  1:1094590/1‑150 (MQ=255)
                                                      aGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACattatt                                                                                            >  2:441669/1‑150 (MQ=255)
                                                        tGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTAtttt                                                                                          <  2:1315422/150‑1 (MQ=255)
                                                         gAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTAttttt                                                                                         >  2:1606692/1‑150 (MQ=255)
                                                          aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                        >  2:259044/1‑150 (MQ=255)
                                                          aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                        >  2:214183/1‑150 (MQ=255)
                                                          aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                        <  2:497600/150‑1 (MQ=255)
                                                          aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                        <  1:528847/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:1138495/150‑1 (MQ=255)
                                                                   aaGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgata                                                                               >  1:122236/1‑150 (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:259044/150‑1 (MQ=255)
                                                                                  cATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTAtt                                                                >  1:1663166/1‑150 (MQ=255)
                                                                                   aTTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTc                                                               <  1:284017/150‑1 (MQ=255)
                                                                                   aTTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTc                                                               >  1:52575/1‑150 (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:297038/1‑150 (MQ=255)
                                                                                           cGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGcc                                                       >  1:1311832/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)
                                                                                                  gtgtTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGa                                                >  1:941576/1‑150 (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)
                                                                                                               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               >  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:243239/1‑150 (MQ=255)
                                                                                                                                             tctcGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGc     >  1:243245/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)
                                                                                                                                                    |                                                                                                                                                 
CGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCGT  >  NC_000913/2897765‑2898058

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