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,4410GA92.5% 100.2 / ‑1.3 40E236E (GAG→GAAthrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (15/22);  total (16/24)
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.20e-01
*NC_0009134,4420AG90.0% 96.1 / ‑0.7 40T237A (ACG→GCG) thrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (14/22);  total (16/24)
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.70e-01

TCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGC  >  NC_000913/4200‑4684
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                  
tctAGCCACGAGTCAAAAACAGTCCACTGCAAGAAAAACTGTTCTGTACATTTGGCGGCAATATCGAAGCTTTTGCCTTCGACGGCGAGTTCGATGCCTGTCAGGCGCTGGTGAAGTAGGCGTTTTATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGAGTTGCTACTACGTTGAAGCTTTTGCGCAGCTGCCGCAGGAGACGCGCaa                                                                                                                                                                                                                                             <  2:185463/250‑1 (MQ=255)
          aGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTg                                                                                                                                                                                                                                  <  1:305933/251‑1 (MQ=255)
          aGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTg                                                                                                                                                                                                                                  <  1:81442/251‑1 (MQ=255)
             tacaaTCAGTACACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGCTGCCTGTCAGGCGCTGGCGAAGCAGACGTTTGATGATGAAGAGCTGAAAGAGGCGCTGGGGTAAAACTCAGCTAGCTCGATTAACATCCGCCGTTTGCTGGCGCAGTTTTGCTACTACTTTGAAGCAGTTGCGTAGCTGCCTCCGGAAGCGCGCAACCAGCTGGTTGtcc                                                                                                                                                                                                                               >  2:46827/4‑250 (MQ=255)
              aaaaTCAGTCCACTGCAAGAATAACTGTTCTGTACATGGGGCGGCAATATCGAAACTGTTGCCATCGTCGGCGATTTCGATGCCTGTCATGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAATCTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCTGGAAGCGCGCAACCAGCTGGTTGtctc                                                                                                                                                                                                                              <  1:297756/251‑1 (MQ=255)
                    aGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGcc                                                                                                                                                                                                                        <  1:48573/251‑1 (MQ=255)
                              aaGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAAc                                                                                                                                                                                                              >  1:254807/1‑251 (MQ=255)
                                 aaaaaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCAAAGCGGAAACTTc                                                                                                                                                                                                           <  1:250836/251‑1 (MQ=255)
                                  aaaaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCCGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCg                                                                                                                                                                                                          >  2:130479/1‑251 (MQ=255)
                                        gTTCTGTACATTGGGCGGCAATATCGAAACTGTGGCCAGCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGTCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGAt                                                                                                                                                                                                     <  1:23410/250‑1 (MQ=255)
                                          tCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCTCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTAGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCAGCAGGAAACGCGCAACCAGCTGGCTGACTCGGGGCCAAACggcacatccggcgattt                                                                                                                                                                                                   >  2:303599/1‑235 (MQ=255)
                                                    gggCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCggggt                                                                                                                                                                                        >  1:205800/1‑249 (MQ=255)
                                                     ggcggcAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCTGTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTct                                                                                                                                                                                       <  2:318297/251‑1 (MQ=255)
                                                      gctgaaagttCGAAACTGTTGCATTGGAAGGCGATTTCGATGCCTGTCAGGCGCTGGTGACGCAGTCGTTGGCGGATGAAGAGCTGAAAGTGGCGCTGGGGTTAACCTCTGCTAACGCTTTTAACAGCAGCCGTTTGCTGTCGCATATTTGCGACTCCTCTGAAGCAGTTGCGCATCTGCTGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCTGAAACTTCGGCGATTTGACGGCGGGTctg                                                                                                                                                                                      <  2:286160/242‑1 (MQ=255)
                                                                cGAAACTGTTGCCATCGGCGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCTGCGATTTGACGGCGGGTCTGCTGGCGAAGt                                                                                                                                                                            <  1:110911/251‑1 (MQ=255)
                                                                cGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCACCCAGCTGGTTGTCTCGGTGCAAAGCGGAACCTTCGGCGATTTGACGGCgggggaggcgggcgg                                                                                                                                                                              >  2:104392/1‑237 (MQ=255)
                                                                     cTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCggggcttcttggggaggcactc                                                                                                                                                                       >  1:124867/1‑232 (MQ=255)
                                                                             aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCTTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACGCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCTCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGcc                                                                                                                                                               <  1:344346/251‑1 (MQ=255)
                                                                             aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGcc                                                                                                                                                               <  1:278192/251‑1 (MQ=255)
                                                                               cGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAGGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCgg                                                                                                                                                             <  1:12810/251‑1 (MQ=255)
                                                                               cGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCg                                                                                                                                                              <  1:130479/250‑1 (MQ=255)
                                                                                 aCGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGTGAAGTCACTCGGTCTGCCGGTg                                                                                                                                                           <  1:287844/251‑1 (MQ=255)
                                                                                 aCGGCAATTTAGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCCGCTAACTCGATTAACATCAGCCGTTTTCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTTTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTg                                                                                                                                                           <  2:285941/251‑1 (MQ=255)
                                                                                    gCGATTTCGATGCCGAGCAGGCGGTGGGGAAGCAGGCGTTTGAGGAGTAAGAATTTAAAGTGGCTGTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTGGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCTAATTCACTCGGTCTGCCGGTgtaa                                                                                                                                                        <  2:376435/251‑4 (MQ=255)
                                                                                                                        cGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTg                                                                                                                    <  1:115642/251‑1 (MQ=255)
                                                                                                                                        gcTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCacg                                                                                                    <  1:272998/250‑1 (MQ=255)
                                                                                                                                            aaaGTGTCGCTGGGGTTAAACTCATCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATGTGCTACTACTTGGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCCCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCGGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTAtt                                                                                                                                               <  2:277142/204‑1 (MQ=255)
                                                                                                                                            aaaGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCCGGCTGCCGGGGAAACGGTTTAtt                                                                                                                                               >  1:277142/1‑204 (MQ=255)
                                                                                                                                                  gcgcTGGGGTAAAACTCAGCTAACGCGTTTAACATCAGCGGTTTGCTGGCGCAGTTTTGCTACTACTTTGAAGATGTGGCGGAGGTGCCGCAGGAAGCGCGCAACGAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGTCGATTTGACGGCGGGACTGCTGGCGAAGTCGCTCGGTCTGCCGGTGATACGTTTGATTGCTGCGACCATCGGGAACGATTCCGTGCCACGTGTCCTGCTCGACGGTCAGTg                                                                                          <  2:380182/251‑1 (MQ=255)
                                                                                                                                                                  cagCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTCCTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCGCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAATCGTTTTATTGCTGCGACCAACGTGAACGATGCCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACt                                                                          <  1:167880/249‑1 (MQ=255)
                                                                                                                                                                            gATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGCAGTAACTCGGGTTGCCGGGGAAAAGTTTTATTGCTGCGACCAACGTGAACGATAACGTGCCACGGTTCCTGCACGACGGTCAGTGGTCACCCCAAGCGAATCaagggccgt                                                                >  1:342669/1‑243 (MQ=255)
                                                                                                                                                                                  cATCCGCCGTTTGCTGGCGCAGATTTGCTTCTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCa                                                          <  1:287415/251‑1 (MQ=255)
                                                                                                                                                                                             tGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGGTGGCGAAGACACTCCGGCTGCCGGTGGAACGGTTTTTTTCTGCGACCAACGTGAAAGATACCGTGCCACGTTTCCTGCACGAAGGTCAGTGGTCACCCCAAGCCACTCAGGCGACGCTTTCCCACGCGATGGAc                                               >  1:285498/1‑251 (MQ=255)
                                                                                                                                                                                                                  actTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCCGGGGGAGGCCGTCGGGCTGCCGGGGAAACGTTTTATTTGTGCGACCAACGTGAACGATACCGTGCCACGTTTTCTGCACGACGGTCAGGGGTCACCCAAAACGACTACGGCGGCGTTATCCAACTCGATGGagggggtggaaccggacaaatgg                          >  1:272451/1‑229 (MQ=255)
                                                                                                                                                                                                                          gcagTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGCAGTCACTCCGGCTGCCGGGGAAACGTTTTATTTCTGCGACCAACATGAAAGCTACCGTGCCCCGCTTTCTGCACGACGGGCCGTGGGCACCCCAAACGACTCAGGCGACGTTATCCAACGCGATGGGCGTGGGGCAGCCGAACAACAGGCcgcgggg                  >  1:49562/4‑248 (MQ=255)
                                                                                                                                                                                                                           cagTTGCGCAGCTGCCGCAGGAAGCGCGCATCCAGCTGTTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGCCGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGGTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTAGCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGtgtg                 <  2:285498/249‑1 (MQ=255)
                                                                                                                                                                                                                            agTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTTCTGGGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGGCGGGCAGTGGTCACCCAAAACGACTCAGGCGACGTTATCCCAAGCGATGGACGTGGGGCAGCCGAACAAATGGCCGCgtgggg                >  1:118377/2‑248 (MQ=255)
                                                                                                                                                                                                                                      cTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGGAGGCACTCGGTCTGCCGGAGAAACGGTTTTTTTCTGCCGCCCACGATAACGATACCGTGCCACGGTTCCTGCACGACGGGCAATGGGCCCCCAAAGCGCCTCCAGCGACGTTATCCAACGCGGAGGGCGAGAGGCCACCAAACAACTGGCCGCGCGTGGGAGAgttgtt      >  1:264816/1‑251 (MQ=255)
                                                                                                                                                                                                                                       tGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCCGCTGGCGGAGGAACTCGGTCTGCCGGTGAAACGGTTTATTTCTGCGACCAACGTGGACGGAAACGGGCCACGGTTTCTTCCCGGCGGGCAGGGGGCACCCAAAACCGCCCAGGCGACGTTATTCAAACCGATGGGCGGGAGTCCGGCGCAAAACCGGCCGCgtgggggaggggtggtt     >  1:188350/1‑237 (MQ=255)
                                                                                                                                                                                                                                          cGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGGCACTCGGTCTGCCGGGGAAACGGTTTATTGCTGCGAACAACGTGAAAGATACCGGGCCACGTTTTCTGAACGGCGGGCCGTGGGCACCCAAAGCGGCTCAGGCGACGGTATCCAACGCGAAGGGCGGGAGGCAGCCGGACCACTGGACGCGTGTgggaggggtgtttcgc  >  1:255209/1‑237 (MQ=255)
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                  
TCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGC  >  NC_000913/4200‑4684

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