Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 86,658 A→G P343P (CCA→CCG ilvI → acetolactate synthase/acetohydroxybutanoate synthase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091386,6580AG94.0% 145.3 / ‑0.9 50P343P (CCA→CCGilvIacetolactate synthase/acetohydroxybutanoate synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (23/24);  total (24/26)
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.53e-01

GAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCACTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCGGGCAGCACCAGATGTTTGCTGCACTTTATTATCCATTCGACAAACCGCGTCGCTGGATCAATTCCGGTGGCCTCGGCACGATGGGTTT  >  NC_000913/86422‑86903
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       
gAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAACTGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCcg                                                                                                                                                                                                                                           >  1:422759/1‑251 (MQ=255)
 aaGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCcgc                                                                                                                                                                                                                                          >  1:280088/1‑251 (MQ=255)
    ccAATTTGACGAGGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCGTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGAc                                                                                                                                                                                                                                       <  2:13802/251‑1 (MQ=255)
       atatGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTCAAACCGTGACTGCGGATATCCCGATTGTGTGGGATGCTCGCCAGGTCCTCGAACAACTTCTTGCACTCTGGTCGCAAGAATCCGCCCATCAACCGCTGAATGAGCTCCGCGCCtgg                                                                                                                                                                                                                                    >  2:198421/1‑251 (MQ=255)
         atGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATAGCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACtggt                                                                                                                                                                                                                                   <  2:43713/250‑1 (MQ=255)
                   aTAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAAtt                                                                                                                                                                                                                        >  1:260936/1‑251 (MQ=255)
                         cgGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTGCTGCATATCGATATGGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCAGCAACCGCGGGATGAGATCCGCGACTGGTGGCATCAAATTGAACAg                                                                                                                                                                                                                  <  1:21975/251‑1 (MQ=255)
                               tgATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCACAGAACTGCCCAAATGCCACTGTTCTGCATATCGATATTGCTCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCCACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGgcgc                                                                                                                                                                                                            >  1:34849/1‑251 (MQ=255)
                                    ttCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCACATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCACGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCg                                                                                                                                                                                                       >  1:442950/1‑251 (MQ=255)
                                      cGCCGTCGGGGGACGATTTGATGACCGATCGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTTCATATCGATATTGATCCTACTTCCATTGCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGt                                                                                                                                                                                                      <  2:205669/250‑1 (MQ=255)
                                       gCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTc                                                                                                                                                                                                     >  2:44042/1‑250 (MQ=255)
                                                         gatgaCCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGacac                                                                                                                                                                                  >  1:277764/1‑251 (MQ=255)
                                                               cGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTcac                                                                                                                                                                              >  1:348565/1‑249 (MQ=255)
                                                                       gaACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCCATCGTCAGTGCCTGAAATATGACACTCACAGTaaaa                                                                                                                                                                       >  1:13513/1‑248 (MQ=255)
                                                                          cAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGAtt                                                                                                                                                                  >  1:31068/1‑250 (MQ=255)
                                                                                ggCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGGCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCg                                                                                                                                                            <  1:413800/250‑1 (MQ=255)
                                                                                   aaaGTACTGCCCAAATGCCACTGTTCTGCATATAGATATTGATCCTACTTCGATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAgg                                                                                                                                                        >  1:158500/1‑251 (MQ=255)
                                                                                         cTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTgg                                                                                                                                                  >  2:173135/1‑250 (MQ=255)
                                                                                             caaaaTGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGGGAAAAGATTAAACCGCAGGCGGTGATCga                                                                                                                                              <  1:198421/249‑1 (MQ=255)
                                                                                              cAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCg                                                                                                                                               >  1:467818/1‑249 (MQ=255)
                                                                                               aaaTGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCAACGTGAAAAGATTAAACCGCAGGCGGAGATCggg                                                                                                                                             >  1:280581/1‑248 (MQ=255)
                                                                                                 aTGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAGGATTAAACCGCAGGCGGGGAACgagccc                                                                                                                                          >  2:94924/1‑248 (MQ=255)
                                                                                                            cTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTTACTGCGCATATCCCGATTGTGCGGGATTCTCGCCAGGTCCTCGACCAAATGCTTGAACTCTTGTCGCAAGAATCAGCCAATCAACCGCTGGATGAGATCCGCGACTGGTGTCAGCACATTGAACAGAGGCGCGCTCGTCAGTGCCTGAAATCTGACACTCACAGTGAAAAGCTTAAACCGCGGGCGGGGACCGAGACACATTGGAGGtt                                                                                                                               >  2:285250/1‑251 (MQ=255)
                                                                                                              gCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGtt                                                                                                                               <  2:204239/249‑1 (MQ=255)
                                                                                                                     gctatTGATCCTTCTTACAGTTCTAAAACCGTTGCTGTGGATAGCCCGATTGTGGGGGATGCTCGCCAGGTCTTCGAACAAATGCTTGAACTCTTTTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAggg                                                                                                                      <  1:428303/249‑1 (MQ=255)
                                                                                                                          tGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAACTTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCCCAGTGAAAAGATTAAACCGCAGGCGGAGATCGAGACTCTTTGGCgggtgtggcaggggga                                                                                                                   >  2:213770/1‑235 (MQ=255)
                                                                                                                            aTCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGCAGGGAGACGCt                                                                                                               >  1:235284/1‑251 (MQ=255)
                                                                                                                              ccTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACCCTTTGGCGGTTGACGGAAGGAGACGCttt                                                                                                             >  1:122509/1‑250 (MQ=255)
                                                                                                                                  cTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATTACACTCACAGTGAAAAGATTAAACCGCACGCGGTGATCGAAAATCTTTGGCgggttgggaagagggcagcttgtc                                                                                                           >  2:144853/1‑227 (MQ=255)
                                                                                                                                    tCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCACTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTg                                                                                                         <  1:244905/249‑1 (MQ=255)
                                                                                                                                               aaCCGTGACTGCGGATATCCCGATTGTGGGGGATGCCCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCg                                                                                           <  1:83561/251‑1 (MQ=255)
                                                                                                                                                 ccgtgcgtGCGGATATCCCGATTGGGGGGGATGCTCGACAGGGCCTCGAACAAAGGCTGGAACTCTTGGCTCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCGGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCggt                                                                                         <  2:280088/244‑2 (MQ=255)
                                                                                                                                                   gTGACTGCGGATAGCACGATTGTGGGGGATGCTCGCCATGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCACATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCGGgca                                                                                       <  2:378430/251‑1 (MQ=255)
                                                                                                                                                      gttGCGGATATCCCGAGTGTGGGGGATGCACGTCAGGTCCTAGAACAAAGGCTTTAACTCCTGTCACAAGAATCAGCCCAACAACCACGGGATGAGATCCGCGACTGGTGGCAGCAAATTTAACAGGGGGTCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGAGCGAGACTCTTTTGCGGTTGACGAAGGGAGACGCTTACGTGACGGCCGATGTCGGgcagca                                                                                    <  2:413124/249‑1 (MQ=255)
                                                                                                                                                          cGGATAGCCCGATTGATGGGGATGCTGGCCAGGTCCCCGAACAAATGATTGAACTCTGGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGGCTGGGGGCAGCAAATTGAACAGTGGCGCTCTCGTCAGTGCCTGAAATATGACCCTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGGGACGTCGGATGACGGGCAGCACCa                                                                                 <  2:315896/250‑1 (MQ=255)
                                                                                                                                                                 cccGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCGGGCAGCACCAGATGttt                                                                          <  1:213770/250‑1 (MQ=255)
                                                                                                                                                                      ttGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCGTTGGCGGTTGACGAAGGGAGACGCGTACGTGACGTCCGATGTCGGGCAGCACCAGATGTTtgctgc                                                                     <  2:306903/250‑1 (MQ=255)
                                                                                                                                                                               aTGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGGGGCGCTTACGTGACGGCCGATGGCGGGCAGCACCAGATGTTTGCTGGACTTtattat                                                           >  1:108960/1‑251 (MQ=255)
                                                                                                                                                                                        agggccTCGAACAAATGCTTGAACGCTTGGCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACATTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCGGTCAGCACCAGATGTTTGCTGCACTTTATTATCCATTCGAc                                                 <  1:456146/247‑1 (MQ=255)
                                                                                                                                                                                            ccTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCGGGCAGCACCAGATGTTTGCTGCACTTTATTATCCATTCGACAAAc                                             <  2:31068/251‑1 (MQ=255)
                                                                                                                                                                                                   cAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACg                                                                                                           <  1:289444/184‑1 (MQ=255)
                                                                                                                                                                                                           tGAACTCTTGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCGGGCAGCACCAGATGTTTGCTGCACTTTATTATCCATTCGACAAACCGCGTCGCTGGAt                                <  1:93251/249‑1 (MQ=255)
                                                                                                                                                                                                                   tGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCGGGCAGCACCAGATGTTTGCTGCACTTTATTATCCATTCGACAAACCGCGTCGCTGGATCAATTcc                         <  1:173135/248‑1 (MQ=255)
                                                                                                                                                                                                                   tGTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCGGGCAGCACCAGATGTTTGCTGCACTTTATTATCCATTCGACAAACCGCGTCGCTGGATCAATTCCg                        <  1:180215/249‑1 (MQ=255)
                                                                                                                                                                                                                   gtcagccaGAATCCGCTCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCTCGCTCGTCAGTGCCTGAAATATGACCCTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGATACGCTTACGTGACGTCCGATGTCGGGCAGCACCAGATGTTTGCTGCACTTTATTATCCATTCGACAAACCGGGTCGCTGGATCAAGTCCg                        <  2:78712/242‑1 (MQ=255)
                                                                                                                                                                                                                    gTCGCAAGAATCCGCCCATCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGGGGCGCTTACGTGACGTCCGATGTCGGGCAGCACCAGATGTTTTCTGCACTTTATTTTTCCTTTCCCAAACCGCGTCGTCGGGTCAATTCCgggg                     >  1:365276/1‑249 (MQ=255)
                                                                                                                                                                                                                        caaaaaTCCGCCCACCAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCGGGCAGCACCAGATGTTTGCTGCACTTTATTATCCATTCGACAAACCGCGGCGCTGGATCAATTCCGGTGGCCt                 <  2:297646/247‑1 (MQ=255)
                                                                                                                                                                                                                         aaGAATCCGCCCATCAACCACTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGGGACTCTTTGGCGGTTGACGCAGGGAGGAGCTTACGTGGCGGCCCGTGTCCGGGCGAACCCGATGGTTTCTTCACTTTATTATACATTCGTCAAACCGCGGCGCTGGATAAATTTCGGAGGcacc                >  2:178217/1‑248 (MQ=255)
                                                                                                                                                                                                                                       cAACCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGGGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCGGGCAGCACCAGATGTTTGCTGCACTTTATTATCCATTCGACAAACCGCGTCGCTGGATCAATTCCGGTGGCCTCGGCACGATGGGtt   <  2:467818/250‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaCCGCTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCCCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGAGTAAACCGCAGGCGGTGATCGAGACTCTTTGGCTGTTGACGAAGGGATACGCTTACGTGACGTCCGATGTCGGGCTGCACCAGATGTTTGCTGCACTTTATTATCCATGCGACAAACCGCGTCGCTGGATCAATTCCGGTGGCCTCGTCACGATGGGttt  <  1:50787/250‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       
GAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAATCCGCCCATCAACCACTGGATGAGATCCGCGACTGGTGGCAGCAAATTGAACAGTGGCGCGCTCGTCAGTGCCTGAAATATGACACTCACAGTGAAAAGATTAAACCGCAGGCGGTGATCGAGACTCTTTGGCGGTTGACGAAGGGAGACGCTTACGTGACGTCCGATGTCGGGCAGCACCAGATGTTTGCTGCACTTTATTATCCATTCGACAAACCGCGTCGCTGGATCAATTCCGGTGGCCTCGGCACGATGGGTTT  >  NC_000913/86422‑86903

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