Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,418,251 A→C 5.9% *237S (TAA→TCA)  purC → phosphoribosylaminoimidazole‑succinocarboxamide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,418,2510AC 5.9% 249.2 / 5.9 85*237S (TAA→TCA) purCphosphoribosylaminoimidazole‑succinocarboxamide synthase
Reads supporting (aligned to +/- strand):  ref base A (41/39);  new base C (2/3);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

GCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGTTGTACTGGCAA  >  NC_002947/1418105‑1418398
                                                                                                                                                  |                                                                                                                                                      
gCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAg                                                                                                                                                     >  8:326056/1‑149 (MQ=255)
    ggCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGc                                                                                                                                                 <  6:116269/149‑1 (MQ=255)
     ccgatctgTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGc                                                                                                                                                 <  5:70222/141‑1 (MQ=255)
        acgaGTTCAGCCCGGACGGCTGCTGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTTACAGCGGCACgc                                                                                                                                             <  5:363855/149‑2 (MQ=255)
        acgaGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTTACAGCGGCACgc                                                                                                                                             <  1:329035/149‑2 (MQ=255)
         cgaGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGtt                                                                                                                                            <  5:21506/149‑1 (MQ=255)
         cgaGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGtt                                                                                                                                            >  8:269278/1‑149 (MQ=255)
          gaGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTTACAGCGGCACgccc                                                                                                                                           <  7:155372/149‑4 (MQ=255)
          gaGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTc                                                                                                                                           <  5:326167/149‑1 (MQ=255)
             ttCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGCAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTACGCGGc                                                                                                                                                 >  6:70222/1‑140 (MQ=255)
                aGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCaa                                                                                                                                     <  6:2216/149‑1 (MQ=255)
                   ccGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAgcg                                                                                                                                  >  8:183938/1‑149 (MQ=255)
                     ggACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGcc                                                                                                                                <  7:183938/149‑1 (MQ=255)
                         ggCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGAt                                                                                                                            <  7:324447/149‑1 (MQ=255)
                           cTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATAc                                                                                                                          >  4:297733/1‑149 (MQ=255)
                                       ggACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCCCTGTAACCGGCACGGTGACGCAAGCGCCGGAAACCACGCGAttttt                                                                                                             >  1:173599/1‑149 (MQ=255)
                                        gACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGAtttttt                                                                                                            >  6:138473/1‑149 (MQ=255)
                                         aCAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGAttttttt                                                                                                           <  7:336334/149‑1 (MQ=255)
                                           aaGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTg                                                                                                          >  3:271408/1‑148 (MQ=255)
                                             ggAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGaa                                                                                                        >  3:63686/1‑148 (MQ=255)
                                             ggAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGaa                                                                                                        >  7:66299/1‑148 (MQ=255)
                                                     gCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGAtt                                                                                              >  2:86549/1‑149 (MQ=255)
                                                      cAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGAttt                                                                                             <  3:297733/149‑1 (MQ=255)
                                                      cAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGAttt                                                                                             <  8:66299/149‑1 (MQ=255)
                                                       aagaagATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTg                                                                                            <  4:99563/149‑1 (MQ=255)
                                                       aagaagATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTg                                                                                            >  1:102665/1‑149 (MQ=255)
                                                               ggACAAGGCCCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTa                                                                                     >  6:378194/1‑148 (MQ=255)
                                                                 aCAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGAttt                                                                                             <  7:19468/138‑1 (MQ=255)
                                                                 aCAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGAttt                                                                                             >  8:19468/1‑138 (MQ=255)
                                                                     ggACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGa                                                                                >  3:228658/1‑147 (MQ=255)
                                                                     ggACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAAc                                                                              >  5:203091/1‑149 (MQ=255)
                                                                     ggACCGCATCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAATCGGCACCTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGa                                                                                >  5:203097/1‑147 (MQ=255)
                                                                        ccGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCCCGGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAGGATTTGCGTTTTAAAAGAACGCt                                                                           >  7:221786/1‑149 (MQ=255)
                                                                        ccGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCCCAGGAACGCGCCCCGTTACGCCAGGGCCTGGTTAAACGGCAATTTTTTTGAAAAAAAATTTGGGGTTTAAaacaacacg                                                                           >  2:173663/1‑146 (MQ=255)
                                                                         cGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTg                                                                          <  8:88893/149‑1 (MQ=255)
                                                                          gCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTgg                                                                         <  3:106545/149‑1 (MQ=255)
                                                                           cTTCCGCCAGGGCCTGGTCGACGTGATCGATGCCTACGCAGAAGTTGCCATGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGt                                                                        <  6:203097/149‑1 (MQ=255)
                                                                           cTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGt                                                                        <  6:203091/149‑1 (MQ=255)
                                                                            ttCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTa                                                                       >  2:206821/1‑149 (MQ=255)
                                                                              ccGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGAttttttt                                                                                                           <  3:245556/112‑1 (MQ=255)
                                                                              ccGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGAttttttt                                                                                                           >  4:245556/1‑112 (MQ=255)
                                                                                gCCAGGGCGTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCACGCGCCTGATACCACGCGAGTTTTTTGAAAAAAAGATTTGCGTTTTTAGAGTCCGCTGGTatggt                                                                   >  4:315121/1‑147 (MQ=255)
                                                                                gCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTTGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAACGATTTGCGTTTTAAAAGAACGCTGGTatgat                                                                   <  5:28303/149‑1 (MQ=255)
                                                                                   aGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGAttt                                                                                             >  5:293930/1‑120 (MQ=255)
                                                                                   aGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGAttt                                                                                             <  6:293930/120‑1 (MQ=255)
                                                                                   aGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATgcg                                                                >  6:113647/1‑149 (MQ=255)
                                                                                    gggCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATgcgc                                                               >  5:174722/1‑149 (MQ=255)
                                                                                        cTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCa                                                           <  5:138473/149‑1 (MQ=255)
                                                                                           ggCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCAtt                                                         <  4:271408/148‑1 (MQ=255)
                                                                                           ggCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTg                                                        <  6:174722/149‑1 (MQ=255)
                                                                                            gCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTg                                                        >  7:729/1‑148 (MQ=255)
                                                                                                cGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGgagag                                                   <  7:326056/149‑1 (MQ=255)
                                                                                                   gATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCCCGGTAACGCAAGCGCCGGATACCCCGCGATTTTTTTGAAAAAAAGATTTTCGGTTTTAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATg                                                >  6:39664/1‑149 (MQ=255)
                                                                                                   gATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGGTCACGCAAGCGCCTGATACAACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCGGGGATGATGCGCGGGATTGGAGCGATg                                                >  4:333134/1‑149 (MQ=255)
                                                                                                    aTCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAATGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGAt                                                 <  1:86549/148‑1 (MQ=255)
                                                                                                    aTCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGc                                               >  2:274028/1‑149 (MQ=255)
                                                                                                    aTCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGc                                               >  1:61360/1‑149 (MQ=255)
                                                                                                         aGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAg                                          <  5:276218/149‑1 (MQ=255)
                                                                                                            cTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTgg                                       <  5:39664/149‑1 (MQ=255)
                                                                                                             tACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGga                                                      >  4:331858/1‑133 (MQ=255)
                                                                                                             tACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGga                                                      <  3:331858/133‑1 (MQ=255)
                                                                                                                gaagaaGTTGCCAAGCGCCTGGGCGTGCCGCGGTAAGCGGCCCGTTCACGCAAGCGCCTGGTACCACGCGGTTTTTTTGAAAAAAAGATTTGCGGTTTAAAAGGACGCGGGTATGATGCGGGGGATTGGAGAGCTGGCGGAGGGGtg                                     >  8:146641/1‑146 (MQ=255)
                                                                                                                 aagaagTTGCCAAGCGCCTGGGCGTGCCGGGGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTGTAACAGAACGCTGGTAAGATGCGCGGGATTGGAGAGATGCCGGAGTGGTCg                                    >  6:91131/1‑147 (MQ=255)
                                                                                                                   gaagTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACg                                >  7:229371/1‑149 (MQ=255)
                                                                                                                       ttGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGAc                            <  1:173663/149‑1 (MQ=255)
                                                                                                                        tGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAGGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGAc                            <  8:221786/148‑1 (MQ=255)
                                                                                                                          ccAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGa                         >  7:212409/1‑149 (MQ=255)
                                                                                                                           cAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACg                                                                             >  6:258188/1‑96 (MQ=255)
                                                                                                                           cAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACg                                                                             <  5:258188/96‑1 (MQ=255)
                                                                                                                            aaGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTTGAGAGATGCCGGAGTCGTCGACCGGGACGGCtt                       >  6:115396/1‑149 (MQ=255)
                                                                                                                              gcgcCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTTATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCg                     <  4:63686/149‑1 (MQ=255)
                                                                                                                                  gtgaacGTGCCGCTTACAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAAt                 >  3:213317/6‑149 (MQ=255)
                                                                                                                                  cTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAAt                 >  1:331038/1‑149 (MQ=255)
                                                                                                                                     aacGTGCCGCTTACAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGGAGTGGTGAGATGCCGGAGTGGTCGAAGGGGACGGGGTCGGAATCCg              >  8:142889/3‑149 (MQ=255)
                                                                                                                                      gCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGt             <  4:228658/149‑1 (MQ=255)
                                                                                                                                      gCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGt             <  5:113647/149‑1 (MQ=255)
                                                                                                                                       cGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCTCTGCATTGGTGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGtt            <  5:115396/149‑1 (MQ=255)
                                                                                                                                       cGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGtt            <  5:378194/149‑1 (MQ=255)
                                                                                                                                       cGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGtt            <  2:102665/149‑1 (MQ=255)
                                                                                                                                            cGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGTTGTACt       <  1:206821/149‑1 (MQ=255)
                                                                                                                                            cGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGTTGTACt       >  3:45983/1‑149 (MQ=255)
                                                                                                                                             gCTGTAAGCGGCAGGGTCACGCTAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGCTGCCGGAGTGGTCGAACGGGCCGGATTCGAAATCCGTTGTACTg      <  5:91131/149‑1 (MQ=255)
                                                                                                                                                 tAAGCGGCACGTTCACGCAAGCGCCTGATACCCCGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGTTGTACTGGc    <  7:269278/147‑1 (MQ=255)
                                                                                                                                                 tAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGTTGTACTgttgc  >  6:155029/1‑145 (MQ=255)
                                                                                                                                                 tAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGTTGTACTGGCaa  >  2:76161/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGTTGTACTGGCAA  >  NC_002947/1418105‑1418398

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