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,4410GA93.1% 195.9 / 1.9 72E236E (GAG→GAA) thrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base A (30/37);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,4420AG95.8% 211.0 / ‑1.8 72T237A (ACG→GCG) thrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (30/38);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

CCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCC  >  NC_000913/4205‑4685
                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                                      
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)
                                                     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                                                                                                                                                                                                                       <  2:452215/201‑1 (MQ=255)
                                                                      ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCg                                                                                                                                                                                                                       >  1:452215/1‑201 (MQ=255)
                                                                        aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGcc                                                                                                                                                                  <  1:245566/251‑1 (MQ=255)
                                                                        aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCGGCTGGCggcgacaacccgcgcggc                                                                                                                                                                   >  2:312113/1‑234 (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                                                                                                                                                                                                                                         >  2:156575/4‑124 (MQ=255)
                                                                                                                                  gagcTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTgg                                                                                                                                                                                                                                         <  1:156575/121‑1 (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)
                                                                                                                                                                                                                        agTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTg                                                                                                                       <  1:152619/150‑1 (MQ=255)
                                                                                                                                                                                                                             cgcAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAagag            <  1:392901/251‑1 (MQ=255)
                                                                                                                                                                                                                               cAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTCGGGG‑AGTCACTCGGGCTGCCGGTGGAAAGTTTTATTTTTGCGACCCAAGTGAACGATAACGCGCCCCG‑TTTCCTGCACG‑ACGGGCAGGGGTAACCCAAAGCGACTCAGGAGACGTTATCCAACGCGAGGGACGTGGGTCAGACGGAAAACTGGCCCCgtgggggag              >  1:319905/1‑241 (MQ=255)
                                                                                                                                                                                                                                       gCAGGAGGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTccgcc  <  2:341423/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                                      
CCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGA‑AGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACG‑TTTCCTGCACG‑ACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCC  >  NC_000913/4205‑4685

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