Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,829,887 G→A A293A (GCC→GCT)  astA ← arginine N‑succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,829,8870GA96.4% 145.2 / NA 56A293A (GCC→GCT) astAarginine N‑succinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (24/30);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTG‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGA  >  NC_000913/1829642‑1830131
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
gccttgccATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑CATGGGGATAATTTTCCTTA‑GCGAc                                                                                                                                                                                                                                                  >  1:155612/1‑251 (MQ=255)
        aTAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCCCGGTTGCCGGATCGGGACGCACCAGCACCACGCGGG‑AATGGGGATAATTTTCATTA‑GCGCCccggcagg                                                                                                                                                                                                                                          >  1:259502/1‑251 (MQ=255)
             acCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTggg                                                                                                                                                                                                                                     >  1:123349/1‑251 (MQ=255)
                 cGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGGGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAAt                                                                                                                                                                                                                                 >  1:191175/1‑251 (MQ=255)
                          ccGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTgcg                                                                                                                                                                                                                        <  1:60794/251‑1 (MQ=255)
                                     ttCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGGTGCCGGATCGGGCCGCACCAGCACCACGAGGA‑AATGGGGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGcccc                                                                                                                                                                                                             >  1:536712/1‑250 (MQ=255)
                                               aTGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCCGCACCACGCGGA‑AATGGGGTTAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGACGGCTGCCCTTCTGACACtt                                                                                                                                                                                                   >  1:99519/1‑251 (MQ=255)
                                                        ggCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAg                                                                                                                                                                                            <  1:517895/249‑1 (MQ=255)
                                                             gTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCCCCAGCACCACGCGGA‑AATGGGGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAAACGCCCTGCGCCGGCTGCCCTTCTGCCACTTCCACCACCCGCCtt                                                                                                                                                                                     >  1:464246/1‑251 (MQ=255)
                                                                         ccGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCTGCTCGGGACGCCCCAGCACCACGCGGA‑AAAGGGGGTAAATTTTATTA‑GCGACCCAGCAGGCTGGGGAATCGGCCTGCCCAGGGTTCCCTTCTCCCACTTTCCCCAGCCGCCTTTTTCggcatgtg                                                                                                                                                                         >  2:130494/1‑245 (MQ=255)
                                                                          cGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGgcgc                                                                                                                                                                        <  1:141705/251‑1 (MQ=255)
                                                                              aaTCCATAAAGTCATGCTGTTTTCTCCTCTGCGAACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGGACGGTTGCCGGCTCGGAACGCACCAGCACCCCGCGGA‑ACTGGTGATAATTTTCATTA‑GCGGCCTGGCCAGCCGGGACATCGCCCCGCGCTGGCTGGCCTTCTTCCACTTCCCCCAACCGACTGTTACGGCTggtgccccc                                                                                                                                                                    >  1:228943/1‑244 (MQ=255)
                                                                                 ccATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGGGATAAATTTCATTA‑GCGACCAGGCAGGCTGGGAAAACGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACATTTAAGGCTGGGGCGCAcgcg                                                                                                                                                                 >  1:83011/1‑251 (MQ=255)
                                                                                      aaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCAGCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGa                                                                                                                                                            <  1:30245/251‑1 (MQ=255)
                                                                                                      tcctcTGCGCACAGGCGCACCAGACGAACGCGATCCCAGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTTATTA‑GCGACCAGGCAGGCGGGGAAATCGCCCTGCGCAGGGTGCCCTTCTGCCACTTCCACCAGCCGGCTTTTAAGGATGGCGCGCACGCgggcggattcacacttcagcgt                                                                                                                                            >  1:421439/1‑231 (MQ=255)
                                                                                                      ccctcTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGt                                                                                                                                           <  1:89698/250‑1 (MQ=255)
                                                                                                         tctGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGGCTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCgg                                                                                                                                        >  1:187957/1‑251 (MQ=255)
                                                                                                         tctGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑AGGCGGTTAAAATCCAACGATCGGTTGCCGGATAGGGACGCACCAGCAACACGCGGA‑CATGGTGCAAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCATTCTTCCACTTCCACAAGCCGGCTTTTACGGGTGGCGCGCACGCGGGCGC‑AGTCAAACTCAAGAGTCgg                                                                                                                                        >  1:247698/1‑251 (MQ=255)
                                                                                                         cctGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTTTG‑CGGCGGTTAAAATCAACCGCTCGGTTGCCGGATCGGTACGCACCAGCCCCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCgg                                                                                                                                        <  2:93785/250‑1 (MQ=255)
                                                                                                              gctcaGGCGCCCCAGTCGAACGCGACCCCCGGCTTGTCATTTGAGTGAAGCCAGTTGCG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCTACCAGGCAGGCTGGTAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCAc                                                                                                                                   <  2:439286/248‑1 (MQ=255)
                                                                                                              ccacaGGCGCTCCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCAGCCAGGTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCTCCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCAc                                                                                                                                   <  1:468076/250‑1 (MQ=255)
                                                                                                                        aCCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGAt                                                                                                                         <  1:513169/251‑1 (MQ=255)
                                                                                                                           agaggaatgcgATCCCCGGCTTGTCGTTTGAGGGCCTCCAGTTGTG‑CGGCGGTTAAAAGCAACCGCGCGGTTGCCGGATCGGTACGCACCAGCCCCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTc                                                                                                                      <  2:421439/243‑1 (MQ=255)
                                                                                                                            gACGTACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCg                                                                                                                     <  1:522082/251‑1 (MQ=255)
                                                                                                                                aaggcgATGCCCGGCGTGGCATTTGAGGGCATCAAGTTGGC‑CGGCGTGTAAAATCAAACTCTCGGTTGCCGGATCGGTACTCCCCAGCACCACGCGTA‑AATGGTTATAATTTTC‑TTGAGCGACCAGGCAGGCTGTGAAATCGCCCTGCGCGGGCTGCCCTTCTGCCACTTTCACCAGCCGTCGTTTACGGATGGCTCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGatat                                                                                                                 <  2:326360/248‑1 (MQ=255)
                                                                                                                                aaCGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGATGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCCCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTTGatat                                                                                                                 >  1:186438/1‑251 (MQ=255)
                                                                                                                                   gcgATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGGACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGGAATAGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTTCGGATGGCCCGCACGCGGTCGA‑TGGCACACTCAAGCGGCGGCCCACCGTCAAAGATGTCGCTATAGt                                                                                                              >  1:247722/1‑251 (MQ=255)
                                                                                                                                     gATCCCCGGGGGGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCTCCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTa                                                                                                            <  2:232563/251‑1 (MQ=255)
                                                                                                                                       cccccGGCGTGGCATTGGAGGGCATCCAGTTGTG‑CGGCTGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCCCCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACg                                                                                                          <  2:155612/250‑1 (MQ=255)
                                                                                                                                          ccGGGGTGGCATTTGAGGGCATCCAGTTGGG‑CGGCGGTTAAAATCAAACTCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGTTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTa                                                                                                       <  2:369622/251‑1 (MQ=255)
                                                                                                                                            ggCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCg                                                                                                                                         <  1:170753/215‑1 (MQ=255)
                                                                                                                                               gTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTAGGTTGACGGATCTGTCCGCACCAGTAACACGCGGATAATGAT‑ATCCTTTTCATTA‑GCGACCAGGCAGGGTGGGTAAGCGCCCTGCCGAGACTAACCTTCTGTCACTTCCGCCAGCCGACTTTTACGGAAGGCGCGCACGCGGTCGAATG‑AACCCTCAAGCGTCGGCCCACCCTACAAGATGTCGTTATAggttaggtaacgc                                                                                                   >  1:33779/1‑238 (MQ=255)
                                                                                                                                                 ggtatgtGAGTGCATCCAGT‑GTGGCGGCGGTGAAAATTATACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGTGGATAATTTTCATTA‑GCTCCCAGGCAGGCGGGGAAATCGCTCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGTCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGaaaa                                                                                                <  2:536712/245‑1 (MQ=255)
                                                                                                                                                    aTTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACct                                                                                             <  1:346752/251‑1 (MQ=255)
                                                                                                                                                        gggggCTTCCAGTTTTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATTGGGACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACTCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACcttct                                                                                          <  2:222070/248‑1 (MQ=255)
                                                                                                                                                        ggaggCATTCAGTTGTG‑CGTCGTTTAAAATCAACCGATCTGTTGACGGATCTGTAAGACCCAGCACCACGGGGA‑AATGGTGATTATTTTCATTA‑GCTATCAGGCAGGCTGGGAAATCTCCCTGCGCAGGCTGCCCTGCTGCCACTTCCCCCAGCCGTCTTTTACGGCTGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGCCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACcttctt                                                                                         <  2:284444/248‑1 (MQ=255)
                                                                                                                                                                              ggTTAAAATCAAAAGCTCGGTAGCCGGATCGGTACGCACCAGCACCACGCGGA‑CATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTAACACTCACGCGTCGGCCCACCGTCAAAAATATCGCTATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGc                                                                      >  2:355830/1‑249 (MQ=255)
                                                                                                                                                                                 tAAATTCAAACGCTCGGTTGCCGGATCGGTACTCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGGTGGGAAATCGACCGGCGCAGGCTTCCCTTCTGCCACGTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGgc                                                                  <  2:259502/250‑1 (MQ=255)
                                                                                                                                                                                  aaaaTCAAACGCTCTGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGgcgc                                                                >  1:438820/1‑251 (MQ=255)
                                                                                                                                                                                        aaaCGCTCGGTTGCCGGATCGGTCCTCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCATGCAGGCTGGGAAACCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACTCTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTg                                                          <  1:314283/251‑1 (MQ=255)
                                                                                                                                                                                            gCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCggg                                                      >  1:184281/1‑250 (MQ=255)
                                                                                                                                                                                                 gTTGCCGGTTCGGTACGCCCCAGCACCACGCGGA‑ACTGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCAGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCCCGAGGCCGA‑TGTCACACTCACGAGTCGGCCCGCCGTCAAAGATGTCGGTATAGATACGGTATCGGAAACCTTCTCTCACCAGAGCAGAGCGGGCAGGCGCGGGTTGAGGagttt                                                  >  2:3182/1‑246 (MQ=255)
                                                                                                                                                                                                   ggCCGGATCGGTACGCCCCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCt                                               <  2:147681/250‑1 (MQ=255)
                                                                                                                                                                                                       ggATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGAcc                                           <  1:278005/251‑1 (MQ=255)
                                                                                                                                                                                                            ggTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCCGCACCGCGCGGGCAGGCGCGGTTTGCGGGTGTACCtttacgatg                                       >  2:394265/1‑242 (MQ=255)
                                                                                                                                                                                                                         cccccgcgGA‑AATGGTGATAATTTTCTTTG‑GCGCCCAGGCAGGCTGGGAAATCTCCTTGCGTAGGCGGCCCTTCTGCCACTTCCCCCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGGCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTGCCTGACCGATGACGTCCTGGGCTTc                         <  2:524698/247‑1 (MQ=255)
                                                                                                                                                                                                                         cccacGTGGA‑AATGGTGATTATTTTTATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCCTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCGCCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCCCCATCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTc                         <  1:309845/250‑1 (MQ=255)
                                                                                                                                                                                                                          ccgcgcgGA‑AAGGTTGATAATTTTCATTA‑GCGACCGGGCCGGCGTGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGTCTTTTACGGATGGCGCGCCCGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTcc                        <  2:400929/248‑1 (MQ=255)
                                                                                                                                                                                                                           cacGCGGA‑AATGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCt                       <  1:365992/251‑1 (MQ=255)
                                                                                                                                                                                                                                  a‑aaTGGTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATa                 >  1:346212/1‑250 (MQ=255)
                                                                                                                                                                                                                                       ggTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAg            <  1:345452/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCg                                                        <  2:10544/206‑1 (MQ=255)
                                                                                                                                                                                                                                        gTGATAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCg                                                        >  1:10544/1‑206 (MQ=255)
                                                                                                                                                                                                                                           aTAATTTTCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCCCCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGACCTGGGCTTCCTGGGATAAAAAGTggg        >  2:445161/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aTTTTCATTA‑TCGACCAGGCAGGCGGGGAATTCGCCCTGCGCAGGCTGCCCTTCTGCCACTGCCACCAGCCGCCTTTTACGGATGGCGCGCAAGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCGTCCTGGGATAAAAAGTGGGtat     <  2:464246/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  tCATTA‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGa  <  1:345098/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
GCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTG‑CGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGA‑AATGGTGATAATTTTCATTG‑GCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGA‑TGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGA  >  NC_000913/1829642‑1830131

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