Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,287 G→C P88P (CCC→CCG)  araD ← L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,2870GC97.5% 113.3 / NA 41P88P (CCC→CCG) araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base C (18/21);  minor base A (1/0);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGT  >  NC_000913/66042‑66533
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
gAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATGCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCCCAATGCCGCCAATGGACGGGa                                                                                                                                                                                                                                                     >  1:418595/1‑250 (MQ=255)
    ggACCAGGACGCCGGGCATTTGCGCTTCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTGCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGcc                                                                                                                                                                                                                                                <  1:335307/251‑1 (MQ=255)
               ccGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAAAGCCGCCAATGGACGGGAATGCCTGATAGagc                                                                                                                                                                                                                                       >  1:193923/1‑249 (MQ=255)
                            cTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGGTGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGCCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGCATCGACTGGCCCGCCTGCGCCCAGATGGTGGCGAGGCGGGAGTGTGTGTGCACAATGCCGCCAATGTACGGGCATTCCTGATAGAGCAGTCGGTGAGTTggg                                                                                                                                                                                                                        >  1:259780/1‑250 (MQ=255)
                                aTCGATTCCCTGTTTTTCGAAGGTTTCTACGATGACTTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCgtat                                                                                                                                                                                                                    <  1:128284/251‑3 (MQ=255)
                                aTCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGTTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGGGCACAATGCCCCCAATGGACGGGAATGCCTGATAGAGGCGGCGGTGAGTTggggtat                                                                                                                                                                                                                    >  2:198964/1‑246 (MQ=255)
                                   gATACCCTTTTTGTCGAAGTTTTCTACGATGCCTGTACCGTTTTCACACTCAGCTGCACCGTTGATTGCTGCGTCGGTTATTGTGCGGGTGCAGGGAATGGTGCCGTAGAAAGAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGGGGCGTGGCGCGAGTGTGTGGGCACAATGCCGTCAATGGACTGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCgt                                                                                                                                                                                                                 <  2:36664/251‑2 (MQ=255)
                                            tttttCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGCGGGCt                                                                                                                                                                                                        >  1:429037/1‑251 (MQ=255)
                                                   aaGGTTTCTTCGATGACGTTTCCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGTTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCT‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCgc                                                                                                                                                                                                 <  1:324286/251‑2 (MQ=255)
                                                            aCGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAAc                                                                                                                                                                                        >  1:385792/1‑251 (MQ=255)
                                                                                           caccGTTGATTGCTGCGTCGTTCATTTGGCGGGGGCAGGGAATGGTTCCGTAGAAATAGTCGGCGGTGGG‑GGTGCCGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGCCGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGAc                                                                                                                                                         <  2:170752/248‑1 (MQ=255)
                                                                                            accGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGGGCACAATGCCGCCAATGGACGGGAATGCCTGGTAGAGCAGTCGGTGAGGTGGCGTATCGGGGGAGGGCTTTTTAGCACCTTCAACCCCTT‑CACCGGTTTCGGTGCTAaccaccaccc                                                                                                                                                       >  2:90375/2‑250 (MQ=255)
                                                                                            accGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGGGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATCGAGCAGGCGGTGAGTTGGCGGATCGGGGGGGGGCTTTTTCGCACCTTCAACCACTT‑CCCCGGTTTCGATGCTAACCACGAcc                                                                                                                                                        >  2:420630/2‑250 (MQ=255)
                                                                                             ccGTTGGTTTCTGCGTCAGTCATTTTGCGGGTGCAGGGGATGATGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGTCCGCCTGCGCCCAGATAGAGGCGTGGCGCGCGGGTGTGTGGACCATGCCGCCAATGGACGGGAATGCCTGATAGAGGAGTCGGTGAGGTGGCCTATCGGAGGAGGGGTTTTTCGCACCTTCAACCACTTTCACGGGTTT‑GATGcgaaacccggccc                                                                                                                                                       >  2:80916/1‑238 (MQ=255)
                                                                                             ccGTTGATTTCTGGGTCGGTCATTTTGCGGGGGCAGGGAATGGTGCCTTAGAAAGAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAAGCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCat                                                                                                                                                      <  1:289008/251‑1 (MQ=255)
                                                                                                  gATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAggg                                                                                                                                                                                                          >  1:443060/1‑195 (MQ=255)
                                                                                                                  aTTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGGATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGCTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCt                                                                                                                                 >  1:157203/1‑251 (MQ=255)
                                                                                                                              cAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGcc                                                                                                                     <  1:124492/251‑1 (MQ=255)
                                                                                                                                         gCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGa                                                                                                          <  1:413474/251‑1 (MQ=255)
                                                                                                                                         gCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGa                                                                                                          <  1:198964/251‑1 (MQ=255)
                                                                                                                                               gAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATTGACGCCGGAAGGTTTGATcaca                                                                                                     >  1:133066/1‑250 (MQ=255)
                                                                                                                                                aaaTAGTCGGCT‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACaaa                                                                                                   <  1:379824/251‑1 (MQ=255)
                                                                                                                                                           g‑tgggtggTGCCGGTTGTTGGATGCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGTCGCGAGGGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCCCCTTCAACCACTT‑CACCGGTTTCTATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCt                                                                                        <  1:423751/251‑1 (MQ=255)
                                                                                                                                                                    tgCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGAAGCCGCGCTCGCGATc                                                                                 >  1:293713/1‑250 (MQ=255)
                                                                                                                                                                               gaaaccgcttggcccGCTTGCGCCCGGATAGTGGCGTGGCGTGAGTGTGTGTGCACAATGCCGCCAATGGACGGTAATGCCTGATAGAGCAGTCGGTGAGGTGGCGTATCGGAGGAGGGCTTTTTCGCAGTTTCAACCACTT‑CACCGTGTTCGATGCTAACCACGTCCATAGCGTCATCGCTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGAc                                                                     <  1:427298/239‑1 (MQ=255)
                                                                                                                                                                                       ctggcccGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGcccc                                                             >  2:160563/5‑251 (MQ=255)
                                                                                                                                                                                              gCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAg                                                      <  1:41979/251‑1 (MQ=255)
                                                                                                                                                                                                   cgcCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGa                                                 <  1:160563/251‑1 (MQ=255)
                                                                                                                                                                                                       cAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATAGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGCTCAACGGCGCTGACGTTTCCCCATGTGAGCGTGAcc                                               >  1:88820/1‑249 (MQ=255)
                                                                                                                                                                                                          atagtggcgtggcGCGAGTGTTGGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGtt                                          <  1:438626/248‑1 (MQ=255)
                                                                                                                                                                                                                  gtggcgCGAGGGGGTGTGCACAAGGCCGACAATGGACGGGAATGACTGATAGAGCAGTCGGTGAGTTGGCGTATAGGAGGAGGGCTTTTGCGCACCTTCAACCACTT‑CACCGGTTGCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGTCCCATGTGAGCGTGACCAGGTTATGTTTTg                                  <  1:185370/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggcgCGATTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTgg                                 <  1:114107/250‑1 (MQ=255)
                                                                                                                                                                                                                      cgcgAGTGTGTGTGCCCAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTTAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAg                              <  1:2956/251‑1 (MQ=255)
                                                                                                                                                                                                                      cgagAGTGTGTGTGCACAATGCCGCCAATGGAAGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAGAGACGCAGCGCTCGCGTTCAACGGCGCTGACGTTGCCCCATGGGAGCGTGACAAGGTTATGTTTTGGc                                >  2:180700/4‑249 (MQ=255)
                                                                                                                                                                                                                             ttttgtgctcactGCCGCCAATGGTCGGGAACGCCTGATAATGCAGGCGGTTATTTGGCGTATCGGGGGAGGGCTTTGTCGCACCTTCAACCTCTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATAGACGCCGGAAGGGTTGATCACAAAGACGCCGCGCTCGCGATCAACTGCGCTGCCGTTGCCCCATGTGAGCGTGACCAGGTTATTTTTCGGCAGCGCCAg                        <  2:251870/238‑1 (MQ=255)
                                                                                                                                                                                                                                    cacaATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGcct                >  1:255199/1‑249 (MQ=255)
                                                                                                                                                                                                                                     acaATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGcctc               >  1:110375/1‑248 (MQ=255)
                                                                                                                                                                                                                                          gccgccAATGGACGGGAATGCCTGATAGAGCAGTAGGTGAGTTTGCGGATCGGAGGAGGGCGTTTTCGCACCTTCAACCACTT‑CACCGGGTGCGAGGCTAACCGCGACCATATCGTCATCGGTCATGATGCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACTAGGTTATGTTTTGGCAGCGCCATGTGGGCTTCTAAtt          <  2:426025/251‑2 (MQ=255)
                                                                                                                                                                                                                                              ccAAAGGACGGGAATGCATGATTGAGCAGTCTGTGAGGTGGCGTATCGGGTGTGGGCGTTTTCGCACCTTCAACCACGT‑CACCGGGTTCGATGCTAACCACGACCATCGCGTCAGCGGTCATGATTCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGTGCCGACGTTGCCCCATGTTAGCGTGACCAGGTTATGTTTTGGCAGCTCCAGGTTGGCCTCTAATatct       <  2:235404/250‑4 (MQ=255)
                                                                                                                                                                                                                                                aaTGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCGCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGc    >  1:109344/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGt  >  1:463043/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
GAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGT  >  NC_000913/66042‑66533

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