Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,430,339 A→G 20.8% 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,3390AG20.8% 91.3 / 2.9 47E159G (GAG→GGG) cycAD‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base A (10/24);  new base G (7/3);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 3.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTA  >  NC_000913/4430113‑4430556
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                          
tgcAAATCTTTTAGGGACTTCGCTTCCGATTTGCTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTgaga                                                                                                                                                                                                                          <  1:553016/227‑1 (MQ=255)
     aTCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGGCGCTCAATCTCGCCCCCGTGAACATGTTCGGGGCGATGGa                                                                                                                                                                                                                      >  1:28394/1‑228 (MQ=255)
            tgtggtgtCGCTTCCGTTTTGCTCGGGCCGTGGGCGGGTTATTGCGCCGGGTGGATTTTCGGTTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTGATGCTCAGTGCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGt                                                                                                                                                                                                              <  2:120395/222‑1 (MQ=255)
            aGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGATGGCGGCGATAGCGCGGCTGCTGACGCCCAATCTCGCCCCCGTGGAAATGGTTGGTGGGATggtggtgtggg                                                                                                                                                                                                              >  2:595934/1‑220 (MQ=255)
                 cgtCGCTCCCTATTTTCCGGGGCCTTGGGCAGTATATTTAGCGGGCTGGGTTTCCTGGTTCTGCTTGGTTGGAACCGGTATGGCGGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCg                                                                                                                                                                                                         <  2:170579/227‑1 (MQ=255)
                 cTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGACCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCCATCTCGCCACCCTGACACATTTCTGTGGGATTGGGTTCTGGTTTGCg                                                                                                                                                                                                         >  1:507540/1‑229 (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)
                                                                              ttctgttttttAAGCGGTATGGGTGACGTGTGGGGGGTCAGGGCTTATTCGCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTCCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATg                                                                                                                                             <  2:328238/218‑1 (MQ=255)
                                                                                ctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGCTAGTGCTGCTGCTGGCGCTCAATCTCGCCCCCCGGAAAATGGTCGGTGAGATGGAGGTCTGGTTTGCGCTGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGGGGTGCa                                                                                                                                           >  1:645332/1‑228 (MQ=255)
                                                                                     ttGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGTTTCCCCGATCTCTCCGACTGTGTCGCCTCGCTGGCGGTGATTGTGCTTCTGGTGCCACCCAATCTCGCCCCCCTGGAAATTTTTCGTGGGATGGGGGTCTGGTTTTCGGTGATCAAAATCGTCGCCATTGTGTCGCGGATTGTCGTGGGCCTGGTACTGGGGGCGGTGccatttga                                                                                                                                     >  2:756176/1‑222 (MQ=255)
                                                                                      tGTAACCGTGAGGGCAGGCTGGGTGGCGAGCACGGCTTCTGCGCAGTTCTGGTTCCCCGATCTCCCCGAATGGGGCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTAAATCCCGCCACCGTGAAGATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAg                                                                                                                                   <  2:462867/229‑1 (MQ=255)
                                                                                               tATGGCAGGCGTGGTGGGGTTCACGGTTTTTGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATTGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGAc                                                                                                                           <  2:370649/228‑1 (MQ=255)
                                                                                               tATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTTCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGGTAGTGCTGCCGCTGACGCTCAATCTCCCCACCGTGCAAATGTTCGGTGCGATGGAGGTATGGGTTGTGCTCATCAACATCGTAGCCATTGTGTCGCGGGTTTGCGTCGGCCGGGTTCAGGTGGGGGTGCAC‑GTTCACTCGCCGAcc                                                                                                                          >  2:651333/1‑228 (MQ=255)
                                                                                               tATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGGGATTGTGCTGCAGCTGACGCCCACTCTCGCCCCCGGGAACAAGTTCGGTGGGATGGGGTTCTGGGTTGTGTGAATCAAAATCGGCGCCGTTGTGGTGCGTATTGCTGTCGGCCTGGTTCGGGGGGGGATGGCT‑TTTCAGTCACCGaat                                                                                                                          >  2:587690/1‑227 (MQ=255)
                                                                                                     aGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGtgtg                                                                                                                    <  1:457160/229‑1 (MQ=255)
                                                                                                      ggggtggtgGCGTACGCGGTTTATGCTCGGTTTTGGTTCCCCGAACTCTCCGACTGGTTCGCCTCGCTGGCGGTGATAGTGCTGCTGTTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTgg                                                                                                                   <  2:497643/226‑1 (MQ=255)
                                                                                                         gtggtgGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCCCAATCTCGCCACCGTGAAAATGTTCGGTGTGATGGAGTTCTGGTTTGCGCTGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCGC‑TTTCAGTCACCCCCTGGCGTGGaa                                                                                                                 >  1:686759/1‑228 (MQ=255)
                                                                                                                   tctgggctattGCTCAGGTCTGGTTCCCCGATCCTCCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTTCTGACGAGCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCC‑TTTCAGTCGCCGACTGGTGTGGAAGCGTCATTcg                                                                                                       <  1:754984/218‑1 (MQ=255)
                                                                                                                     aCGGCTTATGCTCAGCTCTGGTTCCCCGATCGCTCCGACTGGGTCGCCTCGCTGGCGGGGATAGTGCTGCTGCTGACGCTCAATCTCGACACCGTGAAAAGGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTcgcg                                                                                                     <  1:471783/228‑1 (MQ=255)
                                                                                                                         cTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGGGATGGGGTTCTGGGTTGCGATGATCCAAATCGCCGCCATTGGTTCGCTGGTTgtcgt                                                                                                                                                                   >  2:347196/1‑163 (MQ=255)
                                                                                                                         cTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcgt                                                                                                                                                                   <  1:347196/163‑1 (MQ=255)
                                                                                                                              gCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGGGATGGAGTTCTGGTTTGCGATGATCCAAATCGTCGCCATTGTGTCGCTGATTGTCGTTGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGGGGAAGCGGCATTTGCGCATTtg                                                                                               >  1:545002/1‑225 (MQ=255)
                                                                                                                                 cAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCTGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGGGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATg                                                                                        >  1:227121/1‑229 (MQ=255)
                                                                                                                                         gTTCCCCGATCTCTCCGGCTGGGTCGCTTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACggcggc                                                                                <  1:236694/229‑1 (MQ=255)
                                                                                                                                          ttCCCCGATCTCCCCGAGTGGGTCGCCTCGCTGGCGGTGATAGGGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACggcggc                                                                                <  1:381323/228‑1 (MQ=255)
                                                                                                                                              ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGTTAGGGCTGCTGCTGACGCTCAATATCGGCCCCGTGAAAATGTTCGGGGAGATGGAGTTCTGGTTTGCGATGATCGACATCGTCGCCATTGTGTCGCTGCTTTGCCGCGGCCCGGTCATGGTGGCGATGCCG‑TTTCCGTCACCCACGGGTGGGGAGGCGTCTTTCGCGCTCTTGTGGAAAGAtcgcgggggggt                                                                           >  2:386895/1‑223 (MQ=255)
                                                                                                                                              ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTTTGTCGCCGATTGGCGTTGGcc                                                                                                                                                              >  2:107636/1‑147 (MQ=255)
                                                                                                                                              ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGcc                                                                                                                                                              <  1:107636/147‑1 (MQ=255)
                                                                                                                                                              ggTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGGTCGGTGAGATGGAGTTCTGGTTTGCGAGGATCAAAATCGTCGCCATTGTGTCGCTGATTTTCGTCGGCCTGGGCACAGTGGCGGTGCCCTTTTCTATCCCGCCCGGGTGGGAGAGCGCCTCTCG‑GCATTTTTGGGATGACGGGGGCTGGGTGCCGAGAGGGTTTAga                                                           >  2:283712/1‑228 (MQ=255)
                                                                                                                                                                         tGGCGGGTATAGTGCTGCTGCTGACGCTCAATCTCGCCAGCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGGCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTg                                                 <  1:317666/228‑1 (MQ=255)
                                                                                                                                                                           gCGGTGGTTGGGCTGGTGCTGACGCTGAAGCCCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGGCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGcc                                               <  1:91242/228‑1 (MQ=255)
                                                                                                                                                                                     tgctgctgcTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGGGGCTTCTTTGCCGGATTCCAGa                                     >  1:615335/1‑228 (MQ=255)
                                                                                                                                                                                          ctgctgACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGGTTGCGATGATCAAAATCGTCGCCAGTGGATGGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCC‑TTTCAGTCACCGCCAGGTGTGGAATCGTCATTGGCGCATTTGTGGAATGACGGGGGCGGGTTACTCTAAGGGTTAATGCGCTTCTTTTCCGGGTTACAGTTatcag                               >  2:85898/1‑225 (MQ=255)
                                                                                                                                                                                             ctgACGCTCACTCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGTCGATGCAC‑TTTCCGTCACCGACTGGTGGGGAAGCGCCATTCGCGCATTTGTGGAAGTACGCCGGATGGGTCCCCAAGGGGGTTGATGGCTTTTTTGCCGGGTACCAGaaatggggtt                            >  2:373434/1‑220 (MQ=255)
                                                                                                                                                                                                       atcccccacccGTGAAAATGGTCGGTGAAATTGAGGTCTGGTTTGCGATGATCAAAATCGGCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGt                  <  1:651334/220‑1 (MQ=255)
                                                                                                                                                                                                          tcGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGTTAGCGGTTTTCGCTTTCGTggg               >  1:84243/1‑229 (MQ=255)
                                                                                                                                                                                                                ccGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCGGGTTCCCGCAAGGTTTAAGTGGCTTCTTTTCCGGATTCCAGAAAGCGGTTTTCGCTTTTGGGGGGttttg         >  2:578275/1‑227 (MQ=255)
                                                                                                                                                                                                                 cGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAg        <  1:544812/229‑1 (MQ=255)
                                                                                                                                                                                                                 cGTGAAAATGGTCGGGGAGATGGATTTCTGGTTTGCAATTATCAAAGTCGTCGCCATTTGGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATGCGCGGATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAg        <  1:97865/229‑1 (MQ=255)
                                                                                                                                                                                                                    gAAAATGTTCGGTGGGATGGGGTTCTGGTTTGCGATGGTCAAAATCGTCGCCATTGTGTCGGTGATTGTCGTCGGCCTGGTCCTGGTGGCGATGCCC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTg     <  2:645331/229‑1 (MQ=255)
                                                                                                                                                                                                                    gAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTg     <  1:10727/229‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTTTTGGGTTGGATGGGGTTTTGGTTTGCGATGATCAAAATCGTGGCCATTGTGGCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt   <  2:84243/229‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTGTTCGGTTAGATGGATGTCTGGTTTGCGATGATCAAAATCGTCGCTATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt   <  1:327577/229‑1 (MQ=255)
                                                                                                                                                                                                                       aaTGTTTGGTGGTGTGGGGTTCTGGTTTGCGATGACCAAAATCGTCGCCATTGTTTCGTTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggta  <  1:587691/229‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                          
TACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAC‑TTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTA  >  NC_000913/4430113‑4430556

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