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

TACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTA  >  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)
            tgtggtgtCGCTTCCGTTTTGCTCGGGCCGTGGGCGGGTTATTGCGCCGGGTGGATTTTCGGTTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTGATGCTCAGTGCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGt                                                                                                                                                                                                              <  2:120395/222‑1 (MQ=255)
            aGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGATGGCGGCGATAGCGCGGCTGCTGACGCCCAATCTCGCCCCCGTGGAAATGGTTGGTGGGATggtggtgtggg                                                                                                                                                                                                               >  2:595934/1‑220 (MQ=255)
                 cgtcgcttgtggtgtgcgcGGGCGGTGGGCGGGAGATGTCCCCGGCGGGACTGTCGGGGTGTGCGGGGTTGGAACCGGTATGGCAGACGTGGTGGGGATCACGGCTTAGGCTCAGGTCTGGGTCCCCGATCTCGCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCg                                                                                                                                                                                                         <  2:474660/211‑1 (MQ=255)
                 cgtCGCTCCCTATTTTCCGGGGCCTTGGGCAGTATATTTAGCGGGCTGGGTTTCCTGGTTCTGCTTGGTTGGAACCGGTATGGCGGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCg                                                                                                                                                                                                         <  2:170579/227‑1 (MQ=255)
                 cTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGACCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCCATCTCGCCACCCTGACACATTTCTGTGGGATTGGGTTCTGGTTTGCg                                                                                                                                                                                                         >  1:507540/1‑229 (MQ=255)
                  ttCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGCCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGCACTCTCAGACCGGGTCGCCTCGCTGGCGGTGATTGTGCGGCTGCCGGCGCTCAATCCCGCCACCATGCACATGTTCTGGGAGAtagtgatcatgggtgcg                                                                                                                                                                                                          >  2:97865/1‑212 (MQ=255)
                                      gCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGAGGTGATAGGGCTGTTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTTGGTGGGATGGAGGTGTGGTTTGCCATGATCaaaaactgtgccctt                                                                                                                                                                                    >  1:564771/1‑218 (MQ=255)
                                                  atttttCACCGGCTGGACTTGCTGGTTCTGCTGGGTTGTAACCGGTATTGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGAt                                                                                                                                                                         <  1:215135/225‑1 (MQ=255)
                                                      ttCCCCGGCTGGATTTTCTGGTTCTTGTGGGTTGTAACCGGTATGGACGGCTTGGTGGCGATCCCGGCTTAGGCTCAGTTCTGGTTCCCCGATCTCTCCGCCGGGGTCGCCTCGCTGGCGGTGATAGGGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcg                                                                                                                                                                    <  2:356156/229‑1 (MQ=255)
                                                        cACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcgt                                                                                                                                                                   <  1:10777/228‑1 (MQ=255)
                                                            cgCTGGACTTACTGGTCCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGc                                                                                                                                                               <  1:532346/227‑1 (MQ=255)
                                                                           tttttttgggtttgtaCCGGTATGGCAGGTGGCGTGGGGATCCGGGGTTGTGCTCAGTTCCGGTCCCCCGTTCCCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGa                                                                                                                                               <  2:28394/214‑1 (MQ=255)
                                                                           ttCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGGTGACGCCCAATCTCGCCACCCTGAACATGTTCGGTGAGATGGGGTTCTGGTTTGTGGTGATCACAATACGCGCCCTTTGGTCGCTGATTgtcgggggaccggtacaggagggga                                                                                                                                               >  2:375125/1‑208 (MQ=255)
                                                                              ttctgttttttAAGCGGTATGGGTGACGTGTGGGGGGTCAGGGCTTATTCGCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTCCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATg                                                                                                                                             <  2:328238/218‑1 (MQ=255)
                                                                               gctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATTGGTCCTCCGCTGATCCCCAATCTTCCCCCCGTGAAAATGTTAGGGGGGATGGTGTGCTTGGTTGTGCTGGTCCAAATAGTAGCCCTTGGTTCGCGGGTTTTTGTTGGCCggggcgaggggggggggca                                                                                                                                           >  2:603258/1‑213 (MQ=255)
                                                                                ctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGCTAGTGCTGCTGCTGGCGCTCAATCTCGCCCCCCGGAAAATGGTCGGTGAGATGGAGGTCTGGTTTGCGCTGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGGGGTGCa                                                                                                                                           >  1:645332/1‑228 (MQ=255)
                                                                                     ttGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGTTTCCCCGATCTCTCCGACTGTGTCGCCTCGCTGGCGGTGATTGTGCTTCTGGTGCCACCCAATCTCGCCCCCCTGGAAATTTTTCGTGGGATGGGGGTCTGGTTTTCGGTGATCAAAATCGTCGCCATTGTGTCGCGGATTGTCGTGGGCCTGGTACTGGGGGCGGTGccatttga                                                                                                                                     >  2:756176/1‑222 (MQ=255)
                                                                                     ttGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCCCCACCGTGGAAAAGTTTGGGGGGGTGGGGTTCTGGGTTTGCATGGTTCAAAACCGACCCCTTGTGGTGGTGGATTTCGGCGGCCTGGTCATGGAGGcagggccatttta                                                                                                                                     >  2:91242/1‑217 (MQ=255)
                                                                                      tGTAACCGTGAGGGCAGGCTGGGTGGCGAGCACGGCTTCTGCGCAGTTCTGGTTCCCCGATCTCCCCGAATGGGGCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTAAATCCCGCCACCGTGAAGATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAg                                                                                                                                   <  2:462867/229‑1 (MQ=255)
                                                                                        tAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGCCTGAGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGTCCCTCAATCTCGCCCCCCGGAAAATGTTCGGTGAGTTGGGGTTTTGGTTTGCGGTGATCCAAATCGTTGCCAGTTTGTCCGTGGTTGGCGCGGGCCGGGTCactggggggggcagctttccttc                                                                                                                                  >  2:317666/1‑207 (MQ=255)
                                                                                              gttttgccgcggtgggtgCGTTCCCGGTTTTTGCGTAGTTCGGGTTCCCCGATCTCCCCGACTGGGTCGCCGCGCTGGCGGTGATAGTGGTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCGGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGAc                                                                                                                           <  2:380842/212‑1 (MQ=255)
                                                                                               tATGGCAGGCGTGGTGGGGTTCACGGTTTTTGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATTGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGAc                                                                                                                           <  2:370649/228‑1 (MQ=255)
                                                                                               tATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTTCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGGTAGTGCTGCCGCTGACGCTCAATCTCCCCACCGTGCAAATGTTCGGTGCGATGGAGGTATGGGTTGTGCTCATCAACATCGTAGCCATTGTGTCGCGGGTTTGCGTCGGCCGGGTTCAGGTGGGGGTGCACGTTCACTCGCCGAcc                                                                                                                          >  2:651333/1‑228 (MQ=255)
                                                                                               tATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGGGATTGTGCTGCAGCTGACGCCCACTCTCGCCCCCGGGAACAAGTTCGGTGGGATGGGGTTCTGGGTTGTGTGAATCAAAATCGGCGCCGTTGTGGTGCGTATTGCTGTCGGCCTGGTTCGGGGGGGGATGGCTTTTCAGTCACCGaat                                                                                                                          >  2:587690/1‑227 (MQ=255)
                                                                                                     aGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGtgtg                                                                                                                    <  1:457160/229‑1 (MQ=255)
                                                                                                      ggggtggtgGCGTACGCGGTTTATGCTCGGTTTTGGTTCCCCGAACTCTCCGACTGGTTCGCCTCGCTGGCGGTGATAGTGCTGCTGTTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTgg                                                                                                                   <  2:497643/226‑1 (MQ=255)
                                                                                                         gtggtgGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCCCAATCTCGCCACCGTGAAAATGTTCGGTGTGATGGAGTTCTGGTTTGCGCTGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCGCTTTCAGTCACCCCCTGGCGTGGaa                                                                                                                 >  1:686759/1‑228 (MQ=255)
                                                                                                         gtggcgGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGGTAGTTCCTCTGCTGACCCTCAATCTCCCCCCCGCGGAAATGTTCGGGGAGTTGGAGGTCTGGTTTTCGATGCTCATAATTCTGGCCATTTTGGCGCTGATTTTCGGCGGCCTGGTCATAGTTGGTATGTACttttcagccccgagcggggggggag                                                                                                                >  2:350834/1‑207 (MQ=255)
                                                                                                                   tctgggctattGCTCAGGTCTGGTTCCCCGATCCTCCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTTCTGACGAGCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCCTTTCAGTCGCCGACTGGTGTGGAAGCGTCATTcg                                                                                                       <  1:754984/218‑1 (MQ=255)
                                                                                                                     aCGGCTTATGCTCAGCTCTGGTTCCCCGATCGCTCCGACTGGGTCGCCTCGCTGGCGGGGATAGTGCTGCTGCTGACGCTCAATCTCGACACCGTGAAAAGGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTcgcg                                                                                                     <  1:471783/228‑1 (MQ=255)
                                                                                                                         cTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGGGATGGGGTTCTGGGTTGCGATGATCCAAATCGCCGCCATTGGTTCGCTGGTTgtcgt                                                                                                                                                                   >  2:347196/1‑163 (MQ=255)
                                                                                                                         cTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcgt                                                                                                                                                                   <  1:347196/163‑1 (MQ=255)
                                                                                                                              gCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGGGATGGAGTTCTGGTTTGCGATGATCCAAATCGTCGCCATTGTGTCGCTGATTGTCGTTGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGGGGAAGCGGCATTTGCGCATTtg                                                                                               >  1:545002/1‑225 (MQ=255)
                                                                                                                                 cAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCTGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGGGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATg                                                                                        >  1:227121/1‑229 (MQ=255)
                                                                                                                                         gTTCCCCGATCTCTCCGGCTGGGTCGCTTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACggcggc                                                                                <  1:236694/229‑1 (MQ=255)
                                                                                                                                          ttCCCCGATCTCCCCGAGTGGGTCGCCTCGCTGGCGGTGATAGGGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACggcggc                                                                                <  1:381323/228‑1 (MQ=255)
                                                                                                                                              ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGTTAGGGCTGCTGCTGACGCTCAATATCGGCCCCGTGAAAATGTTCGGGGAGATGGAGTTCTGGTTTGCGATGATCGACATCGTCGCCATTGTGTCGCTGCTTTGCCGCGGCCCGGTCATGGTGGCGATGCCGTTTCCGTCACCCACGGGTGGGGAGGCGTCTTTCGCGCTCTTGTGGAAAGAtcgcgggggggt                                                                           >  2:386895/1‑223 (MQ=255)
                                                                                                                                              ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTTTGTCGCCGATTGGCGTTGGcc                                                                                                                                                              >  2:107636/1‑147 (MQ=255)
                                                                                                                                              ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGcc                                                                                                                                                              <  1:107636/147‑1 (MQ=255)
                                                                                                                                                              ggTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGGTCGGTGAGATGGAGTTCTGGTTTGCGAGGATCAAAATCGTCGCCATTGTGTCGCTGATTTTCGTCGGCCTGGGCACAGTGGCGGTGCCCTTTTCTATCCCGCCCGGGTGGGAGAGCGCCTCTCGGCATTTTTGGGATGACGGGGGCTGGGTGCCGAGAGGGTTTAga                                                           >  2:283712/1‑228 (MQ=255)
                                                                                                                                                                         tGGCGGGTATAGTGCTGCTGCTGACGCTCAATCTCGCCAGCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGGCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTg                                                 <  1:317666/228‑1 (MQ=255)
                                                                                                                                                                           gCGGTGGTTGGGCTGGTGCTGACGCTGAAGCCCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGGCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGcc                                               <  1:91242/228‑1 (MQ=255)
                                                                                                                                                                                     tgctgctgcTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGGGGCTTCTTTGCCGGATTCCAGa                                     >  1:615335/1‑228 (MQ=255)
                                                                                                                                                                                       gtttttttgccgcttattcctgCCCCCTTGAAAAGTTTGGGTGGGGTGGCGTTCTGGTTGCTGTATATCAAGATCGCCGCGTTTGTGTCGAGTATTGTCGTCGGCCTGGTCATGGTGGGGATTGACTTTCAGTCCCCGACTGGTGTGGACGCGTCTTTCGCTCTTTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATa                                   <  2:486501/207‑1 (MQ=255)
                                                                                                                                                                                          ctgctgACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGGTTGCGATGATCAAAATCGTCGCCAGTGGATGGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCCTTTCAGTCACCGCCAGGTGTGGAATCGTCATTGGCGCATTTGTGGAATGACGGGGGCGGGTTACTCTAAGGGTTAATGCGCTTCTTTTCCGGGTTACAGTTatcag                               >  2:85898/1‑225 (MQ=255)
                                                                                                                                                                                             ctgACGCTCACTCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGTCGATGCACTTTCCGTCACCGACTGGTGGGGAAGCGCCATTCGCGCATTTGTGGAAGTACGCCGGATGGGTCCCCAAGGGGGTTGATGGCTTTTTTGCCGGGTACCAGaaatggggtt                            >  2:373434/1‑220 (MQ=255)
                                                                                                                                                                                               gACACTCAATCTCGCCACCGTGAAAATGTTCGGAGAGATGGCGTTCTGGTTTTCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGACATCGGCCTGGTCATGATGGCGATGCTCGTTCAGTCACCGACTGGTGGGGAAGCGTCATTCGCGCATTTGTGGCATGACGGCAGCCGGATGCCGGAAGGTATAAGCGTCGTCGTTTACGGAATCCaaaaagtgtgttta                          >  2:286974/1‑216 (MQ=255)
                                                                                                                                                                                                       atcccccacccGTGAAAATGGTCGGTGAAATTGAGGTCTGGTTTGCGATGATCAAAATCGGCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGt                  <  1:651334/220‑1 (MQ=255)
                                                                                                                                                                                                          tcGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGTTAGCGGTTTTCGCTTTCGTggg               >  1:84243/1‑229 (MQ=255)
                                                                                                                                                                                                                ccGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCGGGTTCCCGCAAGGTTTAAGTGGCTTCTTTTCCGGATTCCAGAAAGCGGTTTTCGCTTTTGGGGGGttttg         >  2:578275/1‑227 (MQ=255)
                                                                                                                                                                                                                 cGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAg        <  1:544812/229‑1 (MQ=255)
                                                                                                                                                                                                                 cGTGAAAATGGTCGGGGAGATGGATTTCTGGTTTGCAATTATCAAAGTCGTCGCCATTTGGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATGCGCGGATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAg        <  1:97865/229‑1 (MQ=255)
                                                                                                                                                                                                                    gAAAATGTTCGGTGGGATGGGGTTCTGGTTTGCGATGGTCAAAATCGTCGCCATTGTGTCGGTGATTGTCGTCGGCCTGGTCCTGGTGGCGATGCCCTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTg     <  2:645331/229‑1 (MQ=255)
                                                                                                                                                                                                                    gAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTg     <  1:10727/229‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTTTTGGGTTGGATGGGGTTTTGGTTTGCGATGATCAAAATCGTGGCCATTGTGGCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt   <  2:84243/229‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTGTTCGGTTAGATGGATGTCTGGTTTGCGATGATCAAAATCGTCGCTATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt   <  1:327577/229‑1 (MQ=255)
                                                                                                                                                                                                                       aaTGTTTGGTGGTGTGGGGTTCTGGTTTGCGATGACCAAAATCGTCGCCATTGTTTCGTTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggta  <  1:587691/229‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                           
TACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTA  >  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: