Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,430,340 A→G 20.3% E159G (GAG→GGG)  cycA → D‑alanine/D‑serine/glycine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,430,3400AG20.3% 111.4 / 3.4 59E159G (GAG→GGG) cycAD‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base A (15/28);  new base G (9/4);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 5.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTA  >  NC_000913/4430114‑4430557
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                           
tgcAAATCTTTTAGGGACTTCGCTTCCGATTTGCTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTgaga                                                                                                                                                                                                                           <  1:553016/227‑1 (MQ=255)
taaaaaatttttgggggcTTCGCTTCCGATTTACTCGGGCGGTGGGCAGGATATTTCACCGGCTGGCCTTACTGGGTCGGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTgaga                                                                                                                                                                                                                           <  2:114581/212‑1 (MQ=255)
  caattttttttgtgccttgtgttCCGATTTGCGTGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGGGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATg                                                                                                                                                                                                                         <  2:569632/208‑1 (MQ=255)
     aTCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGGCGCTCAATCTCGCCCCCGTGAACATGTTCGGGGCGATGGa                                                                                                                                                                                                                       >  1:28394/1‑228 (MQ=255)
            tgtggtgtCGCTTCCGTTTTGCTCGGGCCGTGGGCGGGTTATTGCGCCGGGTGGATTTTCGGTTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTGATGCTCAGTGCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGt                                                                                                                                                                                                              <  2:120395/222‑1 (MQ=255)
            aGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGATGGCGGCGATAGCGCGGCTGCTGACGCCCAATCTCGCCCCCGTGGAAATGGTTGGTGGGATggtggtgtggg                                                                                                                                                                                                               >  2:595934/1‑220 (MQ=255)
                 cgtcgcttgtggtgtgcgcGGGCGGTGGGCGGGAGATGTCCCCGGCGGGACTGTCGGGGTGTGCGGGGTTGGAACCGGTATGGCAGACGTGGTGGGGATCACGGCTTAGGCTCAGGTCTGGGTCCCCGATCTCGCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCg                                                                                                                                                                                                         <  2:474660/211‑1 (MQ=255)
                 cgtCGCTCCCTATTTTCCGGGGCCTTGGGCAGTATATTTAGCGGGCTGGGTTTCCTGGTTCTGCTTGGTTGGAACCGGTATGGCGGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCg                                                                                                                                                                                                         <  2:170579/227‑1 (MQ=255)
                 cTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGACCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCCATCTCGCCACCCTGACACATTTCTGTGGGATTGGGTTCT‑GGTTTGCg                                                                                                                                                                                                         >  1:507540/1‑229 (MQ=255)
                  ttCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGCCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGCACTCTCAGACCGGGTCGCCTCGCTGGCGGTGATTGTGCGGCTGCCGGCGCTCAATCCCGCCACCATGCACATGTTCTGGGAGAtagtgatcatgggtgcg                                                                                                                                                                                                          >  2:97865/1‑212 (MQ=255)
                                      gCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGAGGTGATAGGGCTGTTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTTGGTGGGATGGAGGTGT‑GGTTTGCCATGATCaaaaactgtgccctt                                                                                                                                                                                    >  1:564771/1‑218 (MQ=255)
                                                  atttttCACCGGCTGGACTTGCTGGTTCTGCTGGGTTGTAACCGGTATTGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGAt                                                                                                                                                                         <  1:215135/225‑1 (MQ=255)
                                                      ttCCCCGGCTGGATTTTCTGGTTCTTGTGGGTTGTAACCGGTATGGACGGCTTGGTGGCGATCCCGGCTTAGGCTCAGTTCTGGTTCCCCGATCTCTCCGCCGGGGTCGCCTCGCTGGCGGTGATAGGGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcg                                                                                                                                                                    <  2:356156/229‑1 (MQ=255)
                                                        cACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcgt                                                                                                                                                                   <  1:10777/228‑1 (MQ=255)
                                                            cgCTGGACTTACTGGTCCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGc                                                                                                                                                               <  1:532346/227‑1 (MQ=255)
                                                                           tttttttgggtttgtaCCGGTATGGCAGGTGGCGTGGGGATCCGGGGTTGTGCTCAGTTCCGGTCCCCCGTTCCCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGa                                                                                                                                               <  2:28394/214‑1 (MQ=255)
                                                                           ttCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGGTGACGCCCAATCTCGCCACCCTGAACATGTTCGGTGAGATGGGGTTCT‑GGTTTGTGGTGATCACAATACGCGCCCTTTGGTCGCTGATTgtcgggggaccggtacaggagggga                                                                                                                                               >  2:375125/1‑208 (MQ=255)
                                                                              ttctgttttttAAGCGGTATGGGTGACGTGTGGGGGGTCAGGGCTTATTCGCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTCCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATg                                                                                                                                             <  2:328238/218‑1 (MQ=255)
                                                                               gctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATTGGTCCTCCGCTGATCCCCAATCTTCCCCCCGTGAAAATGTTAGGGGGGATGGTGTGCT‑TGGTTGTGCTGGTCCAAATAGTAGCCCTTGGTTCGCGGGTTTTTGTTGGCCggggcgaggggggggggca                                                                                                                                           >  2:603258/1‑213 (MQ=255)
                                                                                ctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGCTAGTGCTGCTGCTGGCGCTCAATCTCGCCCCCCGGAAAATGGTCGGTGAGATGGAGGTCT‑GGTTTGCGCTGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGGGGTGCa                                                                                                                                           >  1:645332/1‑228 (MQ=255)
                                                                                     ttGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGTTTCCCCGATCTCTCCGACTGTGTCGCCTCGCTGGCGGTGATTGTGCTTCTGGTGCCACCCAATCTCGCCCCCCTGGAAATTTTTCGTGGGATGGGGGTCT‑GGTTTTCGGTGATCAAAATCGTCGCCATTGTGTCGCGGATTGTCGTGGGCCTGGTACTGGGGGCGGTGccatttga                                                                                                                                     >  2:756176/1‑222 (MQ=255)
                                                                                     ttGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCCCCACCGTGGAAAAGTTTGGGGGGGTGGGGTTCTGGGTTTGCATGGTTCAAAACCGAC‑CCCTTGTGGTGGTGGATTTCGGCGGCCTGGTCATGGAGGcagggccatttta                                                                                                                                     >  2:91242/1‑217 (MQ=255)
                                                                                      tGTAACCGTGAGGGCAGGCTGGGTGGCGAGCACGGCTTCTGCGCAGTTCTGGTTCCCCGATCTCCCCGAATGGGGCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTAAATCCCGCCACCGTGAAGATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAg                                                                                                                                   <  2:462867/229‑1 (MQ=255)
                                                                                        tAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGCCTGAGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGTCCCTCAATCTCGCCCCCCGGAAAATGTTCGGTGAGTTGGGGTTTT‑GGTTTGCGGTGATCCAAATCGTTGCCAGTTTGTCCGTGGTTGGCGCGGGCCGGGTCactggggggggcagctttccttc                                                                                                                                  >  2:317666/1‑207 (MQ=255)
                                                                                              gttttgccgcggtgggtgCGTTCCCGGTTTTTGCGTAGTTCGGGTTCCCCGATCTCCCCGACTGGGTCGCCGCGCTGGCGGTGATAGTGGTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCG‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGAc                                                                                                                           <  2:380842/212‑1 (MQ=255)
                                                                                               tATGGCAGGCGTGGTGGGGTTCACGGTTTTTGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATTGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGAc                                                                                                                           <  2:370649/228‑1 (MQ=255)
                                                                                               tATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTTCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGGTAGTGCTGCCGCTGACGCTCAATCTCCCCACCGTGCAAATGTTCGGTGCGATGGAGGTAT‑GGGTTGTGCTCATCAACATCGTAGCCATTGTGTCGCGGGTTTGCGTCGGCCGGGTTCAGGTGGGGGTGCAC‑GTTCACTCGCCGAcc                                                                                                                          >  2:651333/1‑228 (MQ=255)
                                                                                               tATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGGGATTGTGCTGCAGCTGACGCCCACTCTCGCCCCCGGGAACAAGTTCGGTGGGATGGGGTTCT‑GGGTTGTGTGAATCAAAATCGGCGCCGTTGTGGTGCGTATTGCTGTCGGCCTGGTTCGGGGGGGGATGGCT‑TTTCAGTCACCGaat                                                                                                                          >  2:587690/1‑227 (MQ=255)
                                                                                                     aGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGtgtg                                                                                                                    <  1:457160/229‑1 (MQ=255)
                                                                                                      ggggtggtgGCGTACGCGGTTTATGCTCGGTTTTGGTTCCCCGAACTCTCCGACTGGTTCGCCTCGCTGGCGGTGATAGTGCTGCTGTTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTgg                                                                                                                   <  2:497643/226‑1 (MQ=255)
                                                                                                         gtggtgGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCCCAATCTCGCCACCGTGAAAATGTTCGGTGTGATGGAGTTCT‑GGTTTGCGCTGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCGC‑TTTCAGTCACCCCCTGGCGTGGaa                                                                                                                 >  1:686759/1‑228 (MQ=255)
                                                                                                         gtggcgGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGGTAGTTCCTCTGCTGACCCTCAATCTCCCCCCCGCGGAAATGTTCGGGGAGTTGGAGGTCT‑GGTTTTCGATGCTCATAATTCTGGCCATTTTGGCGCTGATTTTCGGCGGCCTGGTCATAGTTGGTATGTAC‑ttttcagccccgagcggggggggag                                                                                                                >  2:350834/1‑207 (MQ=255)
                                                                                                                   tctgggctattGCTCAGGTCTGGTTCCCCGATCCTCCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTTCTGACGAGCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCC‑TTTCAGTCGCCGACTGGTGTGGAAGCGTCATTcg                                                                                                       <  1:754984/218‑1 (MQ=255)
                                                                                                                     aCGGCTTATGCTCAGCTCTGGTTCCCCGATCGCTCCGACTGGGTCGCCTCGCTGGCGGGGATAGTGCTGCTGCTGACGCTCAATCTCGACACCGTGAAAAGGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTcgcg                                                                                                     <  1:471783/228‑1 (MQ=255)
                                                                                                                         cTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGGGATGGGGTTCT‑GGGTTGCGATGATCCAAATCGCCGCCATTGGTTCGCTGGTTgtcgt                                                                                                                                                                   >  2:347196/1‑163 (MQ=255)
                                                                                                                         cTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcgt                                                                                                                                                                   <  1:347196/163‑1 (MQ=255)
                                                                                                                              gCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGGGATGGAGTTCT‑GGTTTGCGATGATCCAAATCGTCGCCATTGTGTCGCTGATTGTCGTTGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGGGGAAGCGGCATTTGCGCATTtg                                                                                               >  1:545002/1‑225 (MQ=255)
                                                                                                                                 cAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCTGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGGGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATg                                                                                        >  1:227121/1‑229 (MQ=255)
                                                                                                                                         gTTCCCCGATCTCTCCGGCTGGGTCGCTTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACggcggc                                                                                <  1:236694/229‑1 (MQ=255)
                                                                                                                                          ttCCCCGATCTCCCCGAGTGGGTCGCCTCGCTGGCGGTGATAGGGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACggcggc                                                                                <  1:381323/228‑1 (MQ=255)
                                                                                                                                              ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGTTAGGGCTGCTGCTGACGCTCAATATCGGCCCCGTGAAAATGTTCGGGGAGATGGAGTTCT‑GGTTTGCGATGATCGACATCGTCGCCATTGTGTCGCTGCTTTGCCGCGGCCCGGTCATGGTGGCGATGCCG‑TTTCCGTCACCCACGGGTGGGGAGGCGTCTTTCGCGCTCTTGTGGAAAGAtcgcgggggggt                                                                           >  2:386895/1‑223 (MQ=255)
                                                                                                                                              ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGGGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTTTGTCGCCGATTGGCGTTGGcc                                                                                                                                                              >  2:107636/1‑147 (MQ=255)
                                                                                                                                              ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGcc                                                                                                                                                              <  1:107636/147‑1 (MQ=255)
                                                                                                                                                              ggTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGGTCGGTGAGATGGAGTTCT‑GGTTTGCGAGGATCAAAATCGTCGCCATTGTGTCGCTGATTTTCGTCGGCCTGGGCACAGTGGCGGTGCCCTTTTCTATCCCGCCCGGGTGGGAGAGCGCCTCTCG‑GCATTTTTGGGATGACGGGGGCTGGGTGCCGAGAGGGTTTAga                                                           >  2:283712/1‑228 (MQ=255)
                                                                                                                                                                         tGGCGGGTATAGTGCTGCTGCTGACGCTCAATCTCGCCAGCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGGCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTg                                                 <  1:317666/228‑1 (MQ=255)
                                                                                                                                                                           gCGGTGGTTGGGCTGGTGCTGACGCTGAAGCCCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGGCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGcc                                               <  1:91242/228‑1 (MQ=255)
                                                                                                                                                                                     tgctgctgcTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGGGGCTTCTTTGCCGGATTCCAGa                                     >  1:615335/1‑228 (MQ=255)
                                                                                                                                                                                       gtttttttgccgcttattcctgCCCCCTTGAAAAGTTTGGGTGGGGTGGCGTTCT‑GGTTGCTGTATATCAAGATCGCCGCGTTTGTGTCGAGTATTGTCGTCGGCCTGGTCATGGTGGGGATTGAC‑TTTCAGTCCCCGACTGGTGTGGACGCGTCTTTCGCTCTTTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATa                                   <  2:486501/207‑1 (MQ=255)
                                                                                                                                                                                          ctgctgACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGGTTGCGATGATCAAAATCGTCGCCAGTGGATGGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCC‑TTTCAGTCACCGCCAGGTGTGGAATCGTCATTGGCGCATTTGTGGAATGACGGGGGCGGGTTACTCTAAGGGTTAATGCGCTTCTTTTCCGGGTTACAGTTatcag                               >  2:85898/1‑225 (MQ=255)
                                                                                                                                                                                             ctgACGCTCACTCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGTCGATGCAC‑TTTCCGTCACCGACTGGTGGGGAAGCGCCATTCGCGCATTTGTGGAAGTACGCCGGATGGGTCCCCAAGGGGGTTGATGGCTTTTTTGCCGGGTACCAGaaatggggtt                            >  2:373434/1‑220 (MQ=255)
                                                                                                                                                                                               gACACTCAATCTCGCCACCGTGAAAATGTTCGGAGAGATGGCGTTCT‑GGTTTTCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGACATCGGCCTGGTCATGATGGCGATGCTC‑GTTCAGTCACCGACTGGTGGGGAAGCGTCATTCGCGCATTTGTGGCATGACGGCAGCCGGATGCCGGAAGGTATAAGCGTCGTCGTTTACGGAATCCaaaaagtgtgttta                          >  2:286974/1‑216 (MQ=255)
                                                                                                                                                                                                       atcccccacccGTGAAAATGGTCGGTGAAATTGAGGTCT‑GGTTTGCGATGATCAAAATCGGCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGt                  <  1:651334/220‑1 (MQ=255)
                                                                                                                                                                                                          tcGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGTTAGCGGTTTTCGCTTTCGTggg               >  1:84243/1‑229 (MQ=255)
                                                                                                                                                                                                                ccGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCGGGTTCCCGCAAGGTTTAAGTGGCTTCTTTTCCGGATTCCAGAAAGCGGTTTTCGCTTTTGGGGGGttttg         >  2:578275/1‑227 (MQ=255)
                                                                                                                                                                                                                 cGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAg        <  1:544812/229‑1 (MQ=255)
                                                                                                                                                                                                                 cGTGAAAATGGTCGGGGAGATGGATTTCT‑GGTTTGCAATTATCAAAGTCGTCGCCATTTGGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATGCGCGGATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAg        <  1:97865/229‑1 (MQ=255)
                                                                                                                                                                                                                    gAAAATGTTCGGTGGGATGGGGTTCT‑GGTTTGCGATGGTCAAAATCGTCGCCATTGTGTCGGTGATTGTCGTCGGCCTGGTCCTGGTGGCGATGCCC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTg     <  2:645331/229‑1 (MQ=255)
                                                                                                                                                                                                                    gAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTg     <  1:10727/229‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTTTTGGGTTGGATGGGGTTTT‑GGTTTGCGATGATCAAAATCGTGGCCATTGTGGCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt   <  2:84243/229‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTGTTCGGTTAGATGGATGTCT‑GGTTTGCGATGATCAAAATCGTCGCTATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt   <  1:327577/229‑1 (MQ=255)
                                                                                                                                                                                                                       aaTGTTTGGTGGTGTGGGGTTCT‑GGTTTGCGATGACCAAAATCGTCGCCATTGTTTCGTTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggta  <  1:587691/229‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                           
TACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCT‑GGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTA  >  NC_000913/4430114‑4430557

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