Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,441 2 bp→AG coding (708‑709/1287 nt) thrC → threonine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,4410GA96.6% 175.6 / ‑4.3 59E236E (GAG→GAA) thrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (28/29);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
*NC_0009134,4420AG96.6% 177.5 / ‑4.0 59T237A (ACG→GCG) thrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (28/29);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGA‑GACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACT‑CAGGCGACGTTATCCAACGCG‑ATGGACG‑TGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGC  >  NC_000913/4196‑4684
                                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                       
aTCCTCTATCCACGAGGCTAAATAAGTTCACTTCAAGAAAAACTTTACTGTCCATTGGTAGTCAATATCGAAACTGTTGCCATCGACGGCGATTTCTATGCCTGTCAGTCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAATCTCGGCTAACTCGATTAACATCAGCGGGATGCTGGCTCAGATTTCCTACTACTTTGAAGCTGTTTCGCAGCTGCCGCAGGT‑GActcg                                                                                                                                                                                                                                                     <  2:409211/251‑4 (MQ=255)
     ctATCCACGAGGCAAAATCATTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCTATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTGGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCaa                                                                                                                                                                                                                                                  <  2:585389/249‑1 (MQ=255)
     ctATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAAcc                                                                                                                                                                                                                                                >  2:213540/1‑251 (MQ=255)
                gCAAAATCAGTCCACTGCAAGAAAAACTTTTCTGTACATTGGGCTGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGGAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAATCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGtc                                                                                                                                                                                                                                     <  2:357277/251‑1 (MQ=255)
                gCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACGTCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGtc                                                                                                                                                                                                                                     >  1:32305/1‑251 (MQ=255)
                   aaagcAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCg                                                                                                                                                                                                                                  <  1:604965/247‑1 (MQ=255)
                    aaTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGACACTGATGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCgg                                                                                                                                                                                                                                 >  2:149560/1‑251 (MQ=255)
                     aTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCCGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCggg                                                                                                                                                                                                                                >  2:431329/1‑250 (MQ=255)
                            cACTGCAAGAAAAACTGTTCTGTTCAGTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGAGGCCTGTCAGGCGCTGTTGAAGCAGGCGTTTGATGATGTAGAGCTGAAAGTGGCGCTGGTGTTAAACGCAGCTAACTCGATGAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAATCAGTTGCCCAGCTGTCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGc                                                                                                                                                                                                                         <  2:277319/251‑1 (MQ=255)
                                gCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCATCCAGCTGGTTGTCTCGGTGCCAAGCGGaa                                                                                                                                                                                                                     <  2:526402/251‑1 (MQ=255)
                                 cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaa                                                                                                                                                                                                                     >  1:419795/1‑250 (MQ=255)
                                    gAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACtt                                                                                                                                                                                                                 >  2:330529/1‑251 (MQ=255)
                                       aaaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAATCGGAATCTTCgg                                                                                                                                                                                                              <  1:499889/251‑1 (MQ=255)
                                       aaaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCgg                                                                                                                                                                                                              >  2:105044/1‑251 (MQ=255)
                                             ttCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGAttt                                                                                                                                                                                                        >  1:315340/1‑251 (MQ=255)
                                               cTGTCCATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTg                                                                                                                                                                                                       >  2:129887/1‑250 (MQ=255)
                                                       tGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCgggg                                                                                                                                                                                              >  2:511917/1‑250 (MQ=255)
                                                               aaTATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTggg                                                                                                                                                                                      >  1:15404/1‑250 (MQ=255)
                                                                     gAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCggaggc                                                                                                                                                                                >  2:174817/1‑246 (MQ=255)
                                                                           gtggcctacGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTGGATGATGAAGAACTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGAGTTGCTGCTACTGTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACGTCGGCGAGTTGACGGCGGGTCTGCTGGCGATGTCACTc                                                                                                                                                                            <  2:547897/241‑1 (MQ=255)
                                                                           gTTGCCATCGACGGCTATTTAGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACGCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACGACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGTCACTCg                                                                                                                                                                           <  2:451502/250‑1 (MQ=255)
                                                                              gtcaacGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGGTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCACCCAGCTGGTTGTCTCGGTGCCAAGCGGAACCTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCt                                                                                                                                                                       <  2:290957/246‑1 (MQ=255)
                                                                              gCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCgggggtgccggcggagtcccccgggct                                                                                                                                                                       >  2:231244/1‑227 (MQ=255)
                                                                                    gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGt                                                                                                                                                                 <  1:171426/251‑1 (MQ=255)
                                                                                          gATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGGCACCCGGTCTGCCGGGGAAACg                                                                                                                                                           >  1:172013/1‑251 (MQ=255)
                                                                                              tCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGCGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGCAGTCACTCCGTCTGCCGGTGaaaaggttt                                                                                                                                                       >  2:14570/1‑245 (MQ=255)
                                                                                               cGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGtttt                                                                                                                                                       <  1:529948/250‑1 (MQ=255)
                                                                                                                 ggTGAAGCCGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTTTTTCTGCGGCCCACGTgg                                                                                                                                    >  1:200014/1‑250 (MQ=255)
                                                                                                                    gAAGCAGGTGGTTGAAGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAATGCAGCTAACTCGATTAACCTCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGatcg                                                                                                                                 <  2:233944/251‑4 (MQ=255)
                                                                                                                                  atgatgaAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCCGGCGAAAGCACACGGTCGGCCGGTGGAACGTTTTATTGCTGCGACCCAAGGGGACGCAACCGTG‑CCACGt                                                                                                                  >  1:404046/1‑251 (MQ=255)
                                                                                                                                      tgaAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGCAGTCACCCGGTCTGCCGGGGTAACGTGTTTTTTCTGCGCCCAACGTGAACGAGAACGTG‑CCACggtttc                                                                                                              >  2:506223/1‑246 (MQ=255)
                                                                                                                                        aagagCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAATTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTg                                                                                                            <  1:563956/251‑1 (MQ=255)
                                                                                                                                           aaCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACGCGATTAACATCAGCCATTTGATGGCGCAGATTTGCTACTATTTTTAAGCAGTTGCGCAGCTGCCGCAGGAAGA‑GCGCAACCAGCTTGTTGTCTCGGTGCCAAGCGGAAACTTCGGCAATTGGACGGCGGGTCTGCTGTCGACGGCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATGCCTTG‑CCGCGTTTCCTGcgc                                                                                                         <  2:7764/251‑3 (MQ=255)
                                                                                                                                            gcTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCacg                                                                                                        <  1:14570/250‑1 (MQ=255)
                                                                                                                                               gAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGGTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCacgac                                                                                                      <  1:104435/250‑1 (MQ=255)
                                                                                                                                                     ggCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAg                                                                                                <  1:129887/250‑1 (MQ=255)
                                                                                                                                                             ggTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGGCACTCGGTCTGCCGGGGGAACGGTTTTTTGCTGCGACCAACGCGAAAGATACCGTG‑CCACGCTTTCTGCCAGACGGTCAGGGGGCA‑ccc                                                                                      >  2:430617/1‑251 (MQ=255)
                                                                                                                                                                      cagCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACt                                                                             <  1:105044/249‑1 (MQ=255)
                                                                                                                                                                        gCTAACTCGATTAACATCAGCCGTTTGCGGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTTTCGGTGCCAAGCTGAAACTTCGGCTATTTGACGGCGGGTCTTCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATTCCGTG‑GCACGTTTCCTGCACGACGGTCAGTGTTCT‑CCCAAAGCGACT‑Ca                                                                          <  1:538589/251‑1 (MQ=255)
                                                                                                                                                                            aCTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACT‑CAGGCg                                                                      <  1:337115/251‑1 (MQ=255)
                                                                                                                                                                               cGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTTGGGGAGTCACTCGGTCTGCCGGGGGAACGTTTTTTTTCTGCGACCCACGTGAACGATACCGTG‑CCACGTTTTCTGCACGACGGGCCGTGGTCA‑CCCCAAGCGACC‑CAGGCGACg                                                                   >  1:2605/1‑251 (MQ=255)
                                                                                                                                                                                    aaCATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTTCTGGCGAAGTACCTCGGTCTGCCGGGGCAACGTGTTATTGGTGCGACCAACGGGAACGAAACCGTG‑CCCCCTTTCCTGCACGGCGGGCCGTGGTCAACCCAAA‑CGACT‑CAGGccgcggtatc                                                              >  1:364483/1‑242 (MQ=255)
                                                                                                                                                                                    aaCATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACT‑CAGGCGACGTTAt                                                               <  1:290014/250‑1 (MQ=255)
                                                                                                                                                                                     aCATCAGCCGTTTGCTGGCGCAGATTTGCTACTACGTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCGGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACT‑CAGGCGACGTTATcc                                                             <  2:295962/251‑1 (MQ=255)
                                                                                                                                                                                             cGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGGCACTCGGGCTGCCGGTGGAACGTTTTATTGCTGCGACCAACGGGAACGATACCCTG‑CCACGCGTCCTGTCCGACGGGCCGTGGTCA‑CCCCAAACGACT‑CAGGCCACGGTATCCAAagcg                                                       >  2:313235/1‑249 (MQ=255)
                                                                                                                                                                                                     ggCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGTTCA‑CCCAAAGCGACT‑CAGGCGACGTTATCCAACGCG‑ATGGACG‑t                                             <  1:520112/249‑1 (MQ=255)
                                                                                                                                                                                                                gCTACTAATTTGAAGTAGTTGCGCAGCTGCCGCAGTA‑AGCGCGCAACCAGCGGGTTGTCGCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAATCGACT‑CATGCGCCGTTATCCAGCGCG‑ATGTACG‑TTAGTCAGCCGaac                                <  2:623721/251‑1 (MQ=255)
                                                                                                                                                                                                                         ttGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACT‑CAGGCGACGTTATCCAACGCG‑ATGGACG‑TGAGTCAGCCGAACAACTGGcc                        <  2:627517/250‑1 (MQ=255)
                                                                                                                                                                                                                          tGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGACTGCTGGCGAAGGCACTCGGTCTGCCGGTGAAACGGTTTATTGCTGCGACCCACGTGAACGATACCGTG‑CCACGGTTCCTGCCCGCCGGTCAGTGGTCC‑CCCAAAACGACT‑CAGGAGGCGTTATCCCACGCG‑ATGGACG‑CGAGGCAGACGAACAACTGGCcgc                      >  1:257674/1‑251 (MQ=255)
                                                                                                                                                                                                                           gAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGGCACTCGGTCTGCCGGGGAAACGTTTTATTGCTGCGACCAACCGGAACGATAACGTG‑CCACGTTTCCTGCACGACGGGCCGTGGTCA‑CCCCAAGCGACT‑CAAGCGACGGTATCCAACGCG‑GTGGACG‑GGGGTCCACCGAACAACTTGCcgcc                     >  1:502626/1‑250 (MQ=255)
                                                                                                                                                                                                                             agcagTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACT‑CAGGCGACGTTATCCAACGCG‑ATGGACG‑TGAGTCAGCCGAACAACTGGCCGCgtg                   <  1:330529/247‑1 (MQ=255)
                                                                                                                                                                                                                              gcagTTGCGAAGCTGCCGCAGGA‑AGCGAGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGTCACTCGGTCTGCCGGGGAAACGGTTTTTTTCTTCGACCAACGAGTACGGTAACGTGGCCA‑GTTTCCTGCAGGGCGGGCAGGGGTTA‑CCCAACGCAACTTCAGG‑GCCGTTATCCAAAGCGGATGGA‑G‑GGTGGAAGCCGAACAACTAGCcgcgggg                  >  2:357705/4‑248 (MQ=255)
                                                                                                                                                                                                                               cagTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTTTGCTGGAGGAGTCACTCGGTCTGCCGGGGGAACGTTTTATTTTTGCGACCCACGGGGAAGATACCGAG‑GCACGGTTTCTGCACGGCGGGCAGTGGGCA‑ACCAAAACGACC‑CAAGCGACGGTTTCCAAACCG‑ATGGACG‑GGCGTCAGCCGAACAACTGGCcgcggggg                 >  2:588288/3‑247 (MQ=255)
                                                                                                                                                                                                                                    gcgcAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACT‑CAGGCGACGTTATCCAACGCG‑ATGGACG‑TGAGTCAGCCGAACAACTGGCCGCGTGTGGAaga            <  1:370704/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cAGCGGCCGCAGGA‑AGCGCGCAACCAGCTGGTTTTCTCGGTGCCAAGCGGAAACTTCGGGGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACT‑CAGGCGACGTTATCCAACGTG‑ATGGACG‑TGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAgtt         <  1:369766/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGCTTTGACGGCGGGTCTGCTGGCGAAGTAACTCGGTCTTCCGGGGAAAAGTTTTATTTCTGCGAACAACGTGAACCATACCGTG‑CCACGGTTTCTGCACGACGGGCCGTGGTCA‑CCCCAAGAGACT‑ACGGAGACGTTATTCAACGCC‑ATGGACGAGGGGCCAGCC‑AACAACCGGCCGCGTGTGGgagaagtgttt     >  2:155481/1‑244 (MQ=255)
                                                                                                                                                                                                                                            gCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACT‑CAGGCGACGTTATCCAACGCG‑ATGGACG‑TGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGGTGTTcc    <  1:149560/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ccGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACT‑CAGGCGACGTTATCCAACGCG‑ATGGACG‑TGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTccg   <  2:172013/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGGAACGTTTTTTTGCTGCGACCCACCTGCAAGATACCGTG‑ACCCGTTTTCTGCACGACGGGCAGGGGGCA‑CCCAAAGCGACT‑CAGGGGACGGTATCCCACCCG‑AAGGGAG‑GGAGTCAGCCCGACAACTGGCCCCGTGTGGGAGAGTTGTTccgc  >  1:576247/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                       
ATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGA‑GACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTG‑CCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACT‑CAGGCGACGTTATCCAACGCG‑ATGGACG‑TGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGC  >  NC_000913/4196‑4684

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