Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,420 T→A A229A (GCT→GCA)  thrC → threonine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,4200TA91.9% 215.8 / 3.7 74A229A (GCT→GCA) thrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base A (31/37);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGC  >  NC_000913/4173‑4661
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
gTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                                                                                   >  1:4511/1‑251 (MQ=255)
 tttACCGAATGTGAAAGTGGTTATCCGCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                                                                                   <  2:456319/250‑1 (MQ=255)
       gAATGTGAAAGTGGTAATACTCTATCCACGAGGCAAAATCAGTCCAATGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGCAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGAAGGCGTTTGCTGATGA‑AGACCTGAAAGAGGCGCTAGGGTTAAACACGGCTAAATCGATTAAAAACAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTgcgc                                                                                                                                                                                                                                               >  2:593325/1‑248 (MQ=255)
                                ccACGAGGCAAAATCATGCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCAGCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCt                                                                                                                                                                                                                   <  1:114863/251‑1 (MQ=255)
                                        cAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAACGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtct                                                                                                                                                                                                           >  2:279846/1‑251 (MQ=255)
                                        cAAAATCAGGCCCCTGAAAGAAAAACTGGTCTGTACATTGGGAGGCAATATCGAAACTGTTGCCACCGACGGCGATGTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGTCACTAGGGTTAAAATCGGCTAACTCGATGAACATCAGCCGTTTGCTGGCGCAGATTTGCTCCTACTTGGAAGCTGTTGCGCGGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGtct                                                                                                                                                                                                           <  2:51029/251‑1 (MQ=255)
                                            aTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGt                                                                                                                                                                                                       <  1:537250/251‑1 (MQ=255)
                                               aGTCCACTGCAAGAAAAACTGTTCTGGACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGACAGTGGCGCTGGGGTTAAACTCAGCTCACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCATTTGCGCAGCGGCCGCAGAAAGCGCGCACCCACCTGGTTGtcgcgtgtcc                                                                                                                                                                                                    >  2:581461/1‑243 (MQ=255)
                                                        cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaa                                                                                                                                                                                            >  1:539105/1‑250 (MQ=255)
                                                        cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCAGCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaaa                                                                                                                                                                                           <  1:286116/251‑1 (MQ=255)
                                                         aaGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAAATCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAAc                                                                                                                                                                                          >  1:379782/1‑251 (MQ=255)
                                                          aGAACAACAGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCTATCAACGGCGCTTTCGATGTCTGACAAGCCCTGGTGAAGCAGACGCTTGATGATGC‑GGAGCTGAAAGGGGCGGTGGGGGTAAACTCAGCTAACTCGATTACCATCAGCGGTGTGCTGGCGCAGGTTGGATACTACTTTGAAGCAGTTGCGTCGCTCCCGCAGGAAGCGCACAACCAGCTGGATGTACCGGTGcccacaggaaact                                                                                                                                                                                         >  2:13155/1‑240 (MQ=255)
                                                          aGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACt                                                                                                                                                                                         >  2:536611/1‑251 (MQ=255)
                                                                 cTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGCCGGCGATTTCGATGCCTGTAAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCgg                                                                                                                                                                                  >  2:34376/1‑250 (MQ=255)
                                                                         tACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGTGTTGAcgg                                                                                                                                                                          >  2:372779/1‑251 (MQ=255)
                                                                          aCATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtc                                                                                                                                                                                                            >  1:137070/1‑216 (MQ=255)
                                                                          aCATTGGGCGGCAATATCGAAACTGTGGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtc                                                                                                                                                                                                            <  2:137070/216‑1 (MQ=255)
                                                                           cATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGAcggcg                                                                                                                                                                        >  1:397026/1‑251 (MQ=255)
                                                                              tGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGt                                                                                                                                                                                                       >  1:487639/1‑217 (MQ=255)
                                                                               gggCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGt                                                                                                                                                                     <  1:313348/250‑1 (MQ=255)
                                                                               gggCGGCAATAACGAACCTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGTAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTc                                                                                                                                                                    <  2:122206/251‑1 (MQ=255)
                                                                                   ggcAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGACTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCACACGGAAACTTCGGCGATTTGCCGGCgggggtgtt                                                                                                                                                                >  2:319302/1‑245 (MQ=255)
                                                                                     cAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCCGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGCAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGt                                                                                                                                                                                                                                                    >  1:645176/1‑165 (MQ=255)
                                                                                     cAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATAAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGtacgctgg                                                                                                                                                              >  1:380436/1‑244 (MQ=255)
                                                                                        tatCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTggggg                                                                                                                                                           >  1:393309/1‑248 (MQ=255)
                                                                                                  gTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCg                                                                                                                                                 <  1:104744/250‑1 (MQ=255)
                                                                                                     gCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCt                                                                                                                                             <  1:536611/251‑1 (MQ=255)
                                                                                                      ccatggccggagagttcggggCCTGTCAGGCGCTGTTGAAGCGTGCGTTTGAGGA‑GAGAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACGCGAGTAACTTCAGCCGTTGGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCTCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGGCACTCGGTCTg                                                                                                                                            <  2:588811/231‑1 (MQ=255)
                                                                                                      ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGGAGGAGTTTGATGATGA‑AGAGCTGAAAGTGGCGATGGGGTTAAACTCAGCTAACTCGATTAAAATCAGCCGTTTGCTGGCGCAGAAGTGCTACTACTTTAAAGAAGTTGCGAAGATGCAGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAACATTCGGCGATTTGACGGCGGGActgtgggcgggagggacaggggtg                                                                                                                                             >  2:659677/1‑230 (MQ=255)
                                                                                                      ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCCGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCCGGtgggctgggggaggaccgtggatgg                                                                                                                                             >  2:514265/1‑227 (MQ=255)
                                                                                                      ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCg                                                                                                                                                                                               >  1:452215/1‑201 (MQ=255)
                                                                                                      ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCg                                                                                                                                                                                               <  2:452215/201‑1 (MQ=255)
                                                                                                        aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGcc                                                                                                                                          <  1:245566/251‑1 (MQ=255)
                                                                                                        aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCGGCTGGCggcgacaacccgcgcggc                                                                                                                                           >  2:312113/1‑234 (MQ=255)
                                                                                                           gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                                                                                   >  1:146957/1‑144 (MQ=255)
                                                                                                           gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                                                                                   <  2:146957/144‑1 (MQ=255)
                                                                                                            aCGGCGATTTCGATTCCTGTCAGGCGCTGGTGAAGCAGGCGGTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTg                                                                                                                                      <  1:330595/251‑1 (MQ=255)
                                                                                                             gggCGATTTGAATGCCTGTCAGGCGTGGGTTAAGCAGGCGTTGGAGGAGGA‑AGAGCGTAAAGTGGCGCTGGGTTTAAACTCAGCTAACGCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGa                                                                                                                                     <  2:314301/250‑1 (MQ=255)
                                                                                                             cGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTg                                                                                                                                                                                                      >  1:401063/1‑187 (MQ=255)
                                                                                                             cGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTGGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTg                                                                                                                                                                                                      <  2:401063/187‑1 (MQ=255)
                                                                                                              ggCGATTTCGATGCCTGTCAGGCGCTGGCGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGTGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAAACAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGaa                                                                                                                                    <  1:259989/251‑1 (MQ=255)
                                                                                                                cGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGATGCTGGCggaggtccaccgtctgccggggaga                                                                                                                                    >  1:232509/1‑226 (MQ=255)
                                                                                                                  aTTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGt                                                                                                                                <  1:617040/251‑1 (MQ=255)
                                                                                                                     tCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGttt                                                                                                                              <  1:431307/250‑1 (MQ=255)
                                                                                                                           ccTGGCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACGCAGCTAAGCCGATGAACATCAGCCGGTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCTATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATtgct                                                                                                                       <  2:397026/251‑1 (MQ=255)
                                                                                                                                cAGGCGCTGGTGAAGCAGGCGGTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGAc                                                                                                                  <  2:316306/251‑1 (MQ=255)
                                                                                                                                          ggagaccggCGTTTGATGATGA‑AGAGCTGAAAGGGGCGCTGGGGTAAAACTAAGCTAACTCGATTAACAACAGCGGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAGCGTGAAc                                                                                                        <  2:4511/244‑1 (MQ=255)
                                                                                                                                                aGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGTTTAAACTCAGCTAACTCGCTTAACATCGGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCACCCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTa                                                                                                                            <  2:542452/225‑1 (MQ=255)
                                                                                                                                                 ggCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCg                                                                                                 <  1:268960/251‑1 (MQ=255)
                                                                                                                                                       tgatgatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTTTTTCTGCGACCAACGTGAACGATACCGTGcccc                                                                                           >  1:193019/1‑249 (MQ=255)
                                                                                                                                                          ggatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTct                                                                                                                                                                   <  2:19485/176‑1 (MQ=255)
                                                                                                                                                           gatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACCCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCGGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑ttt                                                                                      <  2:193019/251‑1 (MQ=255)
                                                                                                                                                                 agaaCTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGc                                                                                 <  1:312113/251‑1 (MQ=255)
                                                                                                                                                                  gagcTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTgg                                                                                                                                                                                                                 <  1:156575/121‑1 (MQ=255)
                                                                                                                                                                  gagcTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTgg                                                                                                                                                                                                                 >  2:156575/4‑124 (MQ=255)
                                                                                                                                                                  gagcTGAAAGTGGCGCGGGGGTTAACCTCAGCGAACTCGATGAACAGCAGCCGTTTGCTGGCGCAGATTTGGTACCACTTTGAAGCAGTTGCGCAGCTGCCGTAGGAAGCGCGCAAACAGATGGTTGTCTCGGTGCCAAGCGGAAACGTCGGCGATTTTACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCa                                                                                <  2:384063/248‑1 (MQ=255)
                                                                                                                                                                        aaaGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACgg                                                                         <  1:36973/251‑1 (MQ=255)
                                                                                                                                                                            tGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTAACTCGGTCTGCCGGTGGAAAGTTTTTTTGTTGCGACCCACGTGAAAGATAACGTGCCCCG‑TTTTCTGcccgccgggccg                                                                      >  1:428663/1‑240 (MQ=255)
                                                                                                                                                                             ggCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCTGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGGCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCATGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAg                                                                     <  1:222900/250‑1 (MQ=255)
                                                                                                                                                                               cgcTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTgg                                                                  <  1:610749/251‑1 (MQ=255)
                                                                                                                                                                               cgcTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGA‑AGGCACTCGGGCTGCCGGTGGAACGGTTTATTGCTGGGACCAAAGTGAAAGACACCGGGGCACG‑CTTCCTGCCCG‑ACGGTCAGTgg                                                                  >  1:104366/1‑251 (MQ=255)
                                                                                                                                                                                  tggggTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCCGGCGA‑AGTCACCCGGGCCGCCGGTGGAACGTTTTTTTGCTTCGACCAACGGGAAAGAAACCGTGCCAAG‑TTTTCTGCACG‑ACCGTCCGTgggcc                                                               >  1:587991/3‑248 (MQ=255)
                                                                                                                                                                                           aCTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCg                                                      <  1:500122/251‑1 (MQ=255)
                                                                                                                                                                                              catcTAACTCGAGTAACCTCAGCCGTTTGCTGGCGCAGATTGGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCGACt                                                   <  2:137228/248‑1 (MQ=255)
                                                                                                                                                                                                gCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTTGCGA‑AGACACTCCGTCTGCCGGGGAAACGTTTTTTTTCTTCGGCCAACCGGAACGATACCGTGCCACG‑TTTCCTTCACG‑ACGGTCAGTGGTGACCCCAACCGACTcc                                                 >  1:208808/1‑250 (MQ=255)
                                                                                                                                                                                                                 cAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGG‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAAcg                                <  1:106865/251‑1 (MQ=255)
                                                                                                                                                                                                                  aGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGG‑AGGCCCTCGGGCTTCCGGTGGAACGGTTTTTTTCTTCGGCCAACGGGGAAGATAACGTGGCCCG‑TTTTCTTCCCGCACG‑TCCGGGGTCGCCCAAAACGACTCAGGCGAAGGTAACCAAcgc                               >  1:234986/1‑251 (MQ=255)
                                                                                                                                                                                                                       tttGCTTGCGCAGATTTGCTACTACTTTTAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGAt                            <  1:659677/249‑1 (MQ=255)
                                                                                                                                                                                                                       tttGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTCCTGGCGACAG‑CACTCGGTCTCCCGGGGAAACGGGTTATTGCTGGGACCAAAGTGAAAGGTACCGTGCCACGGTTTCC‑GTACGCAGGGCCCG‑GGTCACCCAAAACGACCCAGGCGACGTTTACCAACGCGATgg                          >  1:434441/1‑251 (MQ=255)
                                                                                                                                                                                                                          gCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACg                       <  1:479856/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tttGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGAc                        <  1:265871/239‑1 (MQ=255)
                                                                                                                                                                                                                                         ctactaCTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGGCTGCCGGTGAAACGTTTTTTTTCTGCGCCCCACGTGAACGAAAACGTGCCCCG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCCAAACGACTCAGGCGGCGTTATCCAACGCGATGGAAGTGTGTCAACCGaga         >  1:658559/1‑248 (MQ=255)
                                                                                                                                                                                                                                           actactTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGG‑AGACACCCGGGCTGCCGGGGGAACGTTTTTTTGCTTCGACCAACGTGAACGAAACCGGGCCACG‑GTTTCTGCCCG‑ACGGGCCGGGGTCACCCCAAACGGCTCAGGCGACGTTGTCCCACGCGATGGACGTGGGTCAGCCCaacaac      >  1:532088/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ttttaaGCAGTTGCGCAGTGGCCGCAGGAGTCGCGCAACCAGCGGGTTGTCGCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCGCCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGc  <  2:537433/246‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
GTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGC  >  NC_000913/4173‑4661

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