Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,897,910 G→A 7.7% intergenic (‑40/+54) 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,9100GA 7.7% 313.4 / 10.8 103intergenic (‑40/+54)ygcS/ygcUputative transporter YgcS/putative FAD‑containing dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (55/40);  new base A (4/4);  total (59/44)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 104 total aligned reads displayed.

CGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCGTG  >  NC_000913/2897762‑2898059
                                                                                                                                                    |                                                                                                                                                     
 gCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATGAGAACGGCTTCTCGcc                                                                                                                                                     >  2:1937231/1‑150 (MQ=255)
 gCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGcc                                                                                                                                                     <  2:1540469/150‑1 (MQ=255)
  ccGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCt                                                                                                                                                    >  2:1564450/1‑150 (MQ=255)
   cGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTg                                                                                                                                                   <  1:1754252/150‑1 (MQ=255)
   cGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGATAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTg                                                                                                                                                   <  2:1303918/150‑1 (MQ=255)
            aaCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgt                                                                                                                                          <  2:1228940/150‑2 (MQ=255)
             aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGaa                                                                                                                                         <  2:329740/150‑1 (MQ=255)
               cGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgttct                                                                                                                                       <  1:1154520/150‑5 (MQ=255)
               cGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgttct                                                                                                                                       <  1:335819/150‑5 (MQ=255)
                  gggTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGc                                                                                                                                    <  1:842800/149‑1 (MQ=255)
                   ggTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCa                                                                                                                                   >  2:467231/1‑150 (MQ=255)
                        ggTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGcc                                                                                                                              >  2:1745786/1‑150 (MQ=255)
                          tgtgCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgttcttat                                                                                                                                    >  1:820497/1‑135 (MQ=255)
                          tgtgCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgttcttat                                                                                                                                    <  2:820497/142‑8 (MQ=255)
                           gtgCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCg                                                                                                                           >  1:1919112/1‑150 (MQ=255)
                            tgCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGt                                                                                                                          <  2:1837337/150‑1 (MQ=255)
                            tgCGCCGAAAGTGAGCGAAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGt                                                                                                                          >  1:791465/1‑150 (MQ=255)
                             gcgcCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTc                                                                                                                         >  2:246332/1‑150 (MQ=255)
                              cgcCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCa                                                                                                                        <  1:83250/150‑1 (MQ=255)
                              cgcCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCa                                                                                                                        >  2:1630802/1‑150 (MQ=255)
                                  gAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATcgc                                                                                                                    <  1:440092/150‑1 (MQ=255)
                                       tGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAAc                                                                                                               <  2:127469/150‑1 (MQ=255)
                                         aGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGc                                                                                                             >  2:1243852/1‑150 (MQ=255)
                                          gcgcAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCt                                                                                                            >  2:484157/1‑150 (MQ=255)
                                          gcgcAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCt                                                                                                            >  1:723751/1‑150 (MQ=255)
                                           cgcAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTg                                                                                                           >  1:1845186/1‑150 (MQ=255)
                                             cagcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGcc                                                                                                         <  1:1701475/150‑1 (MQ=255)
                                              agcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCt                                                                                                        <  2:1203392/150‑1 (MQ=255)
                                               gcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTg                                                                                                       >  2:4705/1‑150 (MQ=255)
                                               gcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTg                                                                                                       >  2:565040/1‑150 (MQ=255)
                                                   tGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGacac                                                                                                   <  1:1564450/150‑1 (MQ=255)
                                                     cGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACat                                                                                                 >  2:1122011/1‑150 (MQ=255)
                                                      ggCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACatt                                                                                                >  1:468666/1‑150 (MQ=255)
                                                       gCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACatta                                                                                               >  2:1880334/1‑150 (MQ=255)
                                                       gCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACat                                                                                                 <  2:1789552/148‑1 (MQ=255)
                                                       gCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACat                                                                                                 >  1:1789552/1‑148 (MQ=255)
                                                          gTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTAttt                                                                                            <  2:494013/150‑1 (MQ=255)
                                                           tGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTAtttt                                                                                           <  2:1638163/150‑1 (MQ=255)
                                                           tGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTAtttt                                                                                           <  1:1489264/150‑1 (MQ=255)
                                                            gAAAACGGTTAAGCGGTAAATCCTCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTAttttt                                                                                          >  1:1862691/1‑150 (MQ=255)
                                                            gAAAACCGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCAGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTAttttt                                                                                          >  1:150977/1‑150 (MQ=255)
                                                             aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGCACGCTGCCTGACACATTATTTTTc                                                                                         >  2:1635651/1‑150 (MQ=255)
                                                             aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                         >  1:308764/1‑150 (MQ=255)
                                                             aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                         >  2:126940/1‑150 (MQ=255)
                                                                 cGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATGAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACGCATTATTTTTCAATc                                                                                     <  1:1937231/150‑1 (MQ=255)
                                                                 cGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATc                                                                                     <  1:952458/150‑1 (MQ=255)
                                                                       aGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatag                                                                               >  1:1097518/1‑150 (MQ=255)
                                                                        gCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgataga                                                                              >  1:1851568/1‑150 (MQ=255)
                                                                        gCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgataga                                                                              >  1:1505308/1‑150 (MQ=255)
                                                                        gCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgataga                                                                              <  1:1630802/150‑1 (MQ=255)
                                                                        gCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCAGCCTGACACATTATTTTTCAATCGgataga                                                                              <  2:568765/150‑1 (MQ=255)
                                                                          ggTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatagata                                                                            <  2:1078004/150‑1 (MQ=255)
                                                                             aaaTCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTTACACATTATTTTTCAATCGGATAGATAGTa                                                                         >  2:852632/1‑150 (MQ=255)
                                                                             aaaTCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTa                                                                         <  2:1189275/150‑1 (MQ=255)
                                                                                 catcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCg                                                                     >  1:1445838/1‑150 (MQ=255)
                                                                                   tcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTa                                                                   <  2:308764/150‑1 (MQ=255)
                                                                                       ttCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCa                                                               <  2:636859/150‑1 (MQ=255)
                                                                                               ggTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATgaaa                                                       >  2:1373737/1‑147 (MQ=255)
                                                                                               ggTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCa                                                       <  1:467231/150‑1 (MQ=255)
                                                                                                gTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCAt                                                      <  2:1097518/150‑1 (MQ=255)
                                                                                                   aaGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAg                                                   <  2:1919112/150‑1 (MQ=255)
                                                                                                    aGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAgg                                                  >  2:1319308/1‑150 (MQ=255)
                                                                                                    aGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAgg                                                  >  1:1267139/1‑150 (MQ=255)
                                                                                                      tgtTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGAt                                                >  2:419633/1‑150 (MQ=255)
                                                                                                       gtTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATc                                               >  1:953565/1‑150 (MQ=255)
                                                                                                             cTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATa                                                             <  1:652389/130‑1 (MQ=255)
                                                                                                             cTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATa                                                             >  2:652389/1‑130 (MQ=255)
                                                                                                               ttACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCt                                       <  1:228673/150‑1 (MQ=255)
                                                                                                                  ccccTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCtttt                                    >  1:66909/1‑150 (MQ=255)
                                                                                                                  ccccTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCtttt                                    >  1:1877746/1‑150 (MQ=255)
                                                                                                                  ccccTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCtttt                                    >  1:1791530/1‑150 (MQ=255)
                                                                                                                    ccTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTc                                  <  1:4705/150‑1 (MQ=255)
                                                                                                                      ttCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAg                                <  1:126940/150‑1 (MQ=255)
                                                                                                                        cATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGcc                              <  2:1443100/150‑1 (MQ=255)
                                                                                                                          ttGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCt                            >  1:622228/1‑150 (MQ=255)
                                                                                                                              ccGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCa                        >  2:1110162/1‑150 (MQ=255)
                                                                                                                                ggTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGc                      <  1:1319308/150‑1 (MQ=255)
                                                                                                                                    ataaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACg                  <  1:246332/150‑1 (MQ=255)
                                                                                                                                     taaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGc                 >  2:87467/1‑150 (MQ=255)
                                                                                                                                     ccttttcGGCTTCTCAGGCGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGc                 >  2:445842/7‑150 (MQ=255)
                                                                                                                                     ccttttcGGCTTCTCAGGCGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGc                 >  1:272232/7‑150 (MQ=255)
                                                                                                                                      aaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCt                            <  2:381558/138‑1 (MQ=255)
                                                                                                                                      aaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCt                            >  1:381558/1‑138 (MQ=255)
                                                                                                                                      aaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGcc                >  1:1619388/1‑150 (MQ=255)
                                                                                                                                      aaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGcc                >  2:1447332/1‑150 (MQ=255)
                                                                                                                                      aaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGcc                >  2:129273/1‑150 (MQ=255)
                                                                                                                                       aGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGccc               >  1:838414/1‑150 (MQ=255)
                                                                                                                                          tcGGCTTCTCAGGCGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcg            >  2:1470714/2‑150 (MQ=255)
                                                                                                                                          aCGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcg            >  1:1102125/1‑150 (MQ=255)
                                                                                                                                          aCGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcg            >  1:164122/1‑150 (MQ=255)
                                                                                                                                           cGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGTTAGATAGTACCCGTATTCATAATGCCATTTGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcgc           <  2:150977/150‑1 (MQ=255)
                                                                                                                                            ggCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcgcg          >  2:1676631/1‑150 (MQ=255)
                                                                                                                                            ggCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcgcg          >  1:1631136/1‑150 (MQ=255)
                                                                                                                                            ggCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcgcg          >  1:1747421/1‑150 (MQ=255)
                                                                                                                                              cTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCt        >  1:663156/1‑150 (MQ=255)
                                                                                                                                                tctcGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGc      <  1:565040/150‑1 (MQ=255)
                                                                                                                                                 ctcGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGAGGAACGCTGCCTCACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGcc     >  2:240529/1‑150 (MQ=255)
                                                                                                                                                   cGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCgt   >  1:1700531/1‑150 (MQ=255)
                                                                                                                                                   cGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACACCCACGCGCTGCCgt   >  2:76933/1‑150 (MQ=255)
                                                                                                                                                    gCCTGAGAAGCCGAAAAGGCAACGACCCGTCATCGCGTAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGAGCTGCCgtg  >  2:251856/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCGTG  >  NC_000913/2897762‑2898059

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: ‑