Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,824 C→T 6.2% intergenic (+172/‑261) lysQ → / → nadA tRNA‑Lys/quinolinate synthase, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913781,8240CT 6.2% 307.8 / 7.9 96intergenic (+172/‑261)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base C (48/42);  new base T (4/2);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

TTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGAC  >  NC_000913/781681‑781971
                                                                                                                                               |                                                                                                                                                   
tttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCTCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                <  1:632549/149‑1 (MQ=17)
tttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                <  3:385012/149‑1 (MQ=32)
tttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                <  3:749806/149‑1 (MQ=32)
                cTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAtt                                                                                                                                <  2:408121/149‑1 (MQ=11)
                 tCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTc                                                                                                                               >  1:755572/1‑149 (MQ=31)
                 tCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTc                                                                                                                               >  1:282720/1‑149 (MQ=31)
                 tCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTc                                                                                                                               >  4:178531/1‑149 (MQ=31)
                 tCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTc                                                                                                                               >  4:39993/1‑149 (MQ=31)
                   aaaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAg                                                                                                                             >  1:177239/1‑149 (MQ=32)
                    aaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAgg                                                                                                                            >  2:371246/1‑149 (MQ=33)
                        ttCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATc                                                                                                                        >  1:337060/1‑149 (MQ=37)
                            aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCa                                                                                                                    >  2:889200/1‑149 (MQ=40)
                            aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCa                                                                                                                    >  3:589180/1‑149 (MQ=40)
                             aCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCaa                                                                                                                   <  1:872634/149‑1 (MQ=40)
                                 gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTc                                                                                                               >  1:680889/1‑149 (MQ=255)
                                  cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCa                                                                                                              <  1:316379/149‑1 (MQ=255)
                                  cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCa                                                                                                              >  4:772995/1‑149 (MQ=255)
                                  cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCa                                                                                                              <  1:2590/149‑1 (MQ=255)
                                   agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCAc                                                                                                             >  3:14521/1‑149 (MQ=255)
                                    gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACt                                                                                                            >  4:504499/1‑149 (MQ=255)
                                    gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACt                                                                                                            <  2:824988/149‑1 (MQ=255)
                                      agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACttt                                                                                                          >  4:72001/1‑149 (MQ=255)
                                      agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACttt                                                                                                          >  1:512688/1‑149 (MQ=255)
                                      agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACttt                                                                                                          <  4:369432/149‑1 (MQ=255)
                                      agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACttt                                                                                                          <  3:178531/149‑1 (MQ=255)
                                         tttGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTa                                                                                                       >  4:688199/1‑149 (MQ=255)
                                            gagttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAc                                                                                                    <  3:435054/149‑1 (MQ=255)
                                              gttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCt                                                                                                  <  2:337060/149‑1 (MQ=255)
                                                    cgcAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAc                                                                                                    <  1:459947/141‑1 (MQ=255)
                                                    cgcAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAc                                                                                                    >  2:459947/1‑141 (MQ=255)
                                                     gcAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACa                                                                                           >  3:373095/1‑149 (MQ=255)
                                                     gcAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACa                                                                                           >  1:87532/1‑149 (MQ=255)
                                                       aGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATc                                                                                         <  2:353199/149‑1 (MQ=255)
                                                        gCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCg                                                                                        >  2:321730/1‑149 (MQ=255)
                                                         cGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGt                                                                                       <  3:163155/149‑1 (MQ=255)
                                                             acaacGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCaga                                                                                   <  2:282720/149‑1 (MQ=255)
                                                              caacGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCagag                                                                                  <  4:157507/149‑1 (MQ=255)
                                                                 cGGAGCCGCTCGCGGCGACGGACCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCtt                                                                               <  2:755572/149‑1 (MQ=255)
                                                                      ccGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACtgt                                                                          <  1:437580/149‑1 (MQ=255)
                                                                              ggcgACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaagcaa                                                                  >  4:844719/1‑149 (MQ=255)
                                                                                 gACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAActc                                                               <  3:775599/149‑1 (MQ=255)
                                                                                 gACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAActc                                                               <  2:177239/149‑1 (MQ=255)
                                                                                 gACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAActc                                                               >  3:268385/1‑149 (MQ=255)
                                                                                    ggCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTAt                                                            <  2:387767/149‑1 (MQ=255)
                                                                                      cccGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGt                                                          >  1:192041/1‑149 (MQ=255)
                                                                                      cccGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGt                                                          <  4:14521/149‑1 (MQ=255)
                                                                                            gggCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTgg                                                    <  4:268385/149‑1 (MQ=255)
                                                                                              gcgagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGaa                                                  <  2:852618/149‑1 (MQ=255)
                                                                                                  gcgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAg                                              <  4:779782/149‑1 (MQ=255)
                                                                                                  gcgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAg                                              >  4:370679/1‑149 (MQ=255)
                                                                                                    gaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGc                                            <  2:30279/149‑1 (MQ=255)
                                                                                                    gaagcgaGGCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGCCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGc                                            >  1:642326/1‑149 (MQ=255)
                                                                                                     aagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtgat                                                                                               >  4:198979/1‑97 (MQ=255)
                                                                                                     aagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtgat                                                                                               <  3:198979/97‑1 (MQ=255)
                                                                                                       gcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTa                                         >  1:322759/1‑149 (MQ=255)
                                                                                                       gcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTa                                         >  1:285644/1‑149 (MQ=255)
                                                                                                         gaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTaaa                                       <  2:447913/149‑1 (MQ=255)
                                                                                                          aGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTaaaa                                      >  4:77846/1‑149 (MQ=255)
                                                                                                          aGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTaaaa                                      <  2:322759/149‑1 (MQ=255)
                                                                                                           gTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAAt                                     >  2:283404/1‑149 (MQ=255)
                                                                                                              aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCg                                                        >  4:481793/1‑127 (MQ=255)
                                                                                                              aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCg                                                        <  3:481793/127‑1 (MQ=255)
                                                                                                              aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGAc                                  <  1:370800/149‑1 (MQ=255)
                                                                                                              aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGAc                                  >  1:637257/1‑149 (MQ=255)
                                                                                                              aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGAc                                  >  2:74721/1‑149 (MQ=255)
                                                                                                              aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGAc                                  <  3:504499/149‑1 (MQ=255)
                                                                                                              aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGAc                                  <  1:223084/149‑1 (MQ=255)
                                                                                                                ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGc                                <  2:680889/149‑1 (MQ=255)
                                                                                                                 cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCa                               >  3:611215/1‑149 (MQ=255)
                                                                                                                  tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCAt                              >  4:31957/1‑149 (MQ=255)
                                                                                                                    cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATcc                            <  2:642326/149‑1 (MQ=255)
                                                                                                                      cGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTg                          <  3:401439/149‑1 (MQ=255)
                                                                                                                       gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCa                                                                                     >  2:564431/1‑89 (MQ=255)
                                                                                                                       gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCa                                                                                     <  1:564431/89‑1 (MQ=255)
                                                                                                                       gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGc                         >  1:633247/1‑149 (MQ=255)
                                                                                                                          ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGcaca                      >  3:566210/1‑149 (MQ=255)
                                                                                                                          ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGcaca                      >  3:215027/1‑149 (MQ=255)
                                                                                                                           caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACAt                     >  3:551200/1‑149 (MQ=255)
                                                                                                                            accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAAt                                                 >  4:590321/1‑120 (MQ=255)
                                                                                                                            accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAAt                                                 <  3:590321/120‑1 (MQ=255)
                                                                                                                            accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACAtt                    >  2:460145/1‑149 (MQ=255)
                                                                                                                            accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACAtt                    >  4:431193/1‑149 (MQ=255)
                                                                                                                               aaTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAgg                 >  2:430097/1‑149 (MQ=255)
                                                                                                                                aTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGc                >  4:73846/1‑149 (MQ=255)
                                                                                                                                  gcaaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGTGCTTTCTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGt              >  3:833770/3‑149 (MQ=255)
                                                                                                                                    aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGc                                >  1:351722/1‑129 (MQ=255)
                                                                                                                                    aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGc                                <  2:351722/129‑1 (MQ=255)
                                                                                                                                     aaaaaGCGCCTTTAGGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAAt           <  4:833770/149‑1 (MQ=255)
                                                                                                                                      aaaaGCGCCTTTAGGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAAtt          >  2:140402/1‑149 (MQ=255)
                                                                                                                                      aaaaGCGCCTTTAGGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAAtt          >  3:528844/1‑149 (MQ=255)
                                                                                                                                      aaaaGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAAtt          >  3:403735/1‑149 (MQ=255)
                                                                                                                                         aGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGa       <  1:283404/149‑1 (MQ=255)
                                                                                                                                         aGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGa       <  3:72001/149‑1 (MQ=255)
                                                                                                                                         aGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGa       <  4:121018/149‑1 (MQ=255)
                                                                                                                                           cgcCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGt     >  4:706601/1‑149 (MQ=255)
                                                                                                                                            gcCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTg    >  2:533143/1‑149 (MQ=255)
                                                                                                                                              ccTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGAc  >  1:429955/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGAC  >  NC_000913/781681‑781971

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