Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 588 G→A G84G (GGG→GGA thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135880GA96.4% 158.5 / ‑5.5 56G84G (GGG→GGAthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/22);  minor base T (2/0);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

AGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGGTTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCG  >  NC_000913/342‑820
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
ggtgtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCACGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTcc                                                                                                                                                                                                                                        >  1:316648/2‑251 (MQ=255)
   gtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCg                                                                                                                                                                                                                                      <  1:641201/250‑1 (MQ=255)
   gtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGc                                                                                                                                                                                                                                     >  1:390834/1‑251 (MQ=255)
    ttGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCt                                                                                                                                                                                                                                    >  1:408513/1‑251 (MQ=255)
                       cAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTtc                                                                                                                                                                                                                 >  1:211837/1‑251 (MQ=255)
                            gCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTTACGGGACTCGCCGCCGCCCAGCCGGGTTTCCCGCTGGCGCAATTGAAAACTTtcgtcgg                                                                                                                                                                                                           >  2:650196/1‑250 (MQ=255)
                                 tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGCTCAgg                                                                                                                                                                                                      >  2:237918/1‑251 (MQ=255)
                                     gAACGTTTTCGGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAAtt                                                                                                                                                                                                  <  1:237918/251‑1 (MQ=255)
                                         gTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTACCCCCGCCAAAATCACAAACCACCTGGTTGCGCTGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGAATTTTTGCCGAACTTTTGTCGGGACTCGCCGCCTCCCAGCCCGGATTCCCGCTGGCGCCATTGAAAACTTTCGTCGTTCAGGACTTTGcc                                                                                                                                                                                              >  1:113422/1‑251 (MQ=255)
                                            ttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCa                                                                                                                                                                                            <  1:130185/250‑1 (MQ=255)
                                                 gggtgtGCCTATATTCTGGAAAGCATTGCCAGGCAGGTGAAGGTGGCCACCGACCCTTCTGCCCCCGCCAAAATCACCAACCACCTGGGGGCGATGATTGAAAAAACCATGAGCGGCCAGGATGCTTTGCCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATaaa                                                                                                                                                                                     <  2:117738/246‑1 (MQ=255)
                                                       gCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATg                                                                                                                                                                                <  1:650196/250‑1 (MQ=255)
                                                         cGATATTCTGGACAGCGATGCCAGTCAGGAGCAGGTGGCCACCGTCATCTCTGCCCCCGCAAAAATCACCAACCACCTGGTGGCGATGATTGAAACAACCATTTGCGGCCACGATGCTTTACCCAATATCCACGATGCCGAACAAATTTTTGCCGAAGTTTTCACGGGGCTATCAGCCGCCCAGCCGGGATTCCAGCTAGAGATATTGTACATTTTCGTCTGACAGGAATTTGCCCAATATAAACATGTcc                                                                                                                                                                             >  2:75027/1‑251 (MQ=255)
                                                         cGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTcc                                                                                                                                                                             <  1:471157/251‑1 (MQ=255)
                                                           atttgcTGGAAAGCAAGGCCAGGCAGGGCCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCCCCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTg                                                                                                                                                                           <  2:625213/246‑1 (MQ=255)
                                                           atatTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTg                                                                                                                                                                           >  1:625213/1‑251 (MQ=255)
                                                                                        cgggTGGCCAACGTCCTCTCTAACCCCGCCCAAATCACCAACCACCTGGCAGAGATGATTGAAAAAACAATTAACGGCCAGGAGGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTTCCGAACTTTTGACGAGCCTCGAAGCAGCCACGCCGGGTTTCCCGCTGGCGAAATTGGAAACTTACGTCGATCCGGAATTTGCCAAAATAGAACATGTCCTGCATGGCTTTAGTTAGTTGAGGCAATGccc                                                                                                                                              >  2:203968/3‑251 (MQ=255)
                                                                                             ggCCACGGGCCTCTCTGCCCCCGCCAAAATCACCAACCGCCTGGTGGCGATGATTGAAAAAACCATGAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGGCGATCAGTAAGTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATa                                                                                                                                         <  1:639136/251‑1 (MQ=255)
                                                                                                                cccGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCATTATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGACATTGCCCAAATAAAACATGTCCTGCATGGCATTCGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGAt                                                                                                                      >  2:151073/1‑251 (MQ=255)
                                                                                                                 ccGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGAtt                                                                                                                     <  1:587779/251‑1 (MQ=255)
                                                                                                                  cGCCAAAATCACCAACCACCCGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTTCCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGGAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGAttt                                                                                                                    <  2:390834/251‑1 (MQ=255)
                                                                                                                           ccccaaccaCCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCg                                                                                                           <  2:12157/249‑1 (MQ=255)
                                                                                                                                accaccTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGCAttt                                                                                                                                                                                                 >  2:213375/1‑161 (MQ=255)
                                                                                                                                accaccTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAAttt                                                                                                                                                                                                 <  1:213375/161‑1 (MQ=255)
                                                                                                                                          ggCGATGATTGAAAAAACCACTAGCGGCCAGGAGGCTTTACCCAATATCAGCGATACCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCTCCCAGACGGGATTACCGCTGGCGCCATTGAAAACTTACGTCGCACAGGAATTTGCCCAAGTAAAACATGTCCTGCATGGCATTCGTTTGATGGTGCAGTGCCCAGATAGAATCAACGCTGCGCGGATTTGCAGTGACGAGaaacgtacgggt                                                                                             >  2:656303/1‑241 (MQ=255)
                                                                                                                                            cGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCCTGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGGTTTGCCGTGGCGAGAAAATGTCGCTCGcc                                                                                          >  2:15485/1‑251 (MQ=255)
                                                                                                                                               tgatTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCAtt                                                                                       >  1:557773/1‑251 (MQ=255)
                                                                                                                                                   tGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAgg                                                                                                                                                                                                      <  2:408152/137‑1 (MQ=255)
                                                                                                                                                   tGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAgg                                                                                                                                                                                                      >  1:408152/1‑137 (MQ=255)
                                                                                                                                                     aaaaaaCCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCCGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCAGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGcc                                                                                 >  1:456269/1‑251 (MQ=255)
                                                                                                                                                     aaaaaaCCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGcc                                                                                 <  1:629068/251‑1 (MQ=255)
                                                                                                                                                         aaCCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGCATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCACTTATGGCCGGCg                                                                             >  2:590559/1‑251 (MQ=255)
                                                                                                                                                          aCCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGt                                                                            >  1:170748/1‑251 (MQ=255)
                                                                                                                                                             aTTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTAtt                                                                         <  2:463083/251‑1 (MQ=255)
                                                                                                                                                               tAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGACGTAtt                                                                         >  1:452408/1‑249 (MQ=255)
                                                                                                                                                                aGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTAtt                                                                         >  1:423346/1‑248 (MQ=255)
                                                                                                                                                                aGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAg                                                                       <  2:557773/250‑1 (MQ=255)
                                                                                                                                                                aGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACCGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTTGTTTGTTGGGGAAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGGCGATCGCCATTATGGCCGGCggttt                                                                         >  2:287221/1‑244 (MQ=255)
                                                                                                                                                                     ccAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGTTTTGCCGTGGCGAGAACATGTCGATCGCCCTTATGGCCgggggctaagagaccg                                                                 >  2:21047/1‑237 (MQ=255)
                                                                                                                                                                      cAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAgcgcg                                                                <  1:236283/251‑1 (MQ=255)
                                                                                                                                                                                tACCCAATAGCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACaa                                                       <  2:175105/250‑1 (MQ=255)
                                                                                                                                                                                        atCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGCGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACACCGTTAccgt                                              >  1:427058/1‑248 (MQ=255)
                                                                                                                                                                                          cAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGtt                                             >  1:617819/1‑250 (MQ=255)
                                                                                                                                                                                           agctaTGCCGAACGTATTTATGCCGACCTTGTGACACACCTCGTCGCCGGCCAGCCGGAATTCCCGCTGGCGCAATTGAAAACTTCCGTCGATCATGAATTTTCCCAAAGAAAACATGTCCTCCATGGCACTAGTTTGTTGGGGCAGTGCCCGGATATCATCAACGCTGCGCTGATTTGCCGTGGCTATAAAATTTCGATCTCCATTCTGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTAt                                           <  1:51269/247‑1 (MQ=255)
                                                                                                                                                                                           aGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGttt                                                                                                                                                             <  2:394489/137‑1 (MQ=255)
                                                                                                                                                                                           aGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGACGCCGACCAGCCGGGATTCCCGCTGGCGCAATTGCAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGttt                                                                                                                                                             >  1:394489/1‑137 (MQ=255)
                                                                                                                                                                                               aTGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAACTAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTAccggtat                                           >  2:142722/1‑241 (MQ=255)
                                                                                                                                                                                                  ccGAACGTATTTTTGCCGAACTTTTGACGAGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCg                                    >  1:657323/1‑251 (MQ=255)
                                                                                                                                                                                                   cGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCgg                                   >  1:515278/1‑251 (MQ=255)
                                                                                                                                                                                                        gTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCCTTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGa                               >  1:467736/1‑250 (MQ=255)
                                                                                                                                                                                                         tATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGaa                              >  1:613826/1‑250 (MQ=255)
                                                                                                                                                                                                           tttttGACGAACTTTTGGCGGGACTAGCCGCCGCCCAGCCGGGCTTCCCGCTGGCGCAATTGAAAACTTTCGTAGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGAAGTGCCCGGCTAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGCTCGCCATTATGGCCGACGTATTAGACGCGGGAGGTCACACCGATACCGTTATCGCTCCgggcgcaaa                            >  2:1209/1‑243 (MQ=255)
                                                                                                                                                                                                                 ccGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAActgctg                     >  1:239395/1‑251 (MQ=255)
                                                                                                                                                                                                                           taCGGGACGCGCCGCCGCCCAGCCGGGAGTCCCGCGGGCTCAATTGAAAACTTGCGTCGAGCAGGAATTTGCCCAAATAAAACAGGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGAGAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTgggg            <  2:355054/249‑1 (MQ=255)
                                                                                                                                                                                                                                gACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTAc      >  1:123437/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cgccgccgcCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCg  <  1:151073/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
AGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGGTTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCG  >  NC_000913/342‑820

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