Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,430,345 A→G 17.1% E161G (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,3450AG17.1% 147.9 / 2.7 72E161G (GAG→GGG) cycAD‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base A (22/35);  new base G (11/3);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCA  >  NC_000913/4430110‑4430578
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
gattACAAATTTTTTAGGGACTTCGCTTCCGATTTTCTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGtt                                                                                                                                                                                                                                          <  2:786454/236‑1 (MQ=255)
          cTTTTAGTGACTTCGTTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCg                                                                                                                                                                                                                                <  1:49592/239‑1 (MQ=255)
             ttAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGCGGTAGTGCTGCTGCTGACGCTCAATCTCGCCCCCCTGCAAAATGTCGGTGAGATGGAGATCTGGTTTGCgctt                                                                                                                                                                                                                             >  1:171723/1‑236 (MQ=255)
                        gtttcggtttttCTGGGGGCGTGGGCAGGATATTTCCCCGGCTGGGCTTACTTGTTCTGCTGGGTTGTAACGGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTATGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAtcgt                                                                                                                                                                                                                  <  2:577549/231‑1 (MQ=255)
                        gcctctgATTTGCTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGGAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAtcgt                                                                                                                                                                                                                  <  1:1035868/233‑1 (MQ=255)
                        gCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAtcgt                                                                                                                                                                                                                  <  1:77247/239‑1 (MQ=255)
                                 ttACTCGGGCCGTGGGCAGGATATTTCACCGGCTTGACTTACTGGTTCTGCTTGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATtgt                                                                                                                                                                                                         <  1:773782/239‑1 (MQ=255)
                                                 cAGGATATTTCCCCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcgt                                                                                                                                                                                          <  1:100527/238‑1 (MQ=255)
                                                 cAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcgt                                                                                                                                                                                          <  1:399947/238‑1 (MQ=255)
                                                         ttCCCCGGCTGGACTTACTGGGTCTGCTGGGTTGTAACCGGTATGGCAGGCGGTGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGt                                                                                                                                                                                 <  1:610035/239‑1 (MQ=255)
                                                              cGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTTGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGTTAGTTCTGCTGCTGACGCTCCATCTCGCCCCCCTGGAATTTTTTGGGGGGGTGGGGTTCTGGTTTGCGATGATCACAACGGTAGCCCTTTTGTCGCTGATTGTTGTTggagtgggacagg                                                                                                                                                                            >  2:55107/1‑228 (MQ=255)
                                                              cGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCCTTGTGTCGCTGATTGTCGTGGGCCTGGTCAtgg                                                                                                                                                                            >  1:56209/1‑239 (MQ=255)
                                                                               tCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGCTCAAAATCGTCGCCATTGTGTCGCTGATTTTCGGCGGCCTGTTCCTGGTGGCGATGCACTTTCa                                                                                                                                                            >  1:55143/1‑238 (MQ=255)
                                                                                ctgctgGGTTGTAACCGGTTGGGCAGAAGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAg                                                                                                                                                           <  1:737701/238‑1 (MQ=255)
                                                                                   ctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCCCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTAGGTTGCGCTGATCAAAAACGTCGCCATTGTGTCGCTGATTTTCGTCCGCCTGGTCCTGGTGGCTGTGCAATTTCAGTCa                                                                                                                                                        >  1:766481/1‑238 (MQ=255)
                                                                                   ctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGGGCTGCTGCTGACGCTCAATCTCGCCCCCGCGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCCTTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGGTGCCATTTCAggaa                                                                                                                                                        >  1:183863/1‑235 (MQ=255)
                                                                                   acaggTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCa                                                                                                                                                        <  1:64296/235‑1 (MQ=255)
                                                                                       gTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAACGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGGTGGACTTTCAGTCACCGAc                                                                                                                                                  >  1:687327/1‑239 (MQ=255)
                                                                                           tAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGATGCTGACGCCCAATCTCGCCCCCGTGAAAATGTTCGGTGAGTTGGGGGTCTGGGTTGCGATGATCAAAATGGTCGCCCTTGTTGCGCGGTCTTTGGCAGGCCGGGGTATGGGGGGGGAGCACTTTCAGTCACCGACgggg                                                                                                                                              >  2:105589/1‑238 (MQ=255)
                                                                                                ggTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGCCTGGGTCGCCTCGCTGGCGGTGAGAGTGCTGCTGCTGCCCCCCAATCTCGCCACCGCGGAAATCTTCGGGGTGGCGGGGTTTTGGTTTGCCGTTATCAAAATCGTCCCCCTTGTTAAGGCGGTTGTCGTTGTCCTGGGCATGGTGGGAGGGGACTTTTCGTCCCCGACTgggggggg                                                                                                                                         >  2:419420/1‑233 (MQ=255)
                                                                                                  tATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGGGAAAATGTTCGGTGAGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTTGTCGGCCTGGTCATGGGGGCGGTGCACTTTTAGTCACCGACTggggggggag                                                                                                                                       >  1:468767/1‑231 (MQ=255)
                                                                                                            gtggtgGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTTTCGTCGGCCTGGTCATGGTGGGGATGCACTTTCATTCACCGACTGGTGTGGAAGGGGCATTggc                                                                                                                             >  1:525683/1‑239 (MQ=255)
                                                                                                            gtggtgGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTTTCGTCGGCCTTGGCATGGTGGCGATGCACTTTCAGTCACCGCCTGGGGTGGAAGCGTCATTccc                                                                                                                             >  1:313112/1‑237 (MQ=255)
                                                                                                              ggtggCGATCCCGGCTTTTGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTcgcg                                                                                                                            <  1:1083437/238‑1 (MQ=255)
                                                                                                                tggCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACCGGTGTGGAAGCGTCATTtgcggct                                                                                                                         >  1:94667/1‑236 (MQ=255)
                                                                                                                    gATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGGAAATGTTCGGTGGGGTGGAGTTCTGGTTTGCGGTGATCAAAATCGTCGCCCTTGTGTCGCTGTTTGTCGTTGGCCTGGTCCTGTGGGCGGTGCGCTTTACGGCCCCGCCTGGGGGGGGAGCGGGCTTCGCCCCTTtgt                                                                                                                     >  2:87999/1‑239 (MQ=255)
                                                                                                                               aTGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACGTTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGAcgg                                                                                                           <  1:116877/238‑1 (MQ=255)
                                                                                                                                   cacgTTCGGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACggcggc                                                                                                       <  1:208606/235‑1 (MQ=255)
                                                                                                                                     ggTTCTTGTTGCCGGATTTCTCCGTCTGGGTCGCCGCGGTGGCGTTGTTAGTGCTGCTGCGGACGCTCAATCTCGCCACCGTGAAAATGGTCGGGGAGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTgg                                                                                                   <  2:510244/238‑1 (MQ=255)
                                                                                                                                     aGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGGGTTCTGGTTTTCGATGATCAAAATCGTCGCCATTTTGTAGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGCCTGGGGTGGAAGCGTCATTCGCGCATTTGTTGAATGACGGGGGCTgg                                                                                                   >  1:295346/1‑239 (MQ=255)
                                                                                                                                     aGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCCGACGCCCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGGGTGGAAGCGTCATTCGCGCATTTGTGGGATGACGGGGGGTgg                                                                                                   >  1:294690/1‑239 (MQ=255)
                                                                                                                                       ttCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTcgcg                                                                                                                            <  1:338395/213‑1 (MQ=255)
                                                                                                                                            ggtCCCCGTCCTCTCCGCATGGGTCGCCTGGTTGGCGGTTATAGTGCTGCTGCTGTCGCTCATGCCCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGTGGCTGGTTCCCGa                                                                                            <  2:473956/237‑1 (MQ=255)
                                                                                                                                            gTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGCTAGTGCTGCGGCCGGCGGCCAATTCCCCCACCGTGAAAAAGTTCGGGGAGATGGGGGTCTTGTTTGCGGTGAGCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATTGGTGCGGTGCTCCTTTAGTTACCGGCTAGTGTGGGACCGTCATTCCCGCCTTTGGGGACTGACGGCGGGCGGTTgccca                                                                                            >  2:719100/1‑237 (MQ=255)
                                                                                                                                               ccccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTTCTGACGCCCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGGTTGTCGGCGGCCTGGTCATGGTGGGCGTGGCCTTTCAGTCACCGCCCGGGTGGGGAAGCGCCTTCGCCCCCTTTGGGAATGAAGGGGGCGGGGTTCCGCAAg                                                                                         >  2:621880/1‑239 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGGGATGGGGTTCTGGTTTTCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGCCTGGTGTGGAAGCGTCATTCCCGCATTTTGTGAATGACGGCGGCTGGGTTCCCGacggg                                                                                        >  1:1070276/1‑235 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGGGATGGGGTTATGGTTTGCGCTGATCAAAATCCTCGCCATTGTGTCGCTGGTTTTCGTCGGCCTGGGCATGGTGGCGATGCCACTTCAGTCACCGACTGGTGGGGAGGCGCCATTTGCGCCTTTGGGGAATTAAGGCGGGCGGTTCCCCAGAGGt                                                                                       >  2:751407/1‑239 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGGGATGGAGTTCTGGTTTGCGATGCTCAAAATCGTCGCCATTGTGTCGCTGGTTTTCGTCGGCCTGGGCATGGTGGGGGGCGCCTTTCCTTCACCGGCTCGGGTGGGAAGGGATTTTGCGCCTTTGTGGAAAAACGGCGGGCGGTGCCCCAAAGGt                                                                                       >  2:76792/1‑239 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGGGGTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATTGTGGCGATGCCCTTTCCGTCACCCACTGGTGTGGAAGCGTCATTCGCGCATTTTTGGAAAGACGGGGGCGGGGTCCCGAAAggg                                                                                       >  2:917220/1‑238 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGt                                                                                       >  1:502069/1‑239 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGGTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGCGCTGGGCATGGTGGCGATGCACTTTTACGCACCCACTGGGGGTGAAGCGTCCTTTCTGCCTGTTTGGAATGACGGCGGCTGGGTCCCGaaaagg                                                                                       >  2:181175/1‑236 (MQ=255)
                                                                                                                                                          ccGCCTGGGTCCCCTCGCTGGGGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGc                                                                              <  2:687327/239‑1 (MQ=255)
                                                                                                                                                                 gggccccTCGCTGGGGGTGATAGTGCTGGTGCTGACGCGCAATCTCGCCACCGTGAAAATTGTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTg                                                                       <  2:56209/234‑1 (MQ=255)
                                                                                                                                                                  gTCGCCTGGCTGGGGGTTATTGTGCTTCTGCTTGCGCTCAATCTCGCCACCGTGAAAATGGTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGGCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTg                                                                       <  1:21939/238‑1 (MQ=255)
                                                                                                                                                                    ccgcTCGCTTGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGcc                                                                     <  1:917220/235‑1 (MQ=255)
                                                                                                                                                                                      aGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGGCTGGTGTGGAAGCGCCATTCGCGCATTTTCGGAATGACGGCGGCGGGTTCCCGGAAGGTTTAAAGGGGTTCTTTTTCGGATCCCAGAAAGCGGgtt                                                  >  2:468806/1‑239 (MQ=255)
                                                                                                                                                                                       gtggtgCTGGTGGCGCTCAATCCTCCCACCGTGAAAATTTTCGGTGGGGTGGATTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGtttt                                                 <  2:468767/239‑1 (MQ=255)
                                                                                                                                                                                       gTGCTGCTGCTGACGCTTAATCCCGCCACCGTGAAAATTTTCGGTGGGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGtttt                                                 <  2:502069/239‑1 (MQ=255)
                                                                                                                                                                                        tgctgctgcTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGtttt                                                 >  1:373999/1‑238 (MQ=255)
                                                                                                                                                                                        tgctgctgcTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGAGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCCTTTGTGGAATGACGGCGGCTGGTTCCCGGAAGGTTTAAGTGGCTTCTTTGCCGGATTCCTGATAGCGGgttt                                                 >  1:603649/1‑238 (MQ=255)
                                                                                                                                                                                           tgctgctgACGCTCAATCGCGCCAACGTGAAAATGTTCGGTGAGATGGAGGTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCt                                             <  1:112619/239‑1 (MQ=255)
                                                                                                                                                                                            gctgctgACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCTTTCGCGCATTTGGGGAATGACGGGGGCTGGTTCCCGAAAGGTTTAAGGGGGTTCTTTTGCGGATTCAAGATAGCGGGTTTCGGtt                                            >  2:109245/1‑239 (MQ=255)
                                                                                                                                                                                            gctgctgACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGGGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCTTTCGCGCATTTGTGGAATGACGGCGGCCGGTTTCCGGCAGGTTTAAGGGGTCTCTTTGCCGGGTTCCCATTAACGGttttgcctt                                            >  2:325567/1‑234 (MQ=255)
                                                                                                                                                                                                 tgACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGt                                        <  1:588868/238‑1 (MQ=255)
                                                                                                                                                                                                     ggccAAGCTCGCCCACGTGAAAATGTTCGGGGAGATGGATTTTTGGTTTGCGATGATCAAAATCGGCGCCATTGTTGCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCCTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGa                                   <  2:603649/236‑1 (MQ=255)
                                                                                                                                                                                                     gccaaaTCGCCCCCCCGTGAAAATGTTCGGGGGGATGGGGTTTTGGTTTGGGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGa                                   <  2:94667/235‑1 (MQ=255)
                                                                                                                                                                                                     gCTCAATCTCGCCACCGTGAAAATTTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGa                                   <  1:381298/239‑1 (MQ=255)
                                                                                                                                                                                                          aTCTCGCCTACGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGGCAGTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAg                              <  1:55107/239‑1 (MQ=255)
                                                                                                                                                                                                             tctccacgttGAAAATTTTCGGTGAATTGGTGTTCTGGTTTGCGGTGATCAAAATCGGCGCCTTTGGGTCGCTGATTGTCGTCGGTCTGGTCATGGTGGCGATGCCCTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTg                           <  2:183863/230‑1 (MQ=255)
                                                                                                                                                                                                             tccccACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGAACAAAAACGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTg                           <  1:633466/236‑1 (MQ=255)
                                                                                                                                                                                                               gCCCCCGTAAAAATGGTCGGTTAGATGGAGTTCTGGGTTGCGATGATCAAAATCGTCGCCATTGGGTCGCTGATTGTCGTCGGCCGGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt                         <  1:156424/239‑1 (MQ=255)
                                                                                                                                                                                                               gCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAACGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAACGGTTTTCGCTTTCGTGGGGATTGGCCTgga                         >  1:497491/1‑238 (MQ=255)
                                                                                                                                                                                                               gCCACCGGTAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGAACAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt                         <  1:791993/239‑1 (MQ=255)
                                                                                                                                                                                                                tcctcGTTAAATTTTTCGGTTAGTTGGATTTTTGTTTTGCTAGGAACAAAGTCGTCGTCTTTGTGTGCCTGATTGGTGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGGCATGTGCGTATTTGTGGAAGGGAGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTCCTTTGCCGGATTCCAGATAGCGGTGTTCGCTTTCGTGGGGATTGAGCTggga                        <  2:55143/235‑3 (MQ=255)
                                                                                                                                                                                                                   ccGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGGTTAAGGGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGGGGGGGTTGAGGTgggaggg                     >  2:791993/1‑238 (MQ=255)
                                                                                                                                                                                                                                tctgTGATATGGAGTTCTGGTTTGCGATGATCAAAAGCGTCGGCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCTCTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCg         <  1:712609/236‑1 (MQ=255)
                                                                                                                                                                                                                                 cGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGCTAGCGGTTTTCGCTTTCTGGGGGATTGAGCTGGTAGGTACACCAGATTCgg        >  1:490158/1‑239 (MQ=255)
                                                                                                                                                                                                                                  ggTGAGGTGGATTTCTGGTTTGCGATGATCAAACTCGGCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGa       <  1:640199/239‑1 (MQ=255)
                                                                                                                                                                                                                                   gTGAGTTGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGGCGCTGATTGGCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGaa      <  2:1055163/239‑1 (MQ=255)
                                                                                                                                                                                                                                        gtGGAGTTCTGGTTTGCGGGTAACAAAATCGTCGCTATTGTGTCGCGTATTGTCGGCGGCCGGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGGATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCa  <  2:1070276/237‑1 (MQ=255)
                                                                                                                                                                                                                                        aTGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcg                                                                                                                                                                                           <  2:561491/55‑1 (MQ=255)
                                                                                                                                                                                                                                        aTGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcg                                                                                                                                                                                           >  1:561491/1‑55 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
GAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCA  >  NC_000913/4430110‑4430578

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: