Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,599,947 G→T L14F (TTG→TTT gcvR → putative transcriptional regulator GcvR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,599,9470GT88.7% 125.6 / 12.8 53L14F (TTG→TTTgcvRputative transcriptional regulator GcvR
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base T (30/17);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-02

CATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTGGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAACTGGATCTTTTAATCGTGATGAAGCGCACGACGGCGCGTCCGCGTCCGCCAATGCCAGCAT  >  NC_000913/2599706‑2600188
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
caTTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGAt                                                                                                                                                                                                                                             >  1:233155/1‑251 (MQ=255)
               ccatcTGAGTAACAATCGCGACAATACTTCCCGTGAACTTGTGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAATCAGACAGAACCTTTCTTATTGTTTTATGCATGTTTTTGTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACTCTGTCATCGCAACATTATCTGGTGATCACTTCGTTTGGTGCCGATCGCCCTGGAATtgt                                                                                                                                                                                                                               <  2:542104/250‑1 (MQ=255)
               ccatcGGAGTAATAACCGCGACAATACGTCCCGTGAACATTGGCCATCCTCTGTGCAAACAAGTGTCACAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTGTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATAACTGCGTTTGGTGCCGTTCGCCCTGGAATtgt                                                                                                                                                                                                                               <  2:113507/250‑1 (MQ=255)
                 atcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGa                                                                                                                                                                                                                             >  2:90275/1‑250 (MQ=255)
                            aaTCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCAccc                                                                                                                                                                                                                 <  1:462074/251‑1 (MQ=255)
                                    cAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTGGGTGCCGATCGCCCTGGAAt                                                                                                                                                                                                                                  >  1:130573/1‑226 (MQ=255)
                                    cAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTGGGTGCCGATCGCCCTGGAAt                                                                                                                                                                                                                                  <  2:130573/226‑1 (MQ=255)
                                     aaTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCACTGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACAGATCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGCACACCATCCCCCGTCATGTcc                                                                                                                                                                                                        >  2:191976/1‑250 (MQ=255)
                                                           cagccTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGTTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGTAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGTCTGTAATATTGa                                                                                                                                                                                  <  2:200117/248‑1 (MQ=255)
                                                            aTCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTGGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGaa                                                                                                                                                                                 <  1:27911/251‑1 (MQ=255)
                                                              ccTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAg                                                                                                                                                                                <  1:341013/250‑1 (MQ=255)
                                                                       aaaCAAGTGTCTCAATGGTACGGTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTATTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCt                                                                                                                                                                    <  2:5653/251‑1 (MQ=255)
                                                                                    aaTGGTTCGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGga                                                                                                                                                       <  1:372354/251‑1 (MQ=255)
                                                                                    aaTGGTACGTTTGGTATGGCATTAAAAGCAATCAGACAGAACCGTTCTGCTTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTTTCATCGCAACATTATCTGGTGATCACTGCGATTGGTGCCGATCGCCCTGGAATTGTGAACACCCTCACCCGTCATGTCAGTAGTTGCGGCTGTAATTTTGAAGACAGTCGCCTGGCGCTGCTgggg                                                                                                                                                       >  2:243129/1‑250 (MQ=255)
                                                                                       ggTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGgaga                                                                                                                                                     >  1:263915/1‑250 (MQ=255)
                                                                                            gTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGtt                                                                                                                                              >  1:331306/1‑251 (MQ=255)
                                                                                                                caagcaGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCTCCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCttt                                                                                                                          <  1:422911/251‑1 (MQ=255)
                                                                                                                caagcaGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCttt                                                                                                                          <  1:90275/251‑1 (MQ=255)
                                                                                                                     agacagaACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATAACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGACTGGCGAAGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCggg                                                                                                                     >  2:391736/1‑250 (MQ=255)
                                                                                                                               gTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGTAGAAGAGTTCCCGTTTATTTTGCTGCTTTCCGGTTCATGGAATg                                                                                                           >  2:334928/1‑251 (MQ=255)
                                                                                                                                ttCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGc                                                                                                          >  2:334963/1‑251 (MQ=255)
                                                                                                                                         ttgtATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTActct                                                                                                 >  1:214631/1‑251 (MQ=255)
                                                                                                                                         ttgtATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTActct                                                                                                 >  1:13097/1‑251 (MQ=255)
                                                                                                                                          tgtATGAATGTATTTTTTATGCTTTACTTCAGAACAACTCAGCCCTTAAAGGAATAACAAGTTTGACACTGTTATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCAATCGTCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTACTTGCGGCTGTAATATTGAAACAAGTCCGCCGGCAGTGCTGGGGATAGTGTTCACGTTTATTATGCTGCATTCCGGTTctgggatggcacttcctttg                                                                                                >  2:18378/1‑232 (MQ=255)
                                                                                                                                            tATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGAt                                                                                              >  1:186449/1‑251 (MQ=255)
                                                                                                                                                         tttATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGtt                                                                                  <  1:191976/250‑1 (MQ=255)
                                                                                                                                                               cTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGACTCAACGTTACCGttt                                                                           >  2:92132/1‑250 (MQ=255)
                                                                                                                                                                  tCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGGGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGaa                                                                         >  1:318562/1‑250 (MQ=255)
                                                                                                                                                                        aGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGcc                                                                  <  1:47437/251‑1 (MQ=255)
                                                                                                                                                                          aacaacTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCACAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGATTCATGGAATGCCATTACTCTGATTGAATCAAAGTTACCGTTGAAAGGAGCCg                                                                 >  1:380504/1‑250 (MQ=255)
                                                                                                                                                                           acaacTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGa                                                                >  1:28954/1‑250 (MQ=255)
                                                                                                                                                                            caacGCACCCCTTAAAGGAATAACCAGTTTGCCACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGGGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAAc                                                              <  2:4880/251‑1 (MQ=255)
                                                                                                                                                                                  aCCCCTTAAAGCAATAACCAGGTTGACCCTGTCATCGCAACATTATATGATGATCAATGCGTTTGGTGCCAATCGCCCTGGCATTGTGAACACCATAACCCGTCATGTCAGTAGATGCGGCTGTAATGTTGGAGACAGTCGCCTGGCGATGCTGGGAGAAGAGGTCACGGTTATTATGCTGCTTTCCGGTTAATGGAATGCGCTTACTCTGGTTGACTCAACGTTACCGTTGAAACATGCCGTACTGGAtt                                                        >  2:79209/1‑250 (MQ=255)
                                                                                                                                                                                      cTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGACGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCCTTACTCTGATTGAATCAACGTTACAGTTGAAAGGTGCCGATCTGGATCtttt                                                    >  2:369699/1‑251 (MQ=255)
                                                                                                                                                                                      cTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCCTTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAACTGGATCtttt                                                    >  2:383488/1‑251 (MQ=255)
                                                                                                                                                                                       ttAAAGGAATAACCAGTTTGACACTGTAATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACACGTCATGACAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGGGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGGGCCGACCGGGATCtttt                                                    >  2:100858/1‑250 (MQ=255)
                                                                                                                                                                                         aaaGGAATAACCAGTTTGGCACTGTCATCGTAACATTAGCTGGGTATCACTGCGTTTGGTGCCGATCTCCCTGGAATTGTTAACACCATCACGCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCCGGCGATGCTGGGAGAAGATTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAAGGCCATTACTCTGATTGAATCAACGTTACCGTTTAAAGGTGCCGAACCGGATCTCTTAAt                                                 <  2:85572/251‑1 (MQ=255)
                                                                                                                                                                                             gAATAACCAGTTTGACTCTGTCATCGCAACATGATCTGGTGATCACTGCGGTTGGTGCTGATCGACCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCAc                                                                                                                                           <  2:337250/157‑1 (MQ=255)
                                                                                                                                                                                             gAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCAc                                                                                                                                           >  1:337250/1‑157 (MQ=255)
                                                                                                                                                                                                  aCCAGTTTGACACTGTCATCGCAACATTAGCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAACTGGATCTTTTAATCGTGATGaa                                        <  2:263915/251‑1 (MQ=255)
                                                                                                                                                                                                   ccAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAACTGGATCTTTTAATCGtgatga                                         >  1:429910/1‑249 (MQ=255)
                                                                                                                                                                                                     aGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAACTGGATCTTTTAATCGTGATGAAgcg                                     >  1:108729/1‑251 (MQ=255)
                                                                                                                                                                                                         tGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAACTGGATCTTTTAATCGTGATGAAGCGCacg                                 >  1:163971/1‑251 (MQ=255)
                                                                                                                                                                                                            cacTGTCATCGCAACATTATCTTGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATAACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCCCGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGAGCCGAAATGGATCTTTTAATCGTGATGAAgcgccagccg                              >  1:513529/1‑245 (MQ=255)
                                                                                                                                                                                                                        aaCATTATCTGGTGATCACTGCGTTGGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTc                                                                                                                                             >  1:516334/1‑128 (MQ=255)
                                                                                                                                                                                                                        aaCATTATCTGGTGATCACTGCGGTGGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTc                                                                                                                                             <  2:516334/128‑1 (MQ=255)
                                                                                                                                                                                                                          cATTATCTGGTGATCACTGCGGTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGGTCATGGAATGCCATTACTCTGATTGAAGCAACGTTACCGTTGAAAGGGGCCGAACTGGATCTTTTAATCGTGATGAAGCGCACGACGGCgcgtccgcgtcc                <  1:497566/251‑1 (MQ=255)
                                                                                                                                                                                                                           aTTATCTGGTGATCACTGCGTTGGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTACGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTTCCGTTGAAAGGTGCCGAACTGGATCTTTTAATCGTGATGAAGCGCACGACGGCGcgtccgcgtccg               >  2:133423/1‑251 (MQ=255)
                                                                                                                                                                                                                                cTGGTGATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAACTGGATCTTTTAATCGTGATGAAGCGCACGACGGCGCGTCCGCGTCCGCCAAt          >  2:190416/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tGGTGATCACTGCGTTTAGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGGGAAGAGTTCCCGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGAGCCGACATGGGTCTTTTATTCTTGATGAAGCGCCAGACGGcgcgcccactaccgccccg          >  2:204304/1‑235 (MQ=255)
                                                                                                                                                                                                                                     gATCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAACTGGATCTTTTAATCGTGATGAAGCGCACGACGGCGCGTCCGCGTCCGCCAATGCCa      <  1:86301/250‑1 (MQ=255)
                                                                                                                                                                                                                                      aTCACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAACTGGATCTTTTAATCGTGATGAAGCGCACGACGGCGCGTCCGCGTCCGCCAATGCCAGc    >  1:488448/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cACTGCGTTTGGTGCCGATCGCCCTGGAATTGTGAACACCATCACACGTCATGTCAGTAGTTGCGGCTGGAATATTGAAGACAGACGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGCATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAACTGGATCTTTTAATCGTGCTGAAGCGCACGACGGCGCGTCCGCGTCCGCCAATGCCAGCAt  >  1:521900/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
CATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACATTATCTGGTGATCACTGCGTTGGGTGCCGATCGCCCTGGAATTGTGAACACCATCACCCGTCATGTCAGTAGTTGCGGCTGTAATATTGAAGACAGTCGCCTGGCGATGCTGGGAGAAGAGTTCACGTTTATTATGCTGCTTTCCGGTTCATGGAATGCCATTACTCTGATTGAATCAACGTTACCGTTGAAAGGTGCCGAACTGGATCTTTTAATCGTGATGAAGCGCACGACGGCGCGTCCGCGTCCGCCAATGCCAGCAT  >  NC_000913/2599706‑2600188

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