Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 424,873 A→C 6.2% S16A (TCG→GCG) ‡ acpH ← acyl carrier protein phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913424,8730AC 6.2% 336.0 / 10.7 109S16A (TCG→GCG) ‡acpHacyl carrier protein phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (50/52);  new base C (4/3);  total (54/55)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 112 total aligned reads displayed.

AAACCACTCCCGTGCTTCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCATGCGCGTTACC  >  NC_000913/424724‑425022
                                                                                                                                                     |                                                                                                                                                     
aaaCCACTCCCGTGCTTCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGa                                                                                                                                                       >  2:388237/1‑150 (MQ=255)
aaaCCACTCCCGTGCTTCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGa                                                                                                                                                       >  1:185357/1‑150 (MQ=255)
 aaCCACTCCCGTGCTTCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAg                                                                                                                                                      >  1:255404/1‑150 (MQ=255)
        cccGTGCTTCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATTAGCCAGTAAATTGCCGGAAAGCGAGCTTTCcg                                                                                                                                               >  1:1134190/1‑150 (MQ=255)
          cGTGCTTCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCTCGCTTTCcggc                                                                                                                                             <  2:386612/150‑3 (MQ=255)
          cGTGCTTCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCcgcg                                                                                                                                             <  1:621653/150‑1 (MQ=255)
              cTTCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCTCGCTTTCcggcaatt                                                                                                                                         <  1:1073833/150‑7 (MQ=255)
                tCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATgg                                                                                                                                       >  1:617568/1‑150 (MQ=255)
                tCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATgg                                                                                                                                       >  1:780005/1‑150 (MQ=255)
                tCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATgg                                                                                                                                       <  1:1109472/150‑1 (MQ=255)
                 cgcgGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATggg                                                                                                                                      <  1:1222750/150‑1 (MQ=255)
                 cgcgGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATggg                                                                                                                                      <  2:1454821/150‑1 (MQ=255)
                   cgGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCt                                                                                                                                    <  2:523292/150‑1 (MQ=255)
                   cgGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCt                                                                                                                                    <  2:788418/150‑1 (MQ=255)
                   cgGACTTCCGGCAAATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCt                                                                                                                                    <  2:1253882/150‑1 (MQ=255)
                        ttCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATg                                                                                                                               >  1:946544/1‑150 (MQ=255)
                         tCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGc                                                                                                                              >  1:1211792/1‑150 (MQ=255)
                              cAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTGAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTg                                                                                                                         <  2:783306/150‑1 (MQ=255)
                              cAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTg                                                                                                                         >  2:334738/1‑150 (MQ=255)
                              cAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTg                                                                                                                         <  1:193521/150‑1 (MQ=255)
                              cAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTg                                                                                                                         >  1:1004645/1‑150 (MQ=255)
                               aGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGa                                                                                                                        >  1:461609/1‑150 (MQ=255)
                               aGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGATTTCCGCGAGATGGGATAAATGCAGGTGa                                                                                                                        >  2:1162010/1‑150 (MQ=255)
                                  ttGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCt                                                                                                                     <  1:538991/150‑1 (MQ=255)
                                       agtcaATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCTGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTaaaaa                                                                                                                >  2:164416/1‑150 (MQ=255)
                                         tcaATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAAt                                                                                                              >  2:1322905/1‑150 (MQ=255)
                                         tcaATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGGACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAAt                                                                                                              <  2:797110/150‑1 (MQ=255)
                                           aaTACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGTAAAGCGAGCTATCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTc                                                                                                            <  2:1365903/150‑1 (MQ=255)
                                           aaTACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCTCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTc                                                                                                            <  1:1167899/150‑1 (MQ=255)
                                           aaTACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAAtta                                                                                                            >  2:1286523/1‑149 (MQ=255)
                                           aaTACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTc                                                                                                            <  2:949218/150‑1 (MQ=255)
                                           aaTACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTc                                                                                                            >  1:980153/1‑150 (MQ=255)
                                           aaTACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTc                                                                                                            >  1:318183/1‑150 (MQ=255)
                                            atacgtcgatacgtcgATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCa                                                                                                           <  1:151473/150‑1 (MQ=255)
                                            atacgtcgatacgtcgATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCa                                                                                                           <  2:287700/150‑1 (MQ=255)
                                             tacgtcgatacgtcgaTGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCAt                                                                                                          <  2:1188517/150‑1 (MQ=255)
                                               cgtcgatacgtcgatGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTg                                                                                                                         >  2:654700/1‑133 (MQ=255)
                                               cgtcgatacgtcgatGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTg                                                                                                                         <  1:654700/133‑1 (MQ=255)
                                               cgtcgatacgtcgatGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCg                                                                                                        <  2:831472/150‑1 (MQ=255)
                                               cgtcgatacgtcgatGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCg                                                                                                        >  1:1200237/1‑150 (MQ=255)
                                                               cATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAGTTCATCGATTATATTCTATCCaa                                                                                        >  1:301206/1‑150 (MQ=255)
                                                               cATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCaa                                                                                        >  2:311099/1‑150 (MQ=255)
                                                                  atGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAg                                                                                     >  2:238848/1‑150 (MQ=255)
                                                                   tGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAgg                                                                                    >  2:446942/1‑150 (MQ=255)
                                                                   tGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAgg                                                                                    >  1:511185/1‑150 (MQ=255)
                                                                   tGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAgg                                                                                    >  2:569105/1‑150 (MQ=255)
                                                                      aTGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAggggg                                                                                 <  1:1064848/150‑1 (MQ=255)
                                                                      aTGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAggggg                                                                                 <  1:1271176/150‑1 (MQ=255)
                                                                       tGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGt                                                                                <  2:780005/150‑1 (MQ=255)
                                                                          cagccaCGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTaaa                                                                             <  2:255404/150‑1 (MQ=255)
                                                                          cagccaCGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTaaa                                                                             >  2:286198/1‑150 (MQ=255)
                                                                            gccaCGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAgg                                                                           <  2:172180/150‑1 (MQ=255)
                                                                               acgacgTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTg                                                                        >  2:394727/1‑150 (MQ=255)
                                                                                cgacgTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGc                                                                       <  1:164416/150‑1 (MQ=255)
                                                                                cgacgTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGc                                                                       <  2:1238332/150‑1 (MQ=255)
                                                                                 gacgTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCa                                                                      <  1:69432/150‑1 (MQ=255)
                                                                                   cgTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAgg                                                                    <  2:560877/150‑1 (MQ=255)
                                                                                    gTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAggg                                                                   <  1:334738/150‑1 (MQ=255)
                                                                                    gTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAggg                                                                   <  2:318183/150‑1 (MQ=255)
                                                                                       gggCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGgaga                                                                <  2:461609/150‑1 (MQ=255)
                                                                                        ggCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGagag                                                               <  1:311099/150‑1 (MQ=255)
                                                                                                aCTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTTAGGTAAATAATTCTTGGATTATTTTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCgg                                                       <  2:653030/150‑1 (MQ=255)
                                                                                                aCTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCgg                                                       <  1:1286523/150‑1 (MQ=255)
                                                                                                 cTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGc                                                      >  1:1385039/1‑150 (MQ=255)
                                                                                                     cttcAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTa                                                  >  2:1171317/1‑150 (MQ=255)
                                                                                                      ttcAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAg                                                 >  1:734912/1‑150 (MQ=255)
                                                                                                      ttcAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAg                                                 >  1:581332/1‑150 (MQ=255)
                                                                                                      ttcAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAg                                                 <  2:170655/150‑1 (MQ=255)
                                                                                                       tcAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGCGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATCTTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGa                                                >  2:747465/1‑150 (MQ=255)
                                                                                                         aGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACt                                              <  2:222503/150‑1 (MQ=255)
                                                                                                            aTTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTAcc                                           >  1:1220892/1‑150 (MQ=255)
                                                                                                                ccGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTAcc                                           >  2:52015/1‑146 (MQ=255)
                                                                                                                ccGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTAcc                                           <  1:52015/146‑1 (MQ=255)
                                                                                                                ccGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGcc                                       <  2:617568/150‑1 (MQ=255)
                                                                                                                ccGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGcc                                       >  2:630852/1‑150 (MQ=255)
                                                                                                                       cGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTAtt                                >  2:1329557/1‑150 (MQ=255)
                                                                                                                       cGAAATCAGACAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTAtt                                >  2:1170981/1‑150 (MQ=255)
                                                                                                                         aaaTCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTAtttt                              <  1:191979/150‑1 (MQ=255)
                                                                                                                         aaaTCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTAtttt                              <  1:394727/150‑1 (MQ=255)
                                                                                                                         aaaTCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTAtttt                              <  1:1460262/150‑1 (MQ=255)
                                                                                                                           aTCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAg                            >  2:1263540/1‑150 (MQ=255)
                                                                                                                           aTCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAg                            >  1:895108/1‑150 (MQ=255)
                                                                                                                            tCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGt                           >  2:385205/1‑150 (MQ=255)
                                                                                                                              aGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCt                         >  1:1021091/1‑150 (MQ=255)
                                                                                                                                ccAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGa                       <  1:1322905/150‑1 (MQ=255)
                                                                                                                                 cAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAg                      <  1:569105/150‑1 (MQ=255)
                                                                                                                                   gTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTc                    <  1:388237/150‑1 (MQ=255)
                                                                                                                                     aaaTTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAg                  <  1:1170981/150‑1 (MQ=255)
                                                                                                                                        ctcgcGGAAAGCTCGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAg                                                                     >  1:507843/5‑96 (MQ=255)
                                                                                                                                          gCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCa             <  2:1036597/150‑1 (MQ=255)
                                                                                                                                          gCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCa             >  2:1028096/1‑150 (MQ=255)
                                                                                                                                          cgcGGAAAGCTCGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCa             >  1:1168672/3‑150 (MQ=255)
                                                                                                                                           gcGGAAAGCTCGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCCGGGAGAGCGCCCCGGCCCTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCAt            >  2:811111/2‑150 (MQ=255)
                                                                                                                                            cGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCATg           >  1:1350325/1‑150 (MQ=255)
                                                                                                                                            cGGAAAGCGAGCTTTCCGCGAGATGGGCAAAATGCAGGTGAGCTAAAAAATTCATCGATTATAATCTATCCAAAAGCGCCTAAACGTTGCAGCGAGAGCCCACCGGCACTAGACTACCCGCCTCATATTTTAGTCTCAGTCACTGTCATg           >  2:1457047/1‑150 (MQ=255)
                                                                                                                                             ggAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCATgc          <  2:1211792/150‑1 (MQ=255)
                                                                                                                                               aaaGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCATgcgc        >  1:720251/1‑150 (MQ=255)
                                                                                                                                                  gCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCATGCGCGtt     >  1:2964/1‑150 (MQ=255)
                                                                                                                                                  gCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCATGCGCGtt     >  1:920600/1‑150 (MQ=255)
                                                                                                                                                     aGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCATGCGCGTTAcc  <  1:1404164/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACCACTCCCGTGCTTCGCGGACTTCCGGCAGATTGTCAGTCAATACGTCGATACGTCGATGCATATGAATGCCAGCCACGACGTCGGGCGGAAAACTTTCTTCAGGATTTCCGCGTACGAAATCAGCCAGTAAATTGCCGGAAAGCGAGCTTTCCGCGAGATGGGCTAAATGCAGGTGAGCTAAAAAATTCATCGATTATATTCTATCCAAAAGGGGGTAAAGGTTGCAGGGAGAGCGCCCCGGCACTAGACTACCCGCCTCTTATTTTAGTCTGAGTCAGTGTCATGCGCGTTACC  >  NC_000913/424724‑425022

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: