Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,828,890 C→A 6.2% R282L (CGT→CTT)  astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,828,8900CA 6.2% 222.3 / 4.4 65R282L (CGT→CTT) astDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base A (2/2);  total (33/32)
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.17e-01

ACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGG‑GTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACT  >  NC_000913/1828666‑1829113
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
acgtgttCCCTGGTTGAAAT‑GGCCCGGCGCAAACGGGGTGTCGCCCCCCCATGCCTCCGAGTCGCTGCCATGACGAACCCCCCCGCTG‑GCGGGCCTTTCAAAAACCCGCCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAAcggc                                                                                                                                                                                                                               <  2:259192/222‑1 (MQ=255)
        gtgatttGAAATAGCCGGGGCGAAGGCGGGG‑GCCCTCCTCCCATCGCTCCCAGTTGCTGCCATGCAATCACCCCCTGTGGTGCGGCTTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACggc‑gggcg                                                                                                                                                                                                                         <  2:1270919/220‑1 (MQ=255)
           gCTTGTAATAAGCGCGGCGCAAGCAGGG‑GTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACggc‑ggggcccgggg                                                                                                                                                                                                                   >  1:1028527/1‑220 (MQ=255)
            cTTGTAATAAGCGCGGCGCAAGCAGGG‑GTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGCAACCACCTGCTGTGCGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAACCGGC‑GGGCGAAGGAGc                                                                                                                                                                                                                  >  1:460374/1‑228 (MQ=255)
                      gcgcGGCGCAAGCAGGG‑GTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCCGACGTGGAGCGGTGAcc                                                                                                                                                                                                        >  2:155324/1‑228 (MQ=255)
                          ggcgCAAGCAGGG‑GTCGTCCGCCCATCGCTTCCATTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATTGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTCACCGGc                                                                                                                                                                                                     >  1:1727427/1‑227 (MQ=255)
                                                cccATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGGGCAGGGGCAGCGTTGGACGGCTGTGGCAAAAGCCGCATGaaag                                                                                                                                                                               >  2:1626665/1‑226 (MQ=255)
                                                cccATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAG‑ACCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGGGCACGGGCAGCGTTTACCGGCGTTGTCAAACGcccaatgtaag                                                                                                                                                                               >  2:1182988/1‑219 (MQ=255)
                                                            atttGCGGCCATGCATTACCCCGCTGCTGTGGCGTCTGTACAGAAATCCC‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACa                                                                                                                                                                   <  2:1606245/226‑1 (MQ=255)
                                                                             cacccccTGCTGTGCGGCCTGTCCAGAATCCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGc                                                                                                                                                  <  1:1629545/227‑1 (MQ=255)
                                                                                 accTGCTGTGCGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCGCCGTTTATTGCTGAAAAGCG‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCAcc                                                                                                                                              >  1:180440/1‑228 (MQ=255)
                                                                                      ctgTGCGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGt                                                                                                                                                                             >  1:130165/1‑192 (MQ=255)
                                                                                       gtgggcgCCTTTTAAAAAAACAGCCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTc                                                                                                                                                                            <  2:129934/187‑1 (MQ=255)
                                                                                         tgCGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTCACCGGCTGTGCCCAACGCCGACTGAAAGGTCAGATGGCCAGCCGCGCCGCAATCCGCCACCTCCTCgct                                                                                                                                      >  1:699399/1‑226 (MQ=255)
                                                                                          gCGGCCTGTCCAGAAAGCAC‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTc                                                                                                                                                                            <  1:1364056/189‑1 (MQ=255)
                                                                                           cGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGGTGGACAGCCGCGTCGATATCCGCCACCTCATCg                                                                                                                                        >  1:1043786/1‑224 (MQ=255)
                                                                                           cGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGGGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGGCAGATGGACAGCCGCGTCGATATCCGCCACCCCATCGCTaa                                                                                                                                    >  1:1132050/1‑228 (MQ=255)
                                                                                             gCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGGGCAGCGTTGACCGGCTGTGGCAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCCCCTCATCGATAAtt                                                                                                                                  >  1:228308/1‑228 (MQ=255)
                                                                                                       aaaTCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAAcggc                                                                                                                                                                                                                               >  1:610385/1‑126 (MQ=255)
                                                                                                        aaTCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAAcggc‑g                                                                                                                                                                                                                             <  2:609391/126‑1 (MQ=255)
                                                                                                            aG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTg                                                                                                                                                                                   >  2:1476832/1‑164 (MQ=255)
                                                                                                            aG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCCGTGACCAACGCCCGCTGa                                                                                                                                                                                  >  2:41450/1‑165 (MQ=255)
                                                                                                            aG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGGGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGGCAGCCGCGTCGATATCCGCCACCTCATCGATAATTAACGGATTATTTCCg                                                                                                                   >  1:44338/1‑228 (MQ=255)
                                                                                                             g‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGaa                                                                                                                                                                                 <  1:41523/165‑1 (MQ=255)
                                                                                                             agaccGCCAATAAACGGCTGCGGTTCGTCATCAAAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGGGACAAAAGCCGACTGAATGGTCAGATGGAAAGCACGGGCGAGATCCGCGACCGCATCGAGTATTTGCGGATTATTACCg                                                                                                                   >  2:1111894/4‑228 (MQ=255)
                                                                                                               cccGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCGGGCGTTAATCGCTTGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTGTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGcc                                                                                                                 <  1:1956441/228‑1 (MQ=255)
                                                                                                                   ccAATAAGCGGCATCAGTTCGGCATCCCAGTTGCCCGGCGTTAACCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGcc                                                                                                                 <  2:312631/224‑1 (MQ=255)
                                                                                                                            ggCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCg                                                                                                    <  2:45072/228‑1 (MQ=255)
                                                                                                                               tGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGGGCAGCGTTGACCGGCTGTGACAAAAGCCGACTGAAAGGGCAGATGGGCAGCCGCGTCGACATACGCCCCCCCCACGGTAAATAGCGGGGTATTAATCCCCCTATCCAGGGCggga                                                                                                 >  2:1274784/1‑228 (MQ=255)
                                                                                                                                            cccAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAAAGGTCAGATGGGCAGCCGCGTCGATATCCGCCACCCCATCGATAATTTGCgg                                                                                                                             >  2:443122/1‑187 (MQ=255)
                                                                                                                                             ccAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGa                                                                                                                            <  1:443853/187‑1 (MQ=255)
                                                                                                                                                   gCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCGCCGTTTATTGCTGAAAAGCG‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTCTTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGTCCGga                                                                             >  1:1195394/1‑228 (MQ=255)
                                                                                                                                                       ggCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGcc                                                                                                                                                                                        >  2:371786/1‑117 (MQ=255)
                                                                                                                                                        gctttagccGCTGGCTGAGGGCAGCAAGACGAGCAAGAAACGCTCCGCCCTGCGCCCGCCGTTTATTGCTGAAAAGCG‑GGGCGCAGGGGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCt                                                                        <  2:1193252/220‑1 (MQ=255)
                                                                                                                                                        gCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCg                                                                                                                                                                                       <  1:372403/117‑1 (MQ=255)
                                                                                                                                                        gCGTGAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCt                                                                        <  1:600804/228‑1 (MQ=255)
                                                                                                                                                         cGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGagag                                                                          >  2:468406/1‑225 (MQ=255)
                                                                                                                                                              accggttgctgACCGCAACCAGCCGAGCAAGAAACGCCTCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGAt                                                                  <  2:1026660/221‑1 (MQ=255)
                                                                                                                                                               tCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGGGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCCCCTCATCGGTAATTAATGGGTTATTACCGCCCAACCCAAGGGCGAGAAATTTCTCCCGCTGCCCGgagggacggccgg                                                                  >  2:1868317/1‑217 (MQ=255)
                                                                                                                                                                   tggctgACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGGGCAGCGTTGTCCCGCTTGTACAAACGCCGACTGGATGGTCAGATGGCCAGCCGCGTCGGTATACGCCCCCCCCTCGATAATTACCGGGTTATTACCGCCCATCTCCAGCGGTGGAATTTTCTCCGGCCGACCGGCGAAGTGGCGATg                                                                 >  2:1429892/1‑224 (MQ=255)
                                                                                                                                                                     gctgACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCg                                                                                                                                                                                                       >  1:508464/1‑88 (MQ=255)
                                                                                                                                                                     gctgACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGa                                                                                                                                                                                  >  1:1683429/1‑109 (MQ=255)
                                                                                                                                                                      ctgACGGCAACCAGACGAGCAGGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCgg                                                                                                                                                                                                      <  2:507634/88‑1 (MQ=255)
                                                                                                                                                                      ctgACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGaa                                                                                                                                                                                 <  2:1680459/109‑1 (MQ=255)
                                                                                                                                                                                      gAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGgcgc                                          <  1:469173/228‑1 (MQ=255)
                                                                                                                                                                                      gAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGgcgc                                          <  2:1724374/228‑1 (MQ=255)
                                                                                                                                                                                       aGCAGGAACCGCATCGCCCTGC‑CCCCGCCGTTTATTGCTGAAAAGCGGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGgcgc                                          <  2:180129/227‑1 (MQ=255)
                                                                                                                                                                                         cAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTAc                                       <  1:1564088/228‑1 (MQ=255)
                                                                                                                                                                                              atcGCATCGCCCTGCGCCCCGCTGTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                  <  1:411596/226‑1 (MQ=255)
                                                                                                                                                                                               ccGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTa                                        <  1:1871624/220‑1 (MQ=255)
                                                                                                                                                                                                 cgctgcaccTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGGCAGATGGACAGCCGCGGCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAAATTTCTCCGGCTGACCGGAGAGATGACCATGCAACTGGTGGCCTGGATTGTCGCTACCGGTAAAc                               >  2:1449871/8‑228 (MQ=255)
                                                                                                                                                                                                           gcgcCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGGGACAAACGCCGACTGAATGGGCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGCGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATTCAAATGGTAGCCCGTATTGGGGCTACCCGGAAACAGCGAACCgg                     >  2:1057442/1‑227 (MQ=255)
                                                                                                                                                                                                             gcCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCg                   <  1:1277067/228‑1 (MQ=255)
                                                                                                                                                                                                              ccccGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACAGGGAGCCTGGATTGGGGCTACCGGGAAACAGCAAACCGTCg                   >  1:97872/1‑227 (MQ=255)
                                                                                                                                                                                                              ccccGCTTTTCAGCAATAAACGGC‑GGGCGCAGGGGCAGCGTTGACCGGCTGTGACAAACGCCCACTGAATGGTCAGATGGACAGCCGCGTCGATATCCCCCACCCCATCGATAAATa                                                                                                                                 >  2:597180/1‑117 (MQ=255)
                                                                                                                                                                                                              ccccGATTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCCCCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCg                   <  1:348702/227‑1 (MQ=255)
                                                                                                                                                                                                               cccGCTTTTCAGCAGTGAGCGGC‑GGGGGGGGGTGCAACGTTGGCGGGAGGTGACAAACGAGGTCTGAATGGTCAGGTGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAg                                                                                                                                <  1:598159/117‑1 (MQ=255)
                                                                                                                                                                                                                ccGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTa                                                                                                                                 >  1:55373/1‑115 (MQ=255)
                                                                                                                                                                                                                 cGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAg                                                                                                                                <  2:55275/115‑1 (MQ=255)
                                                                                                                                                                                                                 cGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCCAGGGCGAGAATTTTTTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGGGAAAAGCAAACCGTCgcga                >  2:363191/1‑227 (MQ=255)
                                                                                                                                                                                                                   gtttcCAGCAATAAACGGC‑GGGCGCAGGTGCCGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGAtc              <  1:479142/226‑1 (MQ=255)
                                                                                                                                                                                                                   cTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCGGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGACACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGAtc              <  1:1740936/227‑1 (MQ=255)
                                                                                                                                                                                                                    ttttCAGAAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGAtcct            <  2:698228/228‑1 (MQ=255)
                                                                                                                                                                                                                      ttCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTTACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAg                                                       >  2:54621/1‑183 (MQ=255)
                                                                                                                                                                                                                       tCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGc                                                      <  1:54717/183‑1 (MQ=255)
                                                                                                                                                                                                                              tAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTa                                                                                                                                 >  1:1539171/1‑101 (MQ=255)
                                                                                                                                                                                                                              tAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACt  <  2:936185/228‑1 (MQ=255)
                                                                                                                                                                                                                               aaaCGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAg                                                                                                                                <  2:1536431/101‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
ACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGG‑GTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAG‑CCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGC‑GGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACT  >  NC_000913/1828666‑1829113

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