Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 781,757 G→T intergenic (+105/‑328) 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,7570GT97.6% 120.4 / ‑5.5 42intergenic (+105/‑328)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (19/22);  minor base A (0/1);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TGGTAGAGAACGTGAAGGATAACGTTGCGTTAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGT  >  NC_000913/781461‑782051
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                      
tggtgaggaACGGGAAGGATAAGGTTGCGTTAGCAACGGCACGAAGGGCGAGACGAAGTGGAGTCATCCTGTACGACCCACCATCCTGAATGAGTAGGGCAGCATAATCCCGGAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCtcgg                                                                                                                                                                                                                                                                                                     <  2:566132/293‑1 (MQ=18)
         aCGTGAAGGATAACGTTGCGTTAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCg                                                                                                                                                                                                                                                                                           >  1:411165/1‑301 (MQ=255)
                        ttGCGTAAGCAACGGCCCGAAGGGCGAAACGAAGACGAGCCAACCTGCACGAACCACCATCCTGAATGATTAAGGCAGAATAAGCCCGCAAGGGGTAGTTAGCTCCGTTGGGAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAGAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGGTACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGagcgaagc                                                                                                                                                                                                                                                                            <  2:628539/301‑1 (MQ=18)
                                  aaCGGCCCGAAGGGCGAGACGAAGACGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCt                                                                                                                                                                                                                                                                  <  2:395879/301‑1 (MQ=255)
                                   aCGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATc                                                                                                                                                                                                                                                                    >  1:803/1‑298 (MQ=255)
                                                    aCGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGAAAGGTGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATg                                                                                                                                                                                                                                                  <  2:650997/299‑1 (MQ=255)
                                                          tCGAGTCATCCTGCACGACCCACCATACTGATTGAGTAAGGCAGCAGAATTCCGAAAGGGTGCGTTAGCTAAGTTGGTAGTGCAGTTGAATTTTAAGCAATTGGTCGCAGGTTCGAATCCTGAACGACCAACGAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGGTGAGCGCAGCGAAAAAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAg                                                                                                                                                                                                                                          <  2:176764/301‑1 (MQ=14)
                                                                cATCCTGCACGACCCACCATCTGGAATGATTAAGGCAACATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATGGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGAGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCACGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCt                                                                                                                                                                                                                                    <  2:64701/301‑1 (MQ=255)
                                                                   ccTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTTCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                                                                                                 >  1:541955/1‑301 (MQ=255)
                                                                       acaaaccagtcaaccTGAATGATTAAGGCAGCAAAATCCCGCAAGGGGTCGTAGGCACATTTGGAAGAGCAGTGGACTGGTAATCAATTGGTCGCAGGTTCGAAGCCTGCACGACAAACCAAGGTAAAAAAGCGCCCTAAAGGGGCTTTTTGGCGATCGGAGATACTCAAAGATTCGAACCTGAAGCAGGTTTGAGTGGAGCGCAGCGAAAGAACGGAGCCGCTCTCGGGGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                                                                                              <  2:411161/287‑1 (MQ=14)
                                                                             ccaccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAGGCGCttttt                                                                                                                                                                                                                       >  1:594448/1‑301 (MQ=255)
                                                                                        aaTGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGACGGGCGAGCGAAGCGAGTCATACTGCACGACCCACCACTGTAAAACAACGCCCTAAAGGCGCTTTTTTGCTCTTCAgg                                                                                                                                                                                                            >  2:189092/1‑301 (MQ=255)
                                                                                                  gcagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAAt                                                                                                                                                                                                  >  1:688341/1‑301 (MQ=255)
                                                                                                       aTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGACGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCTTCAATTTCa                                                                                                                                                                                              >  1:151426/1‑300 (MQ=255)
                                                                                                       aTAAGCCCGCATGTGGTCGTTAGATCCGATGTTAATGAAGTTGAATTTTAATCACTTTGTGGCAGGTTCGAATCCTGCACGACCCACCACTGTAAAAAAGCGCCTTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCATGTTTGAGTTGAGGGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCAc                                                                                                                                                                                             <  2:488382/301‑1 (MQ=14)
                                                                                                           ccccGCAAGGGGTCGTTAGCGCAGTTGGTAGGGCAGTTGAATTTTAAACAGTTGGTCGCAGGTTCGAATACTGCACGACCAACCAATGTAAGAAAGCGCCCTAAAGTCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTg                                                                                                                                                                                         <  2:306115/300‑1 (MQ=255)
                                                                                                                      gTCGCTAGCTCAGTTGGTAGAGCAGTTGAATTTTAATCACTAGGTCGCCGGGTCGATTCCTGCACGACACACCAATGTAAAAAAGCGCCCTAAAGACGCTTTTTTGCTATCTGCGATGCTCAACGATTCGAACCGGCAGCAGGTTTGAGTTGGGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtgat                                                                                                                                                                               <  2:616136/300‑1 (MQ=14)
                                                                                                                         gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGACGGGCGAGCGACGCGAGTCATCCTGCACGACCCACCAATGTCAGAAAGCGCACTCAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACa                                                                                                                                                                           >  2:652172/1‑301 (MQ=255)
                                                                                                                           tAGCTAAGTTGTTAAAGAAGTTGACATTTGATCAATGGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATc                                                                                                                                                                         <  2:151424/301‑1 (MQ=255)
                                                                                                                              cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTc                                                                                                                                                                      <  1:644046/301‑1 (MQ=255)
                                                                                                                                    tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCagag                                                                                                                                                                  >  1:638942/1‑299 (MQ=255)
                                                                                                                                           gCAGTTGACTTTTAATCAATCGGTCGCAGGTTCGAATCCCGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCCATCCGCGATACTCAAAGACTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGCGCGAAGCGAGTCCTCCTGCCCGACCCATCAATGTAAAACAGCCCTCTAACGGCGCTTTTTTGCTATTCACGCATCCTCAATTTCACTTTGTAAATCTGTTTACATTGTCAGTTCTTACtctt                                                                                                                                                         >  2:466197/1‑298 (MQ=14)
                                                                                                                                            cAGTTGACTTTTAATCAATTGGTCGCAAGTTCGAATCCTGCACGACCCACCACTGTAAAAAAGCGACCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACtgtg                                                                                                                                                         <  2:638935/300‑1 (MQ=255)
                                                                                                                                               ttGACATTAAATGAATTGGTGGAAGGATCGAAGCCTGCAGGACACACAAATGTAAAAAAGCGCCCTAAGGAGGCTTTGTTGCTAAAGGCGCGACTAAAGGATTCAAACCGCCAGAAGGTTAGAGTGGAGGGAAGCGAAAAAACGGAGGCGCTCTGGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaag                                                                                                                                                     <  2:4196/301‑1 (MQ=17)
                                                                                                                                                                                 ctcttccgatctccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATg                                                                                                                    <  1:322046/288‑1 (MQ=255)
                                                                                                                                                                                             ccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGagatcggaagagg                                                                                                       >  2:322043/1‑289 (MQ=255)
                                                                                                                                                                                                   tAAAAAAGCGACTTACAGGCTCTGTTTTGCTATCTGCGATACTCACAGATTCGAACCTTCAGCAGGTTTGATTTAGGCGCCGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGATCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAg                                                                                                  <  2:541950/300‑1 (MQ=21)
                                                                                                                                                                                                              ctgccgtttgctcttccgatctCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAg                                                                                      <  1:32629/280‑1 (MQ=255)
                                                                                                                                                                                                                          tttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGATTGACTTTTccc                                                                           >  2:66071/1‑300 (MQ=255)
                                                                                                                                                                                                                                    cTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCgagagatcggaagagcgtcgtgtag                                                                >  2:32629/1‑282 (MQ=255)
                                                                                                                                                                                                                                     tGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGa                                                                 >  1:264411/1‑299 (MQ=255)
                                                                                                                                                                                                                                              cAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGttt                                                      >  2:72306/1‑301 (MQ=255)
                                                                                                                                                                                                                                              aaaaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGttt                                                      <  1:189094/300‑1 (MQ=255)
                                                                                                                                                                                                                                                               tgcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAAtt                                     <  1:357319/300‑1 (MQ=255)
                                                                                                                                                                                                                                                               agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAAt                                      <  1:72306/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTa                                    >  2:357316/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                    ccccttccgatctgcAGCGAAACAACGGAGCCGCCCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACAcc                                <  2:388480/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tttGAGATGGACGAAGCGAAACAACGGAGCCGGTCTCGGCGACGTCCCGAAGGGCGAGGGAAGCGAGTCATCCTGCAAGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAGACCTGATGACGTCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACAcc                                <  2:594442/300‑1 (MQ=38)
                                                                                                                                                                                                                                                                           ttgagCGCAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAAc                          <  1:66071/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gcAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACAccagatcggaagag                    >  1:388483/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                 gcAGCGAAACAACGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAata                  >  2:164100/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                        aacaacGGAGCCGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGTGCTTACTGTGCAAGCAACTCTATGTCGGTGGAACTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGc            <  1:466202/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                  ccGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGt  >  2:167518/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                      
TGGTAGAGAACGTGAAGGATAACGTTGCGTTAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGT  >  NC_000913/781461‑782051

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