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

ACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACT  >  NC_000913/1828666‑1829113
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
acgtgttCCCTGGTTGAAATGGCCCGGCGCAAACGGGGTGTCGCCCCCCCATGCCTCCGAGTCGCTGCCATGACGAACCCCCCCGCTGGCGGGCCTTTCAAAAACCCGCCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAAcggc                                                                                                                                                                                                                               <  2:259192/222‑1 (MQ=255)
        gtgatttGAAATAGCCGGGGCGAAGGCGGGGGCCCTCCTCCCATCGCTCCCAGTTGCTGCCATGCAATCACCCCCTGTGGTGCGGCTTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACggcgggcg                                                                                                                                                                                                                         <  2:1270919/220‑1 (MQ=255)
           gCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACggcggggcccgggg                                                                                                                                                                                                                   >  1:1028527/1‑220 (MQ=255)
            cTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGCAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAACCGGCGGGCGAAGGAGc                                                                                                                                                                                                                  >  1:460374/1‑228 (MQ=255)
                      gcgcGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCCGACGTGGAGCGGTGAcc                                                                                                                                                                                                        >  2:155324/1‑228 (MQ=255)
                          ggcgCAAGCAGGGGTCGTCCGCCCATCGCTTCCATTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATTGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTCACCGGc                                                                                                                                                                                                     >  1:1727427/1‑227 (MQ=255)
                                                cccATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGGGCAGGGGCAGCGTTGGACGGCTGTGGCAAAAGCCGCATGaaag                                                                                                                                                                               >  2:1626665/1‑226 (MQ=255)
                                                cccATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGACCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGGGCACGGGCAGCGTTTACCGGCGTTGTCAAACGcccaatgtaag                                                                                                                                                                               >  2:1182988/1‑219 (MQ=255)
                                                            atttGCGGCCATGCATTACCCCGCTGCTGTGGCGTCTGTACAGAAATCCCCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACa                                                                                                                                                                   <  2:1606245/226‑1 (MQ=255)
                                                                             cacccccTGCTGTGCGGCCTGTCCAGAATCCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGc                                                                                                                                                  <  1:1629545/227‑1 (MQ=255)
                                                                                 accTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCGCCGTTTATTGCTGAAAAGCGGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCAcc                                                                                                                                              >  1:180440/1‑228 (MQ=255)
                                                                                      ctgTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGt                                                                                                                                                                             >  1:130165/1‑192 (MQ=255)
                                                                                       gtgggcgCCTTTTAAAAAAACAGCCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTc                                                                                                                                                                            <  2:129934/187‑1 (MQ=255)
                                                                                         tgCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTCACCGGCTGTGCCCAACGCCGACTGAAAGGTCAGATGGCCAGCCGCGCCGCAATCCGCCACCTCCTCgct                                                                                                                                      >  1:699399/1‑226 (MQ=255)
                                                                                          gCGGCCTGTCCAGAAAGCACCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTc                                                                                                                                                                            <  1:1364056/189‑1 (MQ=255)
                                                                                           cGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGGTGGACAGCCGCGTCGATATCCGCCACCTCATCg                                                                                                                                        >  1:1043786/1‑224 (MQ=255)
                                                                                           cGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGGGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGGCAGATGGACAGCCGCGTCGATATCCGCCACCCCATCGCTaa                                                                                                                                    >  1:1132050/1‑228 (MQ=255)
                                                                                             gCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGGGCAGCGTTGACCGGCTGTGGCAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCCCCTCATCGATAAtt                                                                                                                                  >  1:228308/1‑228 (MQ=255)
                                                                                                       aaaTCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAAcggc                                                                                                                                                                                                                               >  1:610385/1‑126 (MQ=255)
                                                                                                        aaTCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAAcggcg                                                                                                                                                                                                                             <  2:609391/126‑1 (MQ=255)
                                                                                                            aGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTg                                                                                                                                                                                   >  2:1476832/1‑164 (MQ=255)
                                                                                                            aGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCCGTGACCAACGCCCGCTGa                                                                                                                                                                                  >  2:41450/1‑165 (MQ=255)
                                                                                                            aGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGGGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGGCAGCCGCGTCGATATCCGCCACCTCATCGATAATTAACGGATTATTTCCg                                                                                                                   >  1:44338/1‑228 (MQ=255)
                                                                                                             gCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGaa                                                                                                                                                                                 <  1:41523/165‑1 (MQ=255)
                                                                                                             agaccGCCAATAAACGGCTGCGGTTCGTCATCAAAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGGGACAAAAGCCGACTGAATGGTCAGATGGAAAGCACGGGCGAGATCCGCGACCGCATCGAGTATTTGCGGATTATTACCg                                                                                                                   >  2:1111894/4‑228 (MQ=255)
                                                                                                               cccGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCGGGCGTTAATCGCTTGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTGTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGcc                                                                                                                 <  1:1956441/228‑1 (MQ=255)
                                                                                                                   ccAATAAGCGGCATCAGTTCGGCATCCCAGTTGCCCGGCGTTAACCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGcc                                                                                                                 <  2:312631/224‑1 (MQ=255)
                                                                                                                            ggCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCg                                                                                                    <  2:45072/228‑1 (MQ=255)
                                                                                                                               tGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGGGCAGCGTTGACCGGCTGTGACAAAAGCCGACTGAAAGGGCAGATGGGCAGCCGCGTCGACATACGCCCCCCCCACGGTAAATAGCGGGGTATTAATCCCCCTATCCAGGGCggga                                                                                                 >  2:1274784/1‑228 (MQ=255)
                                                                                                                                            cccAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAAAGGTCAGATGGGCAGCCGCGTCGATATCCGCCACCCCATCGATAATTTGCgg                                                                                                                             >  2:443122/1‑187 (MQ=255)
                                                                                                                                             ccAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGa                                                                                                                            <  1:443853/187‑1 (MQ=255)
                                                                                                                                                   gCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCGCCGTTTATTGCTGAAAAGCGGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTCTTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGTCCGga                                                                             >  1:1195394/1‑228 (MQ=255)
                                                                                                                                                       ggCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGcc                                                                                                                                                                                        >  2:371786/1‑117 (MQ=255)
                                                                                                                                                        gctttagccGCTGGCTGAGGGCAGCAAGACGAGCAAGAAACGCTCCGCCCTGCGCCCGCCGTTTATTGCTGAAAAGCGGGGCGCAGGGGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCt                                                                        <  2:1193252/220‑1 (MQ=255)
                                                                                                                                                        gCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCg                                                                                                                                                                                       <  1:372403/117‑1 (MQ=255)
                                                                                                                                                        gCGTGAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCt                                                                        <  1:600804/228‑1 (MQ=255)
                                                                                                                                                         cGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGagag                                                                          >  2:468406/1‑225 (MQ=255)
                                                                                                                                                              accggttgctgACCGCAACCAGCCGAGCAAGAAACGCCTCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGAt                                                                  <  2:1026660/221‑1 (MQ=255)
                                                                                                                                                               tCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGGGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCCCCTCATCGGTAATTAATGGGTTATTACCGCCCAACCCAAGGGCGAGAAATTTCTCCCGCTGCCCGgagggacggccgg                                                                  >  2:1868317/1‑217 (MQ=255)
                                                                                                                                                                   tggctgACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGGGCAGCGTTGTCCCGCTTGTACAAACGCCGACTGGATGGTCAGATGGCCAGCCGCGTCGGTATACGCCCCCCCCTCGATAATTACCGGGTTATTACCGCCCATCTCCAGCGGTGGAATTTTCTCCGGCCGACCGGCGAAGTGGCGATg                                                                 >  2:1429892/1‑224 (MQ=255)
                                                                                                                                                                     gctgACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCg                                                                                                                                                                                                       >  1:508464/1‑88 (MQ=255)
                                                                                                                                                                     gctgACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGa                                                                                                                                                                                  >  1:1683429/1‑109 (MQ=255)
                                                                                                                                                                      ctgACGGCAACCAGACGAGCAGGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCgg                                                                                                                                                                                                      <  2:507634/88‑1 (MQ=255)
                                                                                                                                                                      ctgACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGaa                                                                                                                                                                                 <  2:1680459/109‑1 (MQ=255)
                                                                                                                                                                                      gAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGgcgc                                          <  1:469173/228‑1 (MQ=255)
                                                                                                                                                                                      gAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGgcgc                                          <  2:1724374/228‑1 (MQ=255)
                                                                                                                                                                                       aGCAGGAACCGCATCGCCCTGCCCCCGCCGTTTATTGCTGAAAAGCGGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGgcgc                                          <  2:180129/227‑1 (MQ=255)
                                                                                                                                                                                         cAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTAc                                       <  1:1564088/228‑1 (MQ=255)
                                                                                                                                                                                              atcGCATCGCCCTGCGCCCCGCTGTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                  <  1:411596/226‑1 (MQ=255)
                                                                                                                                                                                               ccGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTa                                        <  1:1871624/220‑1 (MQ=255)
                                                                                                                                                                                                 cgctgcaccTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGGCAGATGGACAGCCGCGGCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAAATTTCTCCGGCTGACCGGAGAGATGACCATGCAACTGGTGGCCTGGATTGTCGCTACCGGTAAAc                               >  2:1449871/8‑228 (MQ=255)
                                                                                                                                                                                                           gcgcCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGGGACAAACGCCGACTGAATGGGCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGCGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATTCAAATGGTAGCCCGTATTGGGGCTACCCGGAAACAGCGAACCgg                     >  2:1057442/1‑227 (MQ=255)
                                                                                                                                                                                                             gcCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCg                   <  1:1277067/228‑1 (MQ=255)
                                                                                                                                                                                                              ccccGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACAGGGAGCCTGGATTGGGGCTACCGGGAAACAGCAAACCGTCg                   >  1:97872/1‑227 (MQ=255)
                                                                                                                                                                                                              ccccGCTTTTCAGCAATAAACGGCGGGCGCAGGGGCAGCGTTGACCGGCTGTGACAAACGCCCACTGAATGGTCAGATGGACAGCCGCGTCGATATCCCCCACCCCATCGATAAATa                                                                                                                                 >  2:597180/1‑117 (MQ=255)
                                                                                                                                                                                                              ccccGATTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCCCCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCg                   <  1:348702/227‑1 (MQ=255)
                                                                                                                                                                                                               cccGCTTTTCAGCAGTGAGCGGCGGGGGGGGGTGCAACGTTGGCGGGAGGTGACAAACGAGGTCTGAATGGTCAGGTGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAg                                                                                                                                <  1:598159/117‑1 (MQ=255)
                                                                                                                                                                                                                ccGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTa                                                                                                                                 >  1:55373/1‑115 (MQ=255)
                                                                                                                                                                                                                 cGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAg                                                                                                                                <  2:55275/115‑1 (MQ=255)
                                                                                                                                                                                                                 cGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCCAGGGCGAGAATTTTTTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGGGAAAAGCAAACCGTCgcga                >  2:363191/1‑227 (MQ=255)
                                                                                                                                                                                                                   gtttcCAGCAATAAACGGCGGGCGCAGGTGCCGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGAtc              <  1:479142/226‑1 (MQ=255)
                                                                                                                                                                                                                   cTTTTCAGCAATAAACGGCGGGCGCAGGTGCGGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGACACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGAtc              <  1:1740936/227‑1 (MQ=255)
                                                                                                                                                                                                                    ttttCAGAAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGAtcct            <  2:698228/228‑1 (MQ=255)
                                                                                                                                                                                                                      ttCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTTACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAg                                                       >  2:54621/1‑183 (MQ=255)
                                                                                                                                                                                                                       tCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGc                                                      <  1:54717/183‑1 (MQ=255)
                                                                                                                                                                                                                              tAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTa                                                                                                                                 >  1:1539171/1‑101 (MQ=255)
                                                                                                                                                                                                                              tAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACt  <  2:936185/228‑1 (MQ=255)
                                                                                                                                                                                                                               aaaCGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAg                                                                                                                                <  2:1536431/101‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
ACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAATAAACGGCGGGCGCAGGTGCAGCGTTGACCGGCTGTGACAAACGCCGACTGAATGGTCAGATGGACAGCCGCGTCGATATCCGCCACCTCATCGATAATTAGCGGATTATTACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACT  >  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: