Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,790,153 C→T P58S (CCC→TCC)  lhgD → L‑2‑hydroxyglutarate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,790,1530CT98.3% 197.7 / NA 59P58S (CCC→TCC) lhgDL‑2‑hydroxyglutarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (35/23);  minor base G (1/0);  total (36/23)
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.79e-01

TGCGCCAGCGTGGCTATTTCGCTTACGCCTCTAACCACTACCAGACGCATCAGTAAGCGCAAAGGAATTGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGCCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTT  >  NC_000913/2789907‑2790394
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
aggcccaacgaaaagctctattgtttacgttgggcctaaacgcagggagactccccGCGCAAAGGAATTGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAAATGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCTGCCTGTCACCAGACGGGCCACAACAGCGGCGGGATCCATGCCGGGGGCTATTACACGTCCgg                                                                                                                                                                                                                                               >  2:81109‑M2/57‑251 (MQ=255)
           aaagctctattgtttacgttgggcctaaacgcagggagactccccGCGCAAAGGAATTGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCCTGCCGGGGGCTATTAAACGTCCGGCAGCCTGAAgg                                                                                                                                                                                                                                    >  2:72800‑M2/46‑251 (MQ=255)
            aaccgctattgtttacgtggggcctaaacgcagggagactccccGCGCAAAGGAATTGAGCGGA‑TGTATGATTTTGTTATTATTGGCGGCGGCATCA‑TCGGCATGCCGACCGCCATGCAACTGATTGAGGTCTATCCGGACGCCCGCATTCCGTTG‑CTGGAAGAAGATTCCGCCCCGGCCTGTCCCCAGTCGGGCCACAACATCGGCGTGATCCATGCCGGGGTCTATTACACTTCCGGCAGCCTTAAGgc                                                                                                                                                                                                                                   <  2:365162‑M2/207‑1 (MQ=255)
                      gtttacgttgggcctaaacgcagggagactccccGCGCAAAGGAATTGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGc                                                                                                                                                                                                                         >  1:80413‑M2/35‑251 (MQ=255)
                                ggcctaaacgcagggagactccccGCGCAAAGGAATTGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGa                                                                                                                                                                                                                >  2:400410‑M2/25‑250 (MQ=255)
                                ggcctaaacgcagggagactccccGCGCAAAGGAATTGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGa                                                                                                                                                                                                                >  2:167198‑M2/25‑250 (MQ=255)
                                     aaacgcagggagactccccGCGCAAAGGAATTGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACcgc                                                                                                                                                                                                          >  1:405087‑M2/20‑251 (MQ=255)
                                         ccggggatactccccgcgaaaaGGAATTGAGCGGA‑TCTATGATTTTGTGCTTATGGGGGGCTGCATCT‑CCGGCATGTGGACCGCCATGCAACTGATTGATGTCTACCCGGACTCCCGAATTGCGTTT‑TGGGAAAAAGAGTCCGCCCCGGCCTGTCCCCTGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAATCCGCGcca                                                                                                                                                                                                      <  2:444534/232‑1 (MQ=255)
                                                  ctccccGCGCAAAGGAATTGAGCGGA‑TGTATGCTTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCt                                                                                                                                                                                             >  1:53858‑M2/7‑251 (MQ=255)
                                                   cccccgcgaaaaGGAAGTGAG‑GGAGTGTATGGTTTTGTGATGATTGGCGGGGGCA‑CAGTGGTCATGTCGAGCGCCATGCAACTGGTGGATGTCTATCCGGACGCCCGCATTGCGTTG‑GTGGAAAAAGAGCCCGGCCCGGCGTGACGCCGGAGGCGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGGGGGAAACcgcgtctttaaagcctt                                                                                                                                                                                            <  2:152512/242‑14 (MQ=255)
                                                      ccgggcAAAGGAATTGAGCGGG‑TGTATGATTTGGTGATTATGGGCGGCGGCATCA‑TCGGCATGCCGACCGCCATGCAACTGATTGATGTCGATCCGGCCGCCCGCATTGCGTTG‑CTGGAAAAAGAGCCCGCCCCGGCCTGTCCCCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTg                                                                                                                                                                                         <  2:466718/247‑1 (MQ=255)
                                                        gcgcAAAGGAATTGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATTCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCg                                                                                                                                                                                       >  1:16484/1‑251 (MQ=255)
                                                          gcAAAGGAATTGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACCCGTCCGGCCGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCCAAGCCTTTTGCGAt                                                                                                                                                                                     >  2:219921/1‑251 (MQ=255)
                                                           aaagggAATGGAGGGGA‑TGTATGATTTGGTGATTATTGGGGGCGGCATCA‑TCGGAATGGCGACCGCCAGGCAACTGATTGATGTCTACCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATc                                                                                                                                                                                    <  2:307059/247‑1 (MQ=255)
                                                               ggAATTGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCCAAGCCTTTTGCGATCaaa                                                                                                                                                                                 >  2:145110/1‑250 (MQ=255)
                                                                    tGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑CCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGc                                                                                                                                                                            <  1:219921/250‑1 (MQ=255)
                                                                      aGCGGA‑TGTATGATTTTGTGATTATTGGAGGCGGCATCA‑TCGGTATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACCGGCCACAACAGGGGCGTGATCCAGGCCGGGGTCTATTACACTGCCGGCAGCCTGAAGGCGAAGTTTTGCCTGGCGGGGAACCGGGCCACCAAAGCCATTTGCGATCAAAACGGCAtt                                                                                                                                                                         >  2:22084/1‑251 (MQ=255)
                                                                                aTGATTTTGTGATTATGGGCGGCGGCATCA‑TCGGCATGTCGCCCGCCATGCAACTGATTGATGTCTACCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGAc                                                                                                                                                                <  2:151497/251‑1 (MQ=255)
                                                                                  gATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACa                                                                                                                                                               >  1:166475/1‑250 (MQ=255)
                                                                                     tttGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTg                                                                                                                                                           >  2:459296/1‑251 (MQ=255)
                                                                                        gtgATTATGGGCGGCGGCATCA‑CCGGCATGTCGACCGCCATGCAACTGATTGATGTCTAGCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCg                                                                                                                                                         <  2:260990/250‑1 (MQ=255)
                                                                                          gATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCa                                                                                                                                                      <  1:150290/251‑1 (MQ=255)
                                                                                           attattGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTATCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCa                                                                                                                                                      >  1:466409/1‑250 (MQ=255)
                                                                                              attGGCGGCGGCATCA‑TAGGCATGTCGACCGCCATGCAAATGATTGATGTCTATCCGGAAGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCAACAAAGCCTTTTGCGATCACAACGGCCTTCGCCCCGACAACAGCGGCAAGa                                                                                                                                                   >  2:448647/1‑250 (MQ=255)
                                                                                                    ggcggcATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACCACCGCGGCAAGCTGCTGGt                                                                                                                                            >  2:292553/1‑251 (MQ=255)
                                                                                                            ca‑tcGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGt                                                                                                                                    >  1:436953/1‑251 (MQ=255)
                                                                                                                   cATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGAtc                                                                                                                              >  2:12538/1‑251 (MQ=255)
                                                                                                                   cATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGAtc                                                                                                                              <  1:389632/251‑1 (MQ=255)
                                                                                                                         gACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACTACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGGGATCAAAACGGCACTCGCTACGACAACTGCGGCAAGATGCTCGGTGCCACGACCGTTCTCGaa                                                                                                                         >  2:403513/1‑250 (MQ=255)
                                                                                                                             gCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGGTGCTGGTTGCCACGTCCGATCTCGAAATgg                                                                                                                     >  2:15120/1‑250 (MQ=255)
                                                                                                                                     acttattgatGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATgc                                                                                                            >  1:104741/5‑251 (MQ=255)
                                                                                                                                        ggttgatgccgtgcgggccGCCCGCGTTGCGTTG‑CTGAAAAAAGAGTCCGCCCCGGCCTGTCACCAGCCGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATgcg                                                                                                           <  2:73277/231‑1 (MQ=255)
                                                                                                                                         attgatGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATgcgcgc                                                                                                        >  2:398293/1‑251 (MQ=255)
                                                                                                                                             atGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATACTGGTTGCCACGGCCGATCTCGAAATGGAACGGATGCGCGCCtt                                                                                                     >  1:245827/1‑250 (MQ=255)
                                                                                                                                                   agccGGACGCCCGCAGTG‑GTTGGCGGAAAAAAGCGGCCGCCCCGGCATGTACCCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGGCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATggg                                                                                                <  2:196113/247‑1 (MQ=255)
                                                                                                                                                           gCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAgt                                                                                      >  1:414335/1‑250 (MQ=255)
                                                                                                                                                            cccGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAgtg                                                                                     <  1:12538/251‑3 (MQ=255)
                                                                                                                                                                   ttgcggtgtgggaaaaagttgccGCCCCGGCCGGCCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGaa                                                                                <  2:243653/229‑1 (MQ=255)
                                                                                                                                                                    gCGTTG‑CTGGCAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCCCCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGATGCCACGTCCGATCTCGAAATGGAACGGAGGCGCGCCTTATGGGAACGCACAGTGGCGCACgg                                                                             >  2:132610/1‑251 (MQ=255)
                                                                                                                                                                       ttG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTAt                                                                          >  2:327170/1‑251 (MQ=255)
                                                                                                                                                                        tG‑CTGGAAAAAGCGTCCGCCCCGGCCTGTCCCCAGCCGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTAt                                                                          <  2:462170/250‑1 (MQ=255)
                                                                                                                                                                         g‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCg                                                                        >  1:8621/1‑251 (MQ=255)
                                                                                                                                                                                 aaaGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGa                                                                 <  1:250632/251‑1 (MQ=255)
                                                                                                                                                                                  acgagTCCGCCCCGGCCTGTCACCCGACGGGCCACAACAGCGGCGTGATCCCTGCCGGGGGCTATTACCCGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAAACGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAAATGCTGGTTGCCACGTCCGATCTCGGAATGGAAAGGGTGCGCGCCTTATGGGAACGCCAAGAGGCGGACGGTATCGAGCGCGa                                                                 >  2:279353/3‑250 (MQ=255)
                                                                                                                                                                                     agTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATgg                                                             >  1:140791/1‑251 (MQ=255)
                                                                                                                                                                                     agTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATgg                                                             <  2:70427/251‑1 (MQ=255)
                                                                                                                                                                                        ccGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCtt                                                          >  1:88810/1‑251 (MQ=255)
                                                                                                                                                                                          gCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTa                                                         >  1:77546/1‑250 (MQ=255)
                                                                                                                                                                                             ccGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGGGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGc                                                     >  1:354557/1‑251 (MQ=255)
                                                                                                                                                                                              cGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGcc                                                    <  1:6247/251‑1 (MQ=255)
                                                                                                                                                                                               ggCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGTCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTCACGCcg                                                   >  1:91538/1‑251 (MQ=255)
                                                                                                                                                                                               ggCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCcg                                                   <  2:8621/251‑1 (MQ=255)
                                                                                                                                                                                                          aGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTgcgc                                        >  2:126177/1‑251 (MQ=255)
                                                                                                                                                                                                                  cacaACAGAGGAGTGATGAATGCCGGGGTCTATTACACGTCCCACAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGCTCAAAACGGCATTCGCTACGACAACTGTGGCAAGATGCTGGTTGCCACGTCCGATCTCGAACTGGAACGGATGCGCGCCTTATGGGAACGCCCAGTGGCGAATGGTAGCGAGGGCGAATGGCTTAACGCCGACGGACTtcgcgggcccg                                 >  2:465926/1‑244 (MQ=255)
                                                                                                                                                                                                                     aacaGCGGCGTGATCCATGCCGGGGTCTATTCCCCGTCCGGCAGCCTGAAGGCGCAGTTTTGCCGGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGTCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAcc                             <  2:53858/251‑1 (MQ=255)
                                                                                                                                                                                                                                gATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCg                  >  1:186311/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ccGGGGTCTATTACACGTCCGGCCGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACGGCGGCAAGACGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTcggc          <  2:63844/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gggTCTATTCCACGTCCGGCCGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCGACCAAAGCCTTGTGCGAGCAAAACGGCATTCGCTACGACAACTGCGGCAGGATGCTGGTTGCCACGTGCGATCGCGAAATGGAACGGATGCGCGCCGTATGGGAACGCACAGTGGCGAACGGGATCGAGTGCGAAGGGCTTAACGCCGAGGAACAGCGCGAGCGCGACCCGAAGGTCCCCGCGCTcggcggg       <  2:117342/251‑2 (MQ=255)
                                                                                                                                                                                                                                                 aTTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCAtttt  <  1:459296/250‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
TGCGCCAGCGTGGCTATTTCGCTTACGCCTCTAACCACTACCAGACGCATCAGTAAGCGCAAAGGAATTGAGCGGA‑TGTATGATTTTGTGATTATTGGCGGCGGCATCA‑TCGGCATGTCGACCGCCATGCAACTGATTGATGTCTATCCGGACGCCCGCATTGCGTTG‑CTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGCCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTT  >  NC_000913/2789907‑2790394

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: ‑