Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,514,936 T→G 17.6% H41P (CAC→CCC)  rpe ← D‑ribulose‑5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,514,9360TG17.6% 147.5 / 2.5 76H41P (CAC→CCC) rpeD‑ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (21/37);  new base G (11/4);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGACCCGCTTTACAGCCATTTTCTTTAATCAGTTGCAGCGTGCGGTCAACATGCTCGGAGGCTTCTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTAT  >  NC_000913/3514712‑3515161
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
ccggACCCGCTTTTCGGCCTTTTTCTTTAAACAGGTGGAGGGGGCGGGCAACATGCTCGGAGGCTTCTGGAGGAAAGGGAATGATGCTGGCCCCAGCGGCAGCGGAATCAGGCAAAATGCTATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGt                                                                                                                                                                                                                                       <  1:144542/225‑1 (MQ=255)
  aGACCCGCTTTACAGCCATTTTCTTTAATCAGTTGCAGCGTGCGGTCAACATGCTCGGAGGCTTCTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGGGTACGTCGATAGGGGCGGTAAAGCCAGAATTACGCAAGGCTTTTACGACCAATTGCCCCATGTTGAGATTTGGGACAAAGAGGTTa                                                                                                                                                                                                                                     >  1:610626/1‑228 (MQ=255)
        gctgtgcAGGCATTTTCGTTAATAAGGTGCAGCGTGCGGTCAACATGCTAGGAGGCTTCGGGAAGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGAGCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATg                                                                                                                                                                                                                               <  2:219102/222‑1 (MQ=255)
         cTTTACAGCCATTTTCTTTAATCAGTTGCAGCGTGCGGTCAACATGCTCGGAGGCTTCTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGGGTACGTCGATAGGGGCGGTAAAGCCATAGTTACGCAAGGGTTTCAGCACCCTTGGCCCAATCGTCAGATTTGGGACAAAGGGGGTATcccaga                                                                                                                                                                                                                              >  2:252793/1‑224 (MQ=255)
                              tCAGTGGCAGCGTGCGGTCAACCTGGTCGGAGGCTTCTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAg                                                                                                                                                                                                         <  2:475451/228‑1 (MQ=255)
                                  ttGCAGCGTGCGGTCAACATGCTCGGAGGCTTCTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCCAAGTTACGCAAGGGTTTCAGACCCATTGGCCCAATCGTCACGTTGGGGACATCAGGGGTATACATGACGTCcaaaaacccaaaaacaagccc                                                                                                                                                                                                     >  1:556135/1‑211 (MQ=255)
                                        cGTGCGGTCAACATGCTCGGAGGCTTCTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCAAAGTTACGCAAGGATTTCAGCACCATTGGGCCAAACGCAGAATTTGGGACAAAGTGGGTTACCCTGACGGCCAAATggcacacatacagccccacggc                                                                                                                                                                                               >  2:575052/1‑207 (MQ=255)
                                          tGCGGTCAACATGCTCGGAGGCTTCTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCCTAGTTACGCAAGGATTTTCGCACCATTGGCCCAATTGGCGAATTGGGAACATAAGGGTTTTCCATCAACTTAAAAAGCACgcaacccagccccactcccc                                                                                                                                                                                             >  1:485391/1‑209 (MQ=255)
                                           gCGGTCAACATGCTCGGAGGCTTCTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGGGTACGTCGATAGGGGCGGAAATGCCATAGTTAAGCAAGGATTTCTGCACCATTGGCCCCATCGCCAGATTTGGAACAAAAGGGGTTATCCTGACGTACAAAAGAACCACATCAACCCCacatccacg                                                                                                                                                                                           >  1:460688/1‑220 (MQ=255)
                                            cGGTCAACATGCTCGGAGGCTTCTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGGGTACGTCGATAGGGGCGGGAATGCCATAGTTACGCAAGGAGTTCAGCACCAATTGCCCCAACGTCCGATTGGGAACATAAGGGTTATCCCTGACACCAAACAGCACGAAATAAGCGCCCGATGACAgg                                                                                                                                                                                          >  1:258437/1‑228 (MQ=255)
                                                    aTGCTCGGAGGCTTCTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGGGTACGTCGATAGGGGCGGGAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAAACGTCAGATTGGGGACCTATAGGTTATCCATGCAGCCAAAATGCACCACATCCACGCCAACTCCCAGGGGTTTTg                                                                                                                                                                                  >  1:111057/1‑227 (MQ=255)
                                                      gCTCGGAGGGTTCTGGCTGAAAGGTAATGGTGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCg                                                                                                                                                                                <  2:610626/227‑1 (MQ=255)
                                                                  cTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCCTAGTTACGCAAGGATTTCAGCACCATTGGCCCAATAGTCAGATTGGGGACATAGTGGGTATCCATGACGGCAAAATGCACGGAATCAGCGCCAACTGCCAGGGCTTTTGCGGGATCTTCcccc                                                                                                                                                                    >  2:34079/1‑227 (MQ=255)
                                                                    ggttGAAGGGTAATGATGCGGGCCCCGGCGGCAACACAATAAGGCACAATGGGAGCGACGGGTTTCACCATTGGGTTTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAgg                                                                                                                                                                 <  1:403616/225‑1 (MQ=255)
                                                                         aaaGGTAATGATGCTGGCCCCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCggg                                                                                                                                                             <  2:448836/227‑1 (MQ=255)
                                                                          aaGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGGCCCATCGTCAGCTTGGGTACACAGGGGTTATCCATGACGGCAAAATGCACGGCATCAAGGCCCGCCGCCCGGGCTTTTGGGGGATATTTAccacagcgggac                                                                                                                                                           >  1:150368/1‑218 (MQ=255)
                                                                              tAACGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAAt                                                                                                                                                       <  2:488903/228‑1 (MQ=255)
                                                                                 tgctggtGGCTCGAGCGGCAGCGAAATCAGGACCAATGCGATCGACGGGTTTCAACAGAAGGTTGACTGCGATTGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCAACATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCa                                                                                                                                                     <  2:692312/221‑1 (MQ=255)
                                                                                     tccggCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGa                                                                                                                                               <  2:431226/225‑1 (MQ=255)
                                                                                     gtgggCCCCAGGGGCAGCGAAATCAGGCACAATGCGATGGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCAACATTGGCCCAATCGGCAGATTGGGAACATAGGGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCg                                                                                                                                                <  1:317016/224‑1 (MQ=255)
                                                                                        ggCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGGTAGGGGCGGTAAAGCCCAAAGTAACCAAGGGTTTTAGGCCCATTGGGCCCAATGTCAGATTGGGGACATAGTGGTTATCCAAGACGTCACAATGCCACACATAAGCGCCAACAGCCGGGGGCTTTGGGGTATTTTCCCCCCGGGAGGGaaaaaaccgcaaaag                                                                                                                                             >  1:356371/1‑217 (MQ=255)
                                                                                        ggCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAACGATTTCAGCACCCTTGGGCCCATCGTCAGATTGGGAACAAATGGGGTTTCCATGACTGCACAATTCACGACATCAGCCCCCCGCTCCCGGGGCTTTAGGGGATCTTTCCCCCGGCGGGCaaaaaccaccccaaa                                                                                                                                             >  1:599762/1‑217 (MQ=255)
                                                                                                  ggagggAAATCATGCAAAGAGCGAACGACGGGTTTGCCCGACGGGTGTACGTCGATAGGGGCGGTAATGCCAATGTTACGCAAGGATTTCAGCGCCATTGGCCCAATCCTCAGATTGGGAACGTAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTGTCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAggg                                                                                                                                   <  2:142385/222‑1 (MQ=255)
                                                                                                   gggCGAAATCGGGCCAAATGCGATCGTCGGGGTTCCCCATCAGGTGTACGGCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAggggg                                                                                                                                 <  1:575052/226‑1 (MQ=255)
                                                                                                   gaGCGAAATCAGGCACAATGCGATCGGCGGGTTTCCCCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAggggg                                                                                                                                 <  2:485391/227‑1 (MQ=255)
                                                                                                   gaGCGAAAGCAGGCACAATGCGATCGACGGGTTTCGCCATCAGGTGTACGTCGAGAGGGGCGGTAATGCCATAGTTACGCAAGGATGTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAggggg                                                                                                                                 <  2:533824/227‑1 (MQ=255)
                                                                                                   cAGGGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAggggg                                                                                                                                 <  2:114684/228‑1 (MQ=255)
                                                                                                   cAGGGAAATAGGGCACACTGGGATCGCCGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGGCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAggggg                                                                                                                                 <  1:281672/228‑1 (MQ=255)
                                                                                                   aaGCGAAATCAGGCACAATGCGATCGGCGGGTTTCACCATCAGGTTTACGTCGATAGGGGCGGTAAGGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAggggg                                                                                                                                 <  2:150368/227‑1 (MQ=255)
                                                                                                    aGCGAAATCAGGCACGAGGCGATCGACGGGTTTCACCAGCAGGTGTACGTCGATAGGGGCGGTAATGGCATAGTTACGCAAGGATTGCAGCACCAGTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACAGCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAggggg                                                                                                                                 <  2:244271/227‑1 (MQ=255)
                                                                                                    aGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAggggg                                                                                                                                 <  2:111057/227‑1 (MQ=255)
                                                                                                    aGCGAAATCAAGCACAATGCGATCGACGGGTTTCAACATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATTCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAggggg                                                                                                                                 <  2:376522/227‑1 (MQ=255)
                                                                                                        aaaTCAGGCAAAATGCAAACGACGGTTTTTACCATTAGGGTTTCGTCGATAGGGGCGGTAATGCCATAGTTACGTAATGATTTCAGCACCAATGACCCAATCGACAGATTGGGAACATAGGGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCATGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATc                                                                                                                           <  1:549766/228‑1 (MQ=255)
                                                                                                        aaaTCAAGCAAAATGCGATCGCAGGGTTTCACCATCAGGTGTCCGGAGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCAAAGCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATc                                                                                                                           <  2:635969/228‑1 (MQ=255)
                                                                                                                   aaTGCGATCGACGGGTTTAACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGtt                                                                                                                <  2:500380/227‑1 (MQ=255)
                                                                                                                     tGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGGTTTCAGCACCATTGGCCCAATCGGCAGATTGGGAACATAGTGGTTAGCCATGCAGTCAAAATGCACGACATCAGAGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCCAAATCAGCCGACAGCATTGAGGGGGCAATCCAATACTGttt                                                                                                               >  2:281388/1‑226 (MQ=255)
                                                                                                                     tGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGAATTCAGCACCATTGGCCCAATTGGCAGATTGGGAAAATAGTGGGTATCCATGACGGCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTTCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACGGCAGTGAGGGGGCGATCCAATACTGttt                                                                                                               >  2:376552/1‑226 (MQ=255)
                                                                                                                       gcggaatgggtataagagacagggTGTAGGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCAt                                                                                                            <  2:548696/205‑1 (MQ=255)
                                                                                                                           cGACGGGGTTCACCAGCAGGTGTACGTCGGTAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGGCAGATGGGGAACAGAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGc                                                                                                        <  1:372627/227‑1 (MQ=255)
                                                                                                                             aCGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTCCGCAAGGGTTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGGGGTTTTCCATTACGGCCAAAATCCCCAAATCAACGCCCACGGCCAGGGCCTTTTCCGCTACTTTCCCCACGCGGGCAAAACCACCCGCACCAACTAGGGGGGCAATCACACACTGtttacaccgctct                                                                                                     >  1:27453/1‑218 (MQ=255)
                                                                                                                                      aCCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGGCAGATTGGGAACATAGGGGTTATCCATGACGTCAAAAAGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGAGAAATCAGCCGGAAGCATTGAGGGGGGAATCAAAAACCTTTTCAGACCCTTCTCCTTggaa                                                                                            >  1:461032/1‑225 (MQ=255)
                                                                                                                                       ccATCAGGTGCACGTCGATAGGGGCGATAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGaa                                                                                           <  2:129379/228‑1 (MQ=255)
                                                                                                                                                 tACGTCGATAGGGGCTGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTc                                                                                 <  2:461032/228‑1 (MQ=255)
                                                                                                                                                 tACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCCTTGGCCCAATCGTCAGATTGGGAAAAAAGTGGTTATCCATGACGTCAAAATGCACGACCTCAGCGGCAGCTGCCAGGGCCTTTGCGGGATATTACCCCCAGCGGGCAAAATCAGCCCCCAGACTTTGGGGGGCAATACCAAAACGGTTCATCCCCCTTCCCTTGGAAATAAttttttt                                                                                 >  1:120877/1‑227 (MQ=255)
                                                                                                                                                  aCGTCGATAGGGGCGGTAATGCCATAGGTACGCAAGGATTTCAGCCCCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTcg                                                                                <  2:120877/228‑1 (MQ=255)
                                                                                                                                                        aTAGGGGCGGTAATGCCATAGTTACGCCAGGATTTCAGCACCATTGGCCCACTCGTCAGATTGGGAACCTAGGGGTTATCCCTGACCTCAAAATGCACGACATAAGCGCCCGCCTCCCGGGCTTTTGCGGACTCTTCACCCAACCCGGCCAAAACAGCCCACCACATTTAGGGGGGCATCACAACACTTTTCCCACCCCTCTCCCTTAGAATTATTTTTTTCCgg‑‑ggg                                                                        >  1:348665/1‑226 (MQ=255)
                                                                                                                                                        aTAGGGGCGGGAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCCAAATGCACGACATCAGCGCCAGCATCCCAGGCTTTTGCGGTAACTTTACCCAAGCGGGCAAAATCAGCCGGCCAGATTTAGGGGGGAATCACATACTGTTTCATCCGCCTTTCCTCGAGCATTAttttttttcgcggg                                                                          >  1:62158/1‑220 (MQ=255)
                                                                                                                                                                ggTAATGCCATAGGTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCCGATTGGGAACATAGGGGTTAACCAAGACGTCAAAATGCCCGAAATTAGCGCCACATGCCAGGGCTTTTGCGGTAACTTACCCCAAGGCGGCCAAATCAGCCCGACGCATTGGGGGGGCAATCAACTTCTGTTTCATCCGCTTCCCCTTTAGAATTATTTTTTTGCGG‑‑GGTAACCGACt                                                                >  1:538558/1‑228 (MQ=255)
                                                                                                                                                                 gTAATGGCTTAGTTAAGAAAGGAGTTCAGCCCTATTGGCACAATCGTCAGATTGGGAACAGAGGGGTGATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATCCTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTc                                                               <  1:541379/228‑1 (MQ=255)
                                                                                                                                                                          atgtttACGCAAGGATTTCAGCCCCATTGGCCCAAGCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGGGTGTGAAACGACTCCTGGt                                                          <  1:34079/223‑1 (MQ=255)
                                                                                                                                                                             gTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCCGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGCCAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCGCCTTTGGAATTATTTTTTCGCGG‑‑GTGAAACAACTCCTGGTTTGTAc                                                    >  2:560655/1‑228 (MQ=255)
                                                                                                                                                                                 cGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGCGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCgg‑‑g                                                                          >  1:117135/1‑202 (MQ=255)
                                                                                                                                                                                 cGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACaaa                                                 >  2:490454/1‑227 (MQ=255)
                                                                                                                                                                                 cGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTACCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGGTTGTAaaaaa                                                >  1:103121/1‑227 (MQ=255)
                                                                                                                                                                                  gCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑Gt                                                                         <  2:117135/202‑1 (MQ=255)
                                                                                                                                                                                   cAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGTACATAGAGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTCAGGGGGCAATTTAATACTGTTTCATCCGCTTCTCCTTTCGAATTATTTTTTCGCAG‑‑CCGAACCGACTCCTCGTATGTACAAAGc                                               >  2:246731/1‑227 (MQ=255)
                                                                                                                                                                                      ggATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGGCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCcagc                                           <  2:599762/228‑1 (MQ=255)
                                                                                                                                                                                        aTTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACCCCTGGTTTGTACAAAGCcaacag                                         >  2:392026/1‑228 (MQ=255)
                                                                                                                                                                                               cACCATTGGCCCAATCGGCAGATTTGGAACTTAGTGGTTAACCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCTGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTcc                                  <  2:667601/228‑1 (MQ=255)
                                                                                                                                                                                                aCCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACACAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTcc                                  <  2:64703/227‑1 (MQ=255)
                                                                                                                                                                                                  cATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCACCAGTTCGTCCAcc                               >  2:275831/1‑228 (MQ=255)
                                                                                                                                                                                                  cATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GGGAAACGACTCCTGGTTTGTACAAACCCAGCAGttattc                                   >  1:64703/1‑220 (MQ=255)
                                                                                                                                                                                                   aTTGGACAAATCGGCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCAc                                <  2:27453/226‑1 (MQ=255)
                                                                                                                                                                                                    ttGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCCATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑CTGAAACGACTCCTGGTTTGTCCAAAGCCATCAGTTCGTCCACCt                              >  2:372331/1‑227 (MQ=255)
                                                                                                                                                                                                             aTCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCTCGG‑‑GTGAAACGACTCCTGGTTTTTTCAAAGCCAACAGTTCGTCCACCTTTTTAcat                      >  1:372201/1‑225 (MQ=255)
                                                                                                                                                                                                             aTCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAACCAGCAGTTCGTCCACCTTTTTCCgt                      >  2:9815/1‑227 (MQ=255)
                                                                                                                                                                                                                gTCAGTTTGGGAACATAGGTGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTgc                  <  1:379344/228‑1 (MQ=255)
                                                                                                                                                                                                                   aGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCCAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTAGCGG‑‑TGGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTTGTCCACCTTTTTCAGTGGgccgc               >  1:94347/1‑228 (MQ=255)
                                                                                                                                                                                                                      aTTGGGAACATAGTGGTTAACCATGACCTCAAAATGCACGACATCATCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGTCGGGCAACATGAGCCGACAGAATTGATGGGGCAATCAAATACTGTTTCCTCCGCTTCTCCTTGAGAATTATTTTTTCCCGG‑‑TTGAAACACCTCCCGGTTTTTACAAACCCAACATTTCGTACCCCTTTTTACgtttcaccccc             >  1:661514/1‑219 (MQ=255)
                                                                                                                                                                                                                        ttggAACATAGTGGTTAACCATGACCTCAAAATGACAGACAGTAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCAAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGtt           <  2:702086/226‑1 (MQ=255)
                                                                                                                                                                                                                        tGGTAACATAGTGGTTATCCATGACGTCAAAATGCACGGAATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGtt           <  2:46750/228‑1 (MQ=255)
                                                                                                                                                                                                                        tGGGAACATGGTGGTTACCCATGACGTCAAAAGGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCTGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGtt           <  2:460688/228‑1 (MQ=255)
                                                                                                                                                                                                                        tGGAAACATAGTGGTTACCCATGACGTCAAAATGCCCGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGtt           <  2:103121/228‑1 (MQ=255)
                                                                                                                                                                                                                         ttggaaactggTGGTTATCCTTGACGTAAAAATGCCCGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCCCCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGtt           <  1:49467/218‑1 (MQ=255)
                                                                                                                                                                                                                                 tAGTGGGTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCACGG‑‑GTGAAACGCCTCCTGGTTTGTAAAACACCAGGAGTTCGTCCACCTTTTAACGTTAGCCGCCGtttcgggtta   >  1:327746/1‑220 (MQ=255)
                                                                                                                                                                                                                                  aGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCCAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTACCGGTTTGTACAAATCCAACAGTTGGTCCACCCTTCTACGTGTGCCGCCGTGCCTGCTtat  >  2:204487/1‑228 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
CCAGACCCGCTTTACAGCCATTTTCTTTAATCAGTTGCAGCGTGCGGTCAACATGCTCGGAGGCTTCTGGATGAAAGGTAATGATGCTGGCACCAGCGGCAGCGAAATCAGGCACAATGCGATCGACGGGTTTCACCATCAGGTGTACGTCGATAGGGGCGGTAATGCCATAGTTACGCAAGGATTTCAGCACCATTGGCCCAATCGTCAGATTGGGAACATAGTGGTTATCCATGACGTCAAAATGCACGACATCAGCGCCAGCTGCCAGGGCTTTTGCGGTATCTTCACCCAGGCGGGCAAAATCAGCCGACAGAATTGAGGGGGCAATCAAATACTGTTTCATCCGCTTCTCCTTGAGAATTATTTTTTCGCGG‑‑GTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTAT  >  NC_000913/3514712‑3515161

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: