Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,295,700 T→G 20.8% 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,7000TG20.8% 162.1 / 4.7 74V407G (GTG→GGG) gltPglutamate/aspartate:proton symporter
Reads supporting (aligned to +/- strand):  ref base T (11/44);  new base G (17/2);  total (28/46)
Fisher's exact test for biased strand distribution p-value = 1.07e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG  >  NC_000913/4295475‑4295921
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                
gCATTGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTg                                                                                                                                                                                                                                 <  1:475059/227‑1 (MQ=255)
 gatttgccTTTCCATTTGGCAGGAAATCATTGTGGTTCTGACGCTGATGGTGTACTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGAt                                                                                                                                                                                                                               <  2:577937/222‑1 (MQ=255)
 cATTGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGGGAt                                                                                                                                                                                                                               >  1:108057/1‑228 (MQ=255)
   ttGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAACGTGTTTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGCCATGGCGCGTACTGCGCTGAACGTGGTGGGTATTGCGGTGGCGGTGCTGGTGGTTg                                                                                                                                                                                                                            >  1:585305/1‑228 (MQ=255)
   ttGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGGGCGGGGGATTg                                                                                                                                                                                                                            >  1:194733/1‑228 (MQ=255)
           tCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCCGGGGATTGCCaaagtg                                                                                                                                                                                                                    >  1:67554/1‑224 (MQ=255)
           tCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGGATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGGCCTGGGGATTGCCAAggtg                                                                                                                                                                                                                    >  1:577937/1‑225 (MQ=255)
                cTGGCAGGAAATCATTCTGGTTCTGACGCTGGTGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAac                                                                                                                                                                                                                <  1:186552/227‑1 (MQ=255)
                     aGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGCATCCTCCACCTGGCGCGACCTCGGCTGCACGGGGTGGGTAGTCGGCTGCCGGGGAGGGGGATACGACAAGTGGAACACaaaa                                                                                                                                                                                                         >  2:85010/1‑227 (MQ=255)
                       gAAATCTTTCTGGGTCTGACGGTTTGTGTGACCTCGAAAGGGATTGGTGGCGTGCCGGGGGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAAttt                                                                                                                                                                                                       <  2:67554/228‑1 (MQ=255)
                         aaTCTTTCTGGGGTGGACGTGGATGGTGACTGCGAAGGGGAGTGCGTGCGTGGCTGGCGTGTCGTTTGTGGTGTTGCTGGTAACGTTGGGTAGCGTAGGTATCCCGCTGGAAGGTGTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCGGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGa                                                                                                                                                                                                     <  2:394700/228‑1 (MQ=255)
                         aaTCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGGGAGTGGCAATGGGGGACACAAAtttta                                                                                                                                                                                                     >  1:548615/1‑226 (MQ=255)
                                   gTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGGGCTGTGGATTGCCAAGTGGGAAACACAATTTGACCGTAAGaaa                                                                                                                                                                                           >  1:155605/1‑228 (MQ=255)
                                    tGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACCGCGCTGAACGTGGTGGGAATTGCGCTGGCGGGGCGGGGGGTTTCCAAGGGGGACACCAAAATTTACCGGAAGtaaa                                                                                                                                                                                          >  1:303887/1‑227 (MQ=255)
                                     gCTGAGGGTGTTGGTGACCTCGAAGGGGTTTGCTGGCGTGCCTGGCGGTTCTTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAgc                                                                                                                                                                                         <  2:629302/228‑1 (MQ=255)
                                       tGACGCTGATGGTGTCCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAgcgc                                                                                                                                                                                       <  1:396752/228‑1 (MQ=255)
                                       tGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAgcgc                                                                                                                                                                                       <  1:477776/228‑1 (MQ=255)
                                       tGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAgcgc                                                                                                                                                                                       <  1:372534/228‑1 (MQ=255)
                                       tGACGCTGATGGGTACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTGTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAgcgc                                                                                                                                                                                       <  1:337456/228‑1 (MQ=255)
                                        tagcgtgggtgTGCCCTCGAAAGGGATTGCTGGCGTGCCGGGGGTGTCGTGTGTGGTGTTGCTGGCAACGCTGGGTAGCGTATGTATTCCGCTGGAAGGTCTGGCTGTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGACCACAAATTTGACCGTAAGAAAgcgc                                                                                                                                                                                       <  2:159055/217‑1 (MQ=255)
                                         aCGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGGGGTGGGTAATGCGCTGGCGGGGCTGGGGGTTGCCAAGGGGGGACACAAAATTGACCggaaggaaagggag                                                                                                                                                                                     >  2:540430/1‑222 (MQ=255)
                                              gtggggtgccTGGAAAGGGATTTCTGGGGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTa                                                                                                                                                                                <  2:238905/220‑1 (MQ=255)
                                                  ttgtccccGAAAGGGTTGTGTGGGGTGCCGGGGGTGTGGTTTGTGGTTTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGATCGTAAGAAAGCGCTGGCTTATGa                                                                                                                                                                             <  2:103796/220‑1 (MQ=255)
                                                      ccTCGAAAGGGTTTGCGTGCGTGCGTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGTTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTg                                                                                                                                                                        <  2:303887/228‑1 (MQ=255)
                                                         cGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGGGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGGGGGGACACAAATTTGACCGTAAGAAAGGGCTGGCTTATGGGCGTGaaa                                                                                                                                                                     >  1:284980/1‑227 (MQ=255)
                                                                 aTTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGGGATTGCCAAGGGGGGACACAAATTTGACCGTAAGAAAGCGCTGGGTTTTGAGAGTGAAGCGCTGGGc                                                                                                                                                             >  1:211638/1‑228 (MQ=255)
                                                                                  cGTTTCTTTTGTGGTGTTGCTGGCAACGCTGGGTTGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTg                                                                                                                                            <  2:155605/227‑1 (MQ=255)
                                                                                       cGTTTGTGGTGTTGGTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATc                                                                                                                                      <  1:47591/228‑1 (MQ=255)
                                                                                       cGTTTGTGGTGGTGCTGGGAACGTTGGGTAGCGGCGGTATTCCGCGGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGGGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATc                                                                                                                                      <  2:572535/228‑1 (MQ=255)
                                                                                        gTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATc                                                                                                                                      <  1:618741/227‑1 (MQ=255)
                                                                                                tgtTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGGGATTGCCCAGTGGGGACACAAATTTGACCGGAATAAAGCGCTGGCTTCTAAGCGTGAGAGGCTGGGGAAATTTGATAAAAATGGGAGTCaaaaatgta                                                                                                                             >  2:58108/1‑225 (MQ=255)
                                                                                                  ttGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGCCGGGGCTGGGGATTGCCAAGTGGGAACACAAATTTGACCGTGAGGAAGCGCTGGGTGATGAACGTGAAGTGCCGGGCAAAATTGATAAAACTGCGGATCAACAATTAAAg                                                                                                                           >  2:32237/1‑228 (MQ=255)
                                                                                                   tGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCCAGTGGGAACACAAAATTGACCGTGAAAAAGGGCTGGCTTTTTGTAGGGAAGGGCGTGGCAAATTTGGTAAAACTCCGGAGCAAAAATTAAAGa                                                                                                                          >  2:208333/1‑228 (MQ=255)
                                                                                                       ggCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGc                                                                                                                      <  1:329323/228‑1 (MQ=255)
                                                                                                       ggCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGGGGGTAATGCGCTGGCGGTGGTGGGGATTGCCCAGGGGGCACACCAATTTGTCCGTAAGAAGCGGGTGGCTTCGTTGCGGGAAGGCGTGGGCAAATTTGTAAAAAACAGGGGTCAATAATTTAAAAGTGc                                                                                                                      >  2:30101/1‑228 (MQ=255)
                                                                                                              ttGGGTAGCGTGGGTTTTCCGCTGGAAGGTCGGGGGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGat                                                                                                               <  2:211638/227‑1 (MQ=255)
                                                                                                              cTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCTTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGat                                                                                                               <  1:50167/228‑1 (MQ=255)
                                                                                                                  gTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGGTGGGCAAATTTTGTAAAACTGCGGACCAACAATTTAAGATTGCCGGGGgtaaca                                                                                                           >  1:445384/1‑225 (MQ=255)
                                                                                                                     gCGTAGGTTTTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCAc                                                                                                         <  2:449248/227‑1 (MQ=255)
                                                                                                                      cGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGGGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGCCCGTAAGAAAGCGCTGGCTTATGAGAGGGACGTGCTGGGCAAATTTTATAAAACTGCGGATCAATAACTGAAGGATGCCGGGGAATTACAccc                                                                                                       >  2:54240/1‑228 (MQ=255)
                                                                                                                          gttgttCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTAGCCTCGACATGGCGCGTACTGCTCTGAACGTGGTGGGTAATGCGCGGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATAGCCGGGGATATCCACCCGGCa                                                                                                   <  2:585305/224‑1 (MQ=255)
                                                                                                                          ggTTTTCCGGTGGAAGGTCTGGCGTTTATTGCTGGTGTTGGCCGTAACCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCgg                                                                                                     <  1:540430/226‑1 (MQ=255)
                                                                                                                           gTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGGGCTTGGGATTGCCAAGTGGGAAAACAAATTTGACCGGAAGAAAGAGGTGGCTTTTGAGCGGGAAGGGGTGGGCAAAATTGATAAAACTGCGGGTCAATAATTGAAAGTTTCCGGGGGTATCCCCCCGGCaa                                                                                                  >  2:416163/1‑228 (MQ=255)
                                                                                                                           gTATTCCGATGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGAAATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCa                                                                                                   <  1:229331/227‑1 (MQ=255)
                                                                                                                           gTATTCCGATGGAAGGTCTGGCCTTTATTGCTGGGGGTGACCGTATCCTGGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCa                                                                                                   <  1:419486/227‑1 (MQ=255)
                                                                                                                               tCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTTGGGGTTGCCCAGGGGGAAAAAAAATTTGACCGTAAGACAACGCTGGCCTATAAAGGGGAAGTGCAGGGGCAATTTGATAAAACAGTCGTGCAATAATTGAAGAATTGCGGGGATAAACACCCGGCAAagtg                                                                                              >  2:214451/1‑228 (MQ=255)
                                                                                                                                ccGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTg                                                                                                                                                                   <  1:49605/159‑1 (MQ=255)
                                                                                                                                  gCTGGAAGGTCTGGTGTTTATTGGTGGTGTTGACCGTGTCCGCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAtgtgtg                                                                                            <  2:367307/227‑1 (MQ=255)
                                                                                                                                      gAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAAATGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGCCAATGTGGGaaga                                                                                        >  1:517355/1‑225 (MQ=255)
                                                                                                                                                  tttATTGTTGGTTTTGACCGTGGCCTCGTCATGGCGCGGACTGCGCTGATCGGGGGGGGGAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACg                                                                            <  2:49645/227‑1 (MQ=255)
                                                                                                                                                       tgttGGTGTTTGCCGTTTCCTCGGCATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGc                                                                       <  2:289757/224‑1 (MQ=255)
                                                                                                                                                            gttttGACCGTTTCCTCGGCATGGCGGTTACTGCGCTGAACGTGTGGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGt                                                                  <  2:284980/224‑1 (MQ=255)
                                                                                                                                                            gtgtTGACCCTTTCCTCGATATGGCGCGTACTGGGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATAGCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGt                                                                  <  2:548615/227‑1 (MQ=255)
                                                                                                                                                                     gTATCCTCGACATGGCGCGTACTGCGCTGAACGGGGTGGGTAATGCGCTGGCGGTGCTGGTTATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                        <  1:51026/228‑1 (MQ=255)
                                                                                                                                                                      tttCCTCGACTTGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                        <  2:281883/225‑1 (MQ=255)
                                                                                                                                                                      tagccTCGACATGGCGCGTACTGCGCAGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGGGATTGCCAAGGGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCGCCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                        <  1:549526/224‑1 (MQ=255)
                                                                                                                                                                      tATCCTCGACATGGCGCGTACTGGGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                        <  1:626770/227‑1 (MQ=255)
                                                                                                                                                                      tATCCGCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGGAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGCGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                        <  2:31536/227‑1 (MQ=255)
                                                                                                                                                                       aTCCTCGACATGGCGCGTACTGCGCTGAACGTTGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCCGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAt                                                                                                 <  1:22343/185‑1 (MQ=255)
                                                                                                                                                                          ctccgaaTGGGGCGTACTGCGCTGAACGGTGTGGGTAATGCGCTGGGGGTGCTGGGTGTTGCCAAGTGGGAACAAAAATTTGACCGTAAGAAGGCGCTGGCTTATGAGCGTGAAGGGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACg                                                   <  2:612357/222‑1 (MQ=255)
                                                                                                                                                                            cGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATTAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAAGTGGGAATGCCTGATTCGACGCTTGGCGCGTCTTCTCAGGCCTACGcc                                                 >  1:131552/1‑228 (MQ=255)
                                                                                                                                                                               gaTGGGGCGTTGTGTGCTGAACCGGGTGGTTAATTCGTTGGCGGTGCTGGTGATTTGCAAGTTGGAACAAAAATTTGACCGTAAGAAAGCGTGGGCTTAAGAGCGTGAAGGTGTGGGAAAATTTGATAAAACTGCGGATCAATAATTGAAGAATGCCGGGGATAGCCACCCGGCAATGTGTGAATGCCTGATGCGGCGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                              <  2:627801/227‑1 (MQ=255)
                                                                                                                                                                                  ggCGCGTACTGGGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                                            <  1:532872/227‑1 (MQ=255)
                                                                                                                                                                                            gctcTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTTATTGCCAAGTGGGAACACAAATTTGACAGTAAGAGAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                 <  1:427418/225‑1 (MQ=255)
                                                                                                                                                                                             cgcTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGTCTTATGAGGGTGAAGTGCTGGGCAAATTTGGTAAAACTGCGGATCAATAATTGAAAATTGCCTGGGATATCCACCCCGCAATGTGCGAAAGTCTGGTGCGACGCTAGCCGCGTCCTATCAGGCCTACGCCAGACAGCGCACTAGCCt                                >  1:253510/1‑228 (MQ=255)
                                                                                                                                                                                             cgcTGAACGTGGTGGGTAATGCGCTGGCGGTGCCGGGGATTGCCAAGGGGGAACACAAATTTGACCGGAAGAAAACGCTGGCTTATGGACGTGAAGTGCTGGGCAAATTTGGTAAAACTGCGGATCAATAACTGAAGAATGCCGGGGATATCCACACAGCAATGGGTAAAAGCCTGATAGGAACATGTCCCCGCCTTATCAGGCCTCGGCCAAACACCGCCAAAGccc                                >  2:484310/1‑227 (MQ=255)
                                                                                                                                                                                                         tggGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTg                     <  1:563983/227‑1 (MQ=255)
                                                                                                                                                                                                            gTAATGCGCTGGCGGAGCTGGGGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCCGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGGTAAAACTGAGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGGGGGAATGCCTGATGCGACGCTTGCCGCGCCTTATCAGGCCTACCCCCCACAGCGCAATAACCCGATTTAgggggaatt                 >  1:380277/1‑220 (MQ=255)
                                                                                                                                                                                                                tGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGAAATCCACCCGGCAATGTGTGAATGCCTGATTGGACGCTTGCCGCGTCTTTTCAGGCCTACGCCAGACAACGCAAAAGCCCGGTTTACGGGGATTTTTTa             >  1:170208/1‑228 (MQ=255)
                                                                                                                                                                                                                  cgcTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAgg           <  1:54240/228‑1 (MQ=255)
                                                                                                                                                                                                                      ggCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTTAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTTGAGATTTCCGGGGATATTCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGCCTTATCAGGCCTTCCCCATACCGCCCAATAGGCCGATTTTGCGGTGTTTTTTAGGTTgg       >  2:115430/1‑228 (MQ=255)
                                                                                                                                                                                                                       gCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTTTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCggg      >  1:444247/1‑227 (MQ=255)
                                                                                                                                                                                                                          gTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATaa   <  1:115430/228‑1 (MQ=255)
                                                                                                                                                                                                                           tGCTGGTGATTTCCAAGTGGGAACGCAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAg  <  1:567909/228‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                
GCATTGACCTGTCCATCTGGCAGGAAATCATTCTGGTGCTGACGCTGATGGTGACCTCGAAAGGGATTGCTGGCGTGCCTGGCGTGTCGTTTGTGGTGTTGCTGGCAACGCTGGGTAGCGTAGGTATTCCGCTGGAAGGTCTGGCGTTTATTGCTGGTGTTGACCGTATCCTCGACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG  >  NC_000913/4295475‑4295921

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: