Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,951,369 C→T intergenic (+812/‑832) ilvX → / → ilvM uncharacterized protein IlvX/acetolactate synthase II subunit IlvM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,3690CT98.7% 243.4 / NA 75intergenic (+812/‑832)ilvX/ilvMuncharacterized protein IlvX/acetolactate synthase II subunit IlvM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (35/39);  minor base G (0/1);  total (35/40)
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.94e-01

GCGCCAGATGCTGGCAAAAGCGCAAAAACCGATGCTGTACGTTGGCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAC  >  NC_000913/3951126‑3951607
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
gcgcCAGATGCTGGCAAAAGCGCAAAAACCGATGCTGTACGGTGGCGTTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGtt                                                                                                                                                                                                                                           <  2:378857/251‑1 (MQ=255)
    cAGATGCTGGCAAAAGCGCAAAAACCGATGCTGTACGTTGGCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGGGCAGGAGTGCGACTTGCTGAGCGCCGGGGGTGCACGTTTtgg                                                                                                                                                                                                                                       >  2:28529/1‑250 (MQ=255)
     aGATGCTGGCAAAAGCGCAAAACCCGATGCTNTCCGGTGGTGGTGGCGTGGGTTTGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTAGCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAGCTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTtgat                                                                                                                                                                                                                                      <  2:212880/251‑1 (MQ=255)
               aaaaGCGCAAAAACCGATGCTGTACGTTGGCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTAGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTg                                                                                                                                                                                                                             >  2:379149/1‑250 (MQ=255)
                    cgcAAAACCCGATGCTGTACGTTGGCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGc                                                                                                                                                                                                                       <  1:272648/251‑1 (MQ=255)
                      cAAAAACCGATGCTGTAAGTTGGCGGTGGCGCGGGTATGGCGCAGGCAGTTCCGGCTTTGCGAGAATTTCTCGCGGCCACAAAAATGCATGCCACCTGTACGCTGAAAGGGGTGGGCGCAGTAGAAGCAACTTAACCGGCATATCTGCGCATGCTGGGAATGCAGGGCACCACAGCCGCGTACTTCGCGGTCCAGGAGTGCGACGTGCTGATCGCCGCGAGTACAAGTTTTGATGGTCGGGTGAcggggg                                                                                                                                                                                                                      >  2:307395/1‑245 (MQ=255)
                           aCCGATGCTGTTCGTTGGCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGa                                                                                                                                                                                                                <  2:340038/251‑1 (MQ=255)
                             cGATGCTGTACGTTGGCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGGTTTGATGACCGGGTGACCGGCAAACTGAac                                                                                                                                                                                                              <  1:77132/251‑1 (MQ=255)
                               aTGCTGTACGTTGGCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAacac                                                                                                                                                                                                            <  1:401220/251‑1 (MQ=255)
                                    gTACGTTGGGGGGGGCGGGGGTAAGGCGCAGGCAGTTCCGGCTGGGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACGTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTcg                                                                                                                                                                                                       <  2:44407/251‑1 (MQ=255)
                                    gTACGTTGGCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCtt                                                                                                                                                                                                         <  1:450032/249‑1 (MQ=255)
                                    gTACGTTGGCGGGGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAATCAGATTATCCGTACTATCTTGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTcg                                                                                                                                                                                                       <  2:183274/251‑1 (MQ=255)
                                      aCGTTGGCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTcgca                                                                                                                                                                                                     <  1:94851/251‑2 (MQ=255)
                                           ggCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGGTGACCGGGTGACCGGCAAACTGAACACCTTCGcaccac                                                                                                                                                                                                 >  2:279046/1‑250 (MQ=255)
                                            gCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCacac                                                                                                                                                                                               <  1:45853/251‑1 (MQ=255)
                                              ggTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGc                                                                                                                                                                                             >  1:241100/1‑251 (MQ=255)
                                                    gTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAgtgt                                                                                                                                                                                       >  1:112422/1‑251 (MQ=255)
                                                                  ggcagTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGat                                                                                                                                                                          >  2:484942/1‑250 (MQ=255)
                                                                   gcagTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGatat                                                                                                                                                                        >  2:485957/1‑251 (MQ=255)
                                                                         cctccTTTGCGTGAATTGCTCGATGCCACAAAAATGCCTGACACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCAGGCCGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACCCCTTCGCACCTCACGCCAGTGTTATCCATATGGATATCGAcc                                                                                                                                                                   <  2:383824/246‑1 (MQ=255)
                                                                           ggCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCgg                                                                                                                                                                <  1:196703/251‑1 (MQ=255)
                                                                               ttGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAACCTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGc                                                                                                                                                               >  1:177726/1‑248 (MQ=255)
                                                                               ttGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAg                                                                                                                                                             <  1:43791/250‑1 (MQ=255)
                                                                                tGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGa                                                                                                                                                            <  1:379149/250‑1 (MQ=255)
                                                                                 gCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGaaa                                                                                                                                                          <  1:411720/251‑1 (MQ=255)
                                                                                          tctcGCTGCCTCAAAAACGCCTGCCACCTGTACGCTGAAAGGGCTGTGCGCAGTAGAAGCAGATTACCCGTTCTATCTGGGCATGCTGTGAATGCATGGCCCCAAAGCGGCGAACTTCACGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACGGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAg                                                                                                                                                  <  1:154741/250‑1 (MQ=255)
                                                                                                     caAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGca                                                                                                                                      <  1:313265/251‑1 (MQ=255)
                                                                                                       aaaaTGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGCTGACCGGGGGACCGGCAAACTGAACACCTTcgcac                                                                                                                                                                                                    >  2:170569/1‑185 (MQ=255)
                                                                                                       aaaaTGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTcgcac                                                                                                                                                                                                    <  1:170569/187‑3 (MQ=255)
                                                                                                       aaaaTGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGcac                                                                                                                                     >  1:149572/1‑250 (MQ=255)
                                                                                                          aTGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCCTGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACAtgt                                                                                                                                 >  1:55211/1‑251 (MQ=255)
                                                                                                          aTGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACAtgt                                                                                                                                 <  1:279046/251‑1 (MQ=255)
                                                                                                           tgcctgccACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACAtgtg                                                                                                                                >  2:50539/1‑251 (MQ=255)
                                                                                                           tgcctgccACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGACAGGCACAtgtg                                                                                                                                >  2:252089/1‑251 (MQ=255)
                                                                                                            tcctgccACCTGTACGGTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCGGGGCATGCTGGGAAGGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGGACGTGTTGATGACCGGGTGACCGGCAAACTGAACACCTGCGCACCACACGCCAGTTGTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACAtgtg                                                                                                                                <  2:34360/249‑1 (MQ=255)
                                                                                                            gcctgccACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACAtgtg                                                                                                                                <  2:307265/250‑1 (MQ=255)
                                                                                                                       gTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAgg                                                                                                                     <  1:170026/250‑1 (MQ=255)
                                                                                                                       gTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCACGTGTTATCCATATGGATATCGACCCGGCCGAAATGAACANGCTGCGTCAGGCACATGTGGCTTTACAAggg                                                                                                                    >  2:229783/1‑250 (MQ=255)
                                                                                                                         aCGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGCCCGGGTGACCGGCCAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCCGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAggcgc                                                                                                                  >  2:205041/1‑248 (MQ=255)
                                                                                                                              gAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGGTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTa                                                                                                              >  1:266724/1‑250 (MQ=255)
                                                                                                                                 aGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATg                                                                                                          <  1:378443/251‑1 (MQ=255)
                                                                                                                                   gggtGGGCGGAGTAGAAGCAGATTATCCGTACTAGCTGGGAATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGct                                                                                                        <  2:32437/248‑1 (MQ=255)
                                                                                                                                      tGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCCTGGCACCAAAGCGGAGAACTTCGCGGTGCAGGCGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAAATGAACACCTTCGCACCACACGCCCGTGATATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCCCATGTGGCATTACAAGGTGATTTAACTGCTCTg                                                                                                     >  1:293924/1‑251 (MQ=255)
                                                                                                                                          cgcAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCAGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGACATGAACAAGCTGCGTCAGGCACATGTNGCATTACAAGGGGATTTAAATGCTCTGTTAc                                                                                                 >  2:246174/1‑251 (MQ=255)
                                                                                                                                           gcAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTAcc                                                                                                >  2:256352/1‑251 (MQ=255)
                                                                                                                                            cAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCa                                                                                               <  1:210465/251‑1 (MQ=255)
                                                                                                                                                  aaGCATATTATCCGTACTAGCTGGGCATGCTGGGAAGGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGACAGTGTTATCCATATGGATATCGACCCGGCAGAAAGGAACAAGCTGTGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTa                                                                                         <  2:338133/251‑1 (MQ=255)
                                                                                                                                                     cAGATTATCCGTACTATCTGGGCATGATGGGAATGCATGGCACCAAAGCGGCGCACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACAACACGCCAGTGTTATCCATATGGATATCGACCCGGAAGACATGAACAAGCTGCGTCCGGCACAAGTGGCATTCCCAGGGGATTTAAATGCTCTGTTACCAgacttac                                                                                        >  2:455839/1‑243 (MQ=255)
                                                                                                                                                      aGATTATCCGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTAcagc                                                                                     >  1:225279/1‑251 (MQ=255)
                                                                                                                                                           atctgTACTATCTGGGCATGCTGGGAATGCATGGCCCCAACGCTGCGAACTTCGCTGTGCAGGAGTGAGACTTGCGGTTCGCCGTGGGGTCACGTTTTGATGGACGGGTTACCGGCACACTGAACACCGTCGCGCCACCCGCCAGTGGTAGCCATATTGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTTTTACCAGCATTACAGCAGcc                                                                                 <  2:384691/246‑1 (MQ=255)
                                                                                                                                                           atccggtctagcTGGGCATGCTGGGAAGGCATGGCACCAAAGCGGAGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCCCGTTTTGATGACCGGGTGAACGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATTAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCg                                                                                <  2:410450/240‑1 (MQ=255)
                                                                                                                                                             ccGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGtt                                                                              >  1:441323/1‑251 (MQ=255)
                                                                                                                                                             ccGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGtt                                                                              >  1:69403/1‑251 (MQ=255)
                                                                                                                                                             ccGTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGtt                                                                              >  1:227896/1‑251 (MQ=255)
                                                                                                                                                               gTACTATCTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGAGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTa                                                                             >  1:38528/1‑250 (MQ=255)
                                                                                                                                                                      cTGGGCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCACCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑Atatca                                                                    >  1:17388/1‑246 (MQ=255)
                                                                                                                                                                       ggggCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑Atatcaat                                                                  <  2:296729/250‑7 (MQ=255)
                                                                                                                                                                        gggCATGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑Atatcaatg                                                                 <  1:182950/251‑8 (MQ=255)
                                                                                                                                                                            aTGCTGGGAATGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTg                                                               >  1:110310/1‑251 (MQ=255)
                                                                                                                                                                                 gggactGCATGGCACCAAAGCGGCGAACTTCGCGGTGCAGAAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCCGTGTTATCCATATGGATATCGACCCGTCAGCAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTNCAGCACCCGTTAA‑‑Atttcactgactggccgc                                                        >  2:245377/6‑235 (MQ=255)
                                                                                                                                                                                       gcatggCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTacaat                                                    <  1:226072/247‑3 (MQ=255)
                                                                                                                                                                                        catggCACCAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACTg                                                   >  1:265353/4‑251 (MQ=255)
                                                                                                                                                                                               ccAAAGCGGCGAACTTCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACTGCGCGCAg                                            <  1:421781/251‑1 (MQ=255)
                                                                                                                                                                                                   aGCGGCGAACTTCGCGGTGCAGGAGCGCGGCTTGCTGATCGCCGGGGGTGCGAGTTTTGATGACCGGGTGACCGGCAAACTTAACTCCTTCGCACCACCCGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACTGCGCGCAGCTGc                                        <  2:157839/251‑1 (MQ=255)
                                                                                                                                                                                                      ggcgaaCTTCGCGGTGCCGGAGTTCGACTCGCTGATCGCCGTGGGTGCTCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACTGCGCGCAGCTGCGtg                                     <  2:166777/247‑1 (MQ=255)
                                                                                                                                                                                                             ttCGCGGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACTGCGCGCAGCTGCGTGATGAAc                               >  1:256567/1‑250 (MQ=255)
                                                                                                                                                                                                              tCGCGGTGCAGGAGTGCGACTTGCTGATAGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAAATGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATAGACCCGGCAGAAATGACCAAGCTGCGCAAGGCACATGGGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATAAATGACTGGCNGCTACAACGCGCTAAGCTGCGGCATGAACAt                             >  2:241764/1‑251 (MQ=255)
                                                                                                                                                                                                              acgcgGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGCTGACCGGCAACCTGAACACCTTCGCACCACACGCCAGTGTTATCAATATGGATATCGACCCGCAAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACACGGGGCTTTAAATGCTCTGTTACCAGCATTACAGCCGCCGTTAAATATCAATGACTTTCAGCTACACTGCGCGCAGCTGCGtgatggccat                             >  2:31058/2‑246 (MQ=255)
                                                                                                                                                                                                                 cgGTGCAGGAGTGCGACTTGCTGATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACTGCGCGCAGCTGCGTGATGAACAtgcc                          >  1:224874/1‑248 (MQ=255)
                                                                                                                                                                                                                             gcgacttgctgATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACTGCGCGCAGCTGCGTGATGAACATGCCTGGCGTTACGAc              <  1:435686/245‑1 (MQ=255)
                                                                                                                                                                                                                             gcgacttgctgATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACTGCGCGCAGCTGCGTGATGAACATGCCTGGCGTTACGAc              <  1:459837/245‑1 (MQ=255)
                                                                                                                                                                                                                               gacttgctgATCGCCGTGGGTGCACGTTTTGNTGACCGGGTTACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTGACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACGGCGCGCAGCTGCGTGATGAACATGCCTGGCGTTTCGACCa            <  2:241100/247‑1 (MQ=255)
                                                                                                                                                                                                                               gacttgctgATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACTGCGCGCAGCTGCGTGCTGAACATGCCTGGCGTTACGAccc            >  2:119942/5‑250 (MQ=255)
                                                                                                                                                                                                                                       gATCGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACTGCGCGCAGCTGCGTGATGAACATGCCTGGCGTTACGACCATCCCGGt     <  1:258864/250‑1 (MQ=255)
                                                                                                                                                                                                                                          cGCCGTGGGTGCACGTTTTGATGACCGGGTGACCGGCAACCTGAACACCTTCGCACCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGGTACCAGCATTACAGCAGCCGTTAAATATCAATGACTGGCAGCTACACTGCGCGCAGCTGCGTGATGAACATGCCTGGCGTTACGACCATCCCGGTGAc  <  1:307395/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
GCGCCAGATGCTGGCAAAAGCGCAAAAACCGATGCTGTACGTTGGCGGTGGCGTGGGTATGGCGCAGGCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAC  >  NC_000913/3951126‑3951607

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