Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,295,699 T→G 22.3% V407G (GTG→GGG)  gltP → glutamate/aspartate:proton symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,295,6990TG22.3% 164.1 / 5.6 77V407G (GTG→GGG) gltPglutamate/aspartate:proton symporter
Reads supporting (aligned to +/- strand):  ref base T (11/45);  new base G (18/3);  total (29/48)
Fisher's exact test for biased strand distribution p-value = 1.55e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG  >  NC_000913/4295474‑4295920
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                
gCATTGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTg                                                                                                                                                                                                                                 <  2:475059/227‑1 (MQ=255)
 gatttgccTTTCCATTTGGCAGGAAATCATTGTGGTTCTGACGCTGATGGTGTACTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAt                                                                                                                                                                                                                               <  1:577937/222‑1 (MQ=255)
 cATTGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGGGAt                                                                                                                                                                                                                               >  2:108057/1‑228 (MQ=255)
   ttGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAACGTGTTTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGCCATGGCGCGTACTGCGCTGAACGTGGTGGGTATTGCGGTGGCGGTGCTGGTGGT‑Tg                                                                                                                                                                                                                            >  2:585305/1‑228 (MQ=255)
   ttGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGGGCGGGGGAT‑Tg                                                                                                                                                                                                                            >  2:194733/1‑228 (MQ=255)
           tCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCCGGGGAT‑TGCCaaagtg                                                                                                                                                                                                                    >  2:67554/1‑224 (MQ=255)
           tCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGGATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGGCCTGGGGAT‑TGCCAAggtg                                                                                                                                                                                                                    >  2:577937/1‑225 (MQ=255)
                cTGGCAGGAAATCATTCTGGTTCTGACGCTGGTGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAac                                                                                                                                                                                                                <  2:186552/227‑1 (MQ=255)
                     aGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGCATCCTCCACCTGGCGCGACCTCGGCTGCACGGGGTGGGTAGTCGGCTGCCGGGGAGGGGGATACGACAAGT‑GGAACACa‑aaa                                                                                                                                                                                                         >  1:85010/1‑227 (MQ=255)
                       gAAATCTTTCTGGGTCTGACGGTTTGTGTGACCTCGAAAGGGATTGGTGGCGTGCCGGGGGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AAttt                                                                                                                                                                                                       <  1:67554/228‑1 (MQ=255)
                         aaTCTTTCTGGGGTGGACGTGGATGGTGACTGCGAAGGGGAGTGCGTGCGTGGCTGGCGTGTCGTTTGTGGTGTTGCTGGTAACGTTGGGTAGCGTAGGTATCCCGCTGGAAGGTGTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCGGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGa                                                                                                                                                                                                     <  1:394700/228‑1 (MQ=255)
                         aaTCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGGGAG‑TGGCAATGGGGGACACA‑AAtttta                                                                                                                                                                                                     >  2:548615/1‑226 (MQ=255)
                                   gTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGGGCTGTGGAT‑TGCCAAGTGGGAA‑ACACAATTTGACCGTAAGaaa                                                                                                                                                                                           >  2:155605/1‑228 (MQ=255)
                                    tGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACCGCGCTGAACGTGGTGGGAATTGCGCTGGCGGGGCGGGGGGT‑TTCCAAGGGGGACACCA‑AAATTTACCGGAAGtaaa                                                                                                                                                                                          >  2:303887/1‑227 (MQ=255)
                                     gCTGAGGGTGTTGGTGACCTCGAAGGGGTTTGCTGGCGTGCCTGGCGGTTCTTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAgc                                                                                                                                                                                         <  1:629302/228‑1 (MQ=255)
                                       tGACGCTGATGGTGTCCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAgcgc                                                                                                                                                                                       <  2:396752/228‑1 (MQ=255)
                                       tGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAgcgc                                                                                                                                                                                       <  2:372534/228‑1 (MQ=255)
                                       tGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAgcgc                                                                                                                                                                                       <  2:477776/228‑1 (MQ=255)
                                       tGACGCTGATGGGTACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTGTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAgcgc                                                                                                                                                                                       <  2:337456/228‑1 (MQ=255)
                                        tagcgtgggtgTGCCCTCGAAAGGGATTGCTGGCGTGCCGGGGGTGTCGTGTGTGGTGTTGCTGGCAACGCTGGGTAGCGTATGTATTCCGCTGGAAGGTCTGGCTGTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGACCACA‑AATTTGACCGTAAGAAAgcgc                                                                                                                                                                                       <  1:159055/217‑1 (MQ=255)
                                         aCGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGGGGTGGGTAATGCGCTGGCGGGGCTGGGGGT‑TGCCAAGGGGGGACACA‑AAATTGACCggaaggaaagggag                                                                                                                                                                                     >  1:540430/1‑222 (MQ=255)
                                             tGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGGGAT‑TGCCAAGTGGGAGCACA‑AATTTggcggtcatgaatacaggtcgt                                                                                                                                                                                 >  2:548240/1‑207 (MQ=255)
                                              gtggggtgccTGGAAAGGGATTTCTGGGGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTa                                                                                                                                                                                <  1:238905/220‑1 (MQ=255)
                                                  ttgtccccGAAAGGGTTGTGTGGGGTGCCGGGGGTGTGGTTTGTGGTTTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGATCGTAAGAAAGCGCTGGCTTATGa                                                                                                                                                                             <  1:103796/220‑1 (MQ=255)
                                                      ccTCGAAAGGGTTTGCGTGCGTGCGTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGTTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTg                                                                                                                                                                        <  1:303887/228‑1 (MQ=255)
                                                         cGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGGGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGGGGGGACACA‑AATTTGACCGTAAGAAAGGGCTGGCTTATGGGCGTGaaa                                                                                                                                                                     >  2:284980/1‑227 (MQ=255)
                                                                 aTTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGGGAT‑TGCCAAGGGGGGACACA‑AATTTGACCGTAAGAAAGCGCTGGGTTTTGAGAGTGAAGCGCTGGGc                                                                                                                                                             >  2:211638/1‑228 (MQ=255)
                                                                                  cGTTTCTTTTGTGGTGTTGCTGGCAACGCTGGGTTGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTg                                                                                                                                            <  1:155605/227‑1 (MQ=255)
                                                                                       cGTTTGTGGTGTTGGTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATc                                                                                                                                      <  2:47591/228‑1 (MQ=255)
                                                                                       cGTTTGTGGTGGTGCTGGGAACGTTGGGTAGCGGCGGTATTCCGCGGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGGGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATc                                                                                                                                      <  1:572535/228‑1 (MQ=255)
                                                                                        gTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATc                                                                                                                                      <  2:618741/227‑1 (MQ=255)
                                                                                            gggggggtcgttggaaCGCCGGGGAGCGTAGGTATTTCGCTTGAAGGTCTGGGGTTTATTGCTGGAGTTGACCGTATCCGCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCaataa                                                                                                                                 <  2:657299/214‑1 (MQ=255)
                                                                                                tgtTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGGGAT‑TGCCCAGTGGGGACACA‑AATTTGACCGGAATAAAGCGCTGGCTTCTAAGCGTGAGAGGCTGGGGAAATTTG‑ATAAAAATGGGAGTCaaaaatgta                                                                                                                             >  1:58108/1‑225 (MQ=255)
                                                                                                  ttGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGCCGGGGCTGGGGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTGAGGAAGCGCTGGGTGATGAACGTGAAGTGCCGGGCAAAATTG‑ATAAAACTGCGGATCAACAATTAAAg                                                                                                                           >  1:32237/1‑228 (MQ=255)
                                                                                                   tGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCCAGTGGGAACACA‑AAATTGACCGTGAAAAAGGGCTGGCTTTTTGTAGGGAAGGGCGTGGCAAATTTG‑GTAAAACTCCGGAGCAAAAATTAAAGa                                                                                                                          >  1:208333/1‑228 (MQ=255)
                                                                                                       ggCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGc                                                                                                                      <  2:329323/228‑1 (MQ=255)
                                                                                                       ggCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGGGGGTAATGCGCTGGCGGTGGTGGGGAT‑TGCCCAGGGGGCACACC‑AATTTGTCCGTAAGAAGCGGGTGGCTTCGTTGCGGGAAGGCGTGGGCAAATTTGTAAAAAACAG‑GGGTCAATAATTTAAAAGTGc                                                                                                                      >  1:30101/1‑228 (MQ=255)
                                                                                                              ttGGGTAGCGTGGGTTTTCCGCTGGAAGGTCGGGGGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGat                                                                                                               <  1:211638/227‑1 (MQ=255)
                                                                                                              cTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCTTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGat                                                                                                               <  2:50167/228‑1 (MQ=255)
                                                                                                                  gTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGGTGGGCAAATTTT‑GTAAAACTGCGGACCAACAATTTAAGATTGCCGGGGgtaaca                                                                                                           >  2:445384/1‑225 (MQ=255)
                                                                                                                     gCGTAGGTTTTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCAc                                                                                                         <  1:449248/227‑1 (MQ=255)
                                                                                                                      cGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGGGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGCCCGTAAGAAAGCGCTGGCTTATGAGAGGGACGTGCTGGGCAAATTTT‑ATAAAACTGCGGATCAATAACTGAAGGATGCCGGGGAATTACAccc                                                                                                       >  1:54240/1‑228 (MQ=255)
                                                                                                                          gttgttCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTAGCCTCGACATGGCGCGTACTGCTCTGAACGTGGTGGGTAATGCGCGGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATAGCCGGGGATATCCACCCGGCa                                                                                                   <  1:585305/224‑1 (MQ=255)
                                                                                                                          ggTTTTCCGGTGGAAGGTCTGGCGTTTATTGCTGGTGTTGGCCGTAACCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCgg                                                                                                     <  2:540430/226‑1 (MQ=255)
                                                                                                                           gTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGGGCTTGGGAT‑TGCCAAGTGGGAAAACA‑AATTTGACCGGAAGAAAGAGGTGGCTTTTGAGCGGGAAGGGGTGGGCAAAATTG‑ATAAAACTGCGGGTCAATAATTGAAAGTTTCCGGGGGTATCCCCCCGGCaa                                                                                                  >  1:416163/1‑228 (MQ=255)
                                                                                                                           gTATTCCGATGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGAAATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCa                                                                                                   <  2:229331/227‑1 (MQ=255)
                                                                                                                           gTATTCCGATGGAAGGTCTGGCCTTTATTGCTGGGGGTGACCGTATCCTGGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCa                                                                                                   <  2:419486/227‑1 (MQ=255)
                                                                                                                               tCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTTGGGGT‑TGCCCAGGGGGAAAAAA‑AATTTGACCGTAAGACAACGCTGGCCTATAAAGGGGAAGTGCAGGGGCAATTTG‑ATAAAACAGTCGTGCAATAATTGAAGAATTGCGGGGATAAACACCCGGCAAagtg                                                                                              >  1:214451/1‑228 (MQ=255)
                                                                                                                                ccGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTg                                                                                                                                                                   <  2:49605/159‑1 (MQ=255)
                                                                                                                                  gCTGGAAGGTCTGGTGTTTATTGGTGGTGTTGACCGTGTCCGCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAtgtgtg                                                                                            <  1:367307/227‑1 (MQ=255)
                                                                                                                                      gAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAAATGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGCCAATGTGGGaaga                                                                                        >  2:517355/1‑225 (MQ=255)
                                                                                                                                                  tttATTGTTGGTTTTGACCGTGGCCTCGTCATGGCGCGGACTGCGCTGATCGGGGGGGGGAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACg                                                                            <  1:49645/227‑1 (MQ=255)
                                                                                                                                                       tgttGGTGTTTGCCGTTTCCTCGGCATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGc                                                                       <  1:289757/224‑1 (MQ=255)
                                                                                                                                                            gttttGACCGTTTCCTCGGCATGGCGGTTACTGCGCTGAACGTGTGGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGt                                                                  <  1:284980/224‑1 (MQ=255)
                                                                                                                                                            gtgtTGACCCTTTCCTCGATATGGCGCGTACTGGGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATAGCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGt                                                                  <  1:548615/227‑1 (MQ=255)
                                                                                                                                                                     gTATCCTCGACATGGCGCGTACTGCGCTGAACGGGGTGGGTAATGCGCTGGCGGTGCTGGTTAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                        <  2:51026/228‑1 (MQ=255)
                                                                                                                                                                      tttCCTCGACTTGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                        <  1:281883/225‑1 (MQ=255)
                                                                                                                                                                      tagccTCGACATGGCGCGTACTGCGCAGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGGGAT‑TGCCAAGGGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCGCCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                        <  2:549526/224‑1 (MQ=255)
                                                                                                                                                                      tATCCTCGACATGGCGCGTACTGGGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                        <  2:626770/227‑1 (MQ=255)
                                                                                                                                                                      tATCCGCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGGAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGCGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                        <  1:31536/227‑1 (MQ=255)
                                                                                                                                                                       aTCCTCGACATGGCGCGTACTGCGCTGAACGTTGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCCGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAt                                                                                                 <  2:22343/185‑1 (MQ=255)
                                                                                                                                                                          ctccgaaTGGGGCGTACTGCGCTGAACGGTGTGGGTAATGCGCTGGGGGTGCTGGGTGT‑TGCCAAGTGGGAACAAA‑AATTTGACCGTAAGAAGGCGCTGGCTTATGAGCGTGAAGGGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACg                                                   <  1:612357/222‑1 (MQ=255)
                                                                                                                                                                            cGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATTAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAAGTGGGAATGCCTGATTCGACGCTTGGCGCGTCTTCTCAGGCCTACGcc                                                 >  2:131552/1‑228 (MQ=255)
                                                                                                                                                                               gaTGGGGCGTTGTGTGCTGAACCGGGTGGTTAATTCGTTGGCGGTGCTGGTGAT‑TTGCAAGTTGGAACAAA‑AATTTGACCGTAAGAAAGCGTGGGCTTAAGAGCGTGAAGGTGTGGGAAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGAATGCCGGGGATAGCCACCCGGCAATGTGTGAATGCCTGATGCGGCGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                              <  1:627801/227‑1 (MQ=255)
                                                                                                                                                                                  ggCGCGTACTGGGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                                            <  2:532872/227‑1 (MQ=255)
                                                                                                                                                                                           ttcgctgttctgtgtggGGAATGTGGTGGTGGTTGTGGGTAG‑TGCTAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGCGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGGCGGGGATATCAACTCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                  <  1:263419/215‑1 (MQ=37)
                                                                                                                                                                                            gctcTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTTAT‑TGCCAAGTGGGAACACA‑AATTTGACAGTAAGAGAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                 <  2:427418/225‑1 (MQ=255)
                                                                                                                                                                                             cgcTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGTCTTATGAGGGTGAAGTGCTGGGCAAATTTG‑GTAAAACTGCGGATCAATAATTGAAAATTGCCTGGGATATCCACCCCGCAATGTGCGAAAGTCTGGTGCGACGCTAGCCGCGTCCTATCAGGCCTACGCCAGACAGCGCACTAGCCt                                >  2:253510/1‑228 (MQ=255)
                                                                                                                                                                                             cgcTGAACGTGGTGGGTAATGCGCTGGCGGTGCCGGGGAT‑TGCCAAGGGGGAACACA‑AATTTGACCGGAAGAAAACGCTGGCTTATGGACGTGAAGTGCTGGGCAAATTTG‑GTAAAACTGCGGATCAATAACTGAAGAATGCCGGGGATATCCACACAGCAATGGGTAAAAGCCTGATAGGAACATGTCCCCGCCTTATCAGGCCTCGGCCAAACACCGCCAAAGccc                                >  1:484310/1‑227 (MQ=255)
                                                                                                                                                                                                         tggGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTg                     <  2:563983/227‑1 (MQ=255)
                                                                                                                                                                                                            gTAATGCGCTGGCGGAGCTGGGGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCCGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑GTAAAACTGAGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGGGGGAATGCCTGATGCGACGCTTGCCGCGCCTTATCAGGCCTACCCCCCACAGCGCAATAACCCGATTTAgggggaatt                 >  2:380277/1‑220 (MQ=255)
                                                                                                                                                                                                                tGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGAAATCCACCCGGCAATGTGTGAATGCCTGATTGGACGCTTGCCGCGTCTTTTCAGGCCTACGCCAGACAACGCAAAAGCCCGGTTTACGGGGATTTTTTa             >  2:170208/1‑228 (MQ=255)
                                                                                                                                                                                                                  cgcTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAgg           <  2:54240/228‑1 (MQ=255)
                                                                                                                                                                                                                      ggCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTTAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTTGAGATTTCCGGGGATATTCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGCCTTATCAGGCCTTCCCCATACCGCCCAATAGGCCGATTTTGCGGTGTTTTTTAGGTTgg       >  1:115430/1‑228 (MQ=255)
                                                                                                                                                                                                                       gCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTTTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCggg      >  2:444247/1‑227 (MQ=255)
                                                                                                                                                                                                                          gTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATaa   <  2:115430/228‑1 (MQ=255)
                                                                                                                                                                                                                           tGCTGGTGAT‑TTCCAAGTGGGAACGCA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAg  <  2:567909/228‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                
GCATTGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAT‑TGCCAAGTGGGAACACA‑AATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTG‑ATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG  >  NC_000913/4295474‑4295920

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