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

GAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCA  >  NC_000913/4430110‑4430578
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
gattACAAATTTTTTAGGGACTTCGCTTCCGATTTTCTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGtt                                                                                                                                                                                                                                          <  2:786454/236‑1 (MQ=255)
          cTTTTAGTGACTTCGTTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCg                                                                                                                                                                                                                                <  1:49592/239‑1 (MQ=255)
             ttAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGCGGTAGTGCTGCTGCTGACGCTCAATCTCGCCCCCCTGCAAAATG‑TCGGTGAGATGGAGATCTGGTTTGCgctt                                                                                                                                                                                                                             >  1:171723/1‑236 (MQ=255)
                        gtttcggtttttCTGGGGGCGTGGGCAGGATATTTCCCCGGCTGGGCTTACTTGTTCTGCTGGGTTGTAACGGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTATGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAtcgt                                                                                                                                                                                                                  <  2:577549/231‑1 (MQ=255)
                        gcctctgATTTGCTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGGAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAtcgt                                                                                                                                                                                                                  <  1:1035868/233‑1 (MQ=255)
                        gCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAtcgt                                                                                                                                                                                                                  <  1:77247/239‑1 (MQ=255)
                                 ttACTCGGGCCGTGGGCAGGATATTTCACCGGCTTGACTTACTGGTTCTGCTTGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATtgt                                                                                                                                                                                                         <  1:773782/239‑1 (MQ=255)
                                                 cAGGATATTTCCCCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcgt                                                                                                                                                                                          <  1:100527/238‑1 (MQ=255)
                                                 cAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcgt                                                                                                                                                                                          <  1:399947/238‑1 (MQ=255)
                                                         ttCCCCGGCTGGACTTACTGGGTCTGCTGGGTTGTAACCGGTATGGCAGGCGGTGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGt                                                                                                                                                                                 <  1:610035/239‑1 (MQ=255)
                                                              cGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTTGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGTTAGTTCTGCTGCTGACGCTCCATCTCGCCCCCCTG‑GAATTTTTTGGGGGGGTGGGGTTCTGGTTTGCGATGATCACAACGGTAGCCCTTTTGTCGCTGATTGTTGTTggagtgggacagg                                                                                                                                                                            >  2:55107/1‑228 (MQ=255)
                                                              cGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCCTTGTGTCGCTGATTGTCGTGGGCCTGGTCAtgg                                                                                                                                                                            >  1:56209/1‑239 (MQ=255)
                                                                               tCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGCTCAAAATCGTCGCCATTGTGTCGCTGATTTTCGGCGGCCTGTTCCTGGTGGCGATGCACTTTCa                                                                                                                                                            >  1:55143/1‑238 (MQ=255)
                                                                                ctgctgGGTTGTAACCGGTTGGGCAGAAGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAg                                                                                                                                                           <  1:737701/238‑1 (MQ=255)
                                                                                   ctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCCCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTAGGTTGCGCTGATCAAAAACGTCGCCATTGTGTCGCTGATTTTCGTCCGCCTGGTCCTGGTGGCTGTGCAATTTCAGTCa                                                                                                                                                        >  1:766481/1‑238 (MQ=255)
                                                                                   ctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGGGCTGCTGCTGACGCTCAATCTCGCCCCCGCG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCCTTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGGTGCCATTTCAggaa                                                                                                                                                        >  1:183863/1‑235 (MQ=255)
                                                                                   acaggTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCa                                                                                                                                                        <  1:64296/235‑1 (MQ=255)
                                                                                       gTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAACGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGGTGGACTTTCAGTCACCG‑Ac                                                                                                                                                  >  1:687327/1‑239 (MQ=255)
                                                                                           tAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGATGCTGACGCCCAATCTCGCCCCCGTG‑AAAATGTTCGGTGAGTTGGGGGTCTGGGTTGCGATGATCAAAATGGTCGCCCTTGTTGCGCGGTCTTTGGCAGGCCGGGGTATGGGGGGGGAGCACTTTCAGTCACCG‑ACgggg                                                                                                                                              >  2:105589/1‑238 (MQ=255)
                                                                                                ggTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGCCTGGGTCGCCTCGCTGGCGGTGAGAGTGCTGCTGCTGCCCCCCAATCTCGCCACCGCG‑GAAATCTTCGGGGTGGCGGGGTTTTGGTTTGCCGTTATCAAAATCGTCCCCCTTGTTAAGGCGGTTGTCGTTGTCCTGGGCATGGTGGGAGGGGACTTTTCGTCCCCG‑ACTgggggggg                                                                                                                                         >  2:419420/1‑233 (MQ=255)
                                                                                                  tATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGGG‑AAAATGTTCGGTGAGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTTGTCGGCCTGGTCATGGGGGCGGTGCACTTTTAGTCACCG‑ACTggggggggag                                                                                                                                       >  1:468767/1‑231 (MQ=255)
                                                                                                            gtggtgGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTTTCGTCGGCCTGGTCATGGTGGGGATGCACTTTCATTCACCG‑ACTGGTGTGGAAGGGGCATTggc                                                                                                                             >  1:525683/1‑239 (MQ=255)
                                                                                                            gtggtgGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTTTCGTCGGCCTTGGCATGGTGGCGATGCACTTTCAGTCACCG‑CCTGGGGTGGAAGCGTCATTccc                                                                                                                             >  1:313112/1‑237 (MQ=255)
                                                                                                              ggtggCGATCCCGGCTTTTGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTcgcg                                                                                                                            <  1:1083437/238‑1 (MQ=255)
                                                                                                                tggCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACCGGTGTGGAAGCGTCATTtgcggct                                                                                                                         >  1:94667/1‑236 (MQ=255)
                                                                                                                    gATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑GAAATGTTCGGTGGGGTGGAGTTCTGGTTTGCGGTGATCAAAATCGTCGCCCTTGTGTCGCTGTTTGTCGTTGGCCTGGTCCTGTGGGCGGTGCGCTTTACGGCCCCG‑CCTGGGGGGGGAGCGGGCTTCGCCCCTTtgt                                                                                                                     >  2:87999/1‑239 (MQ=255)
                                                                                                                               aTGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACGTTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGAcgg                                                                                                           <  1:116877/238‑1 (MQ=255)
                                                                                                                                   cacgTTCGGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACggcggc                                                                                                       <  1:208606/235‑1 (MQ=255)
                                                                                                                                     ggTTCTTGTTGCCGGATTTCTCCGTCTGGGTCGCCGCGGTGGCGTTGTTAGTGCTGCTGCGGACGCTCAATCTCGCCACCGTG‑AAAATGGTCGGGGAGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑gg                                                                                                   <  2:510244/238‑1 (MQ=255)
                                                                                                                                     aGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGGGTTCTGGTTTTCGATGATCAAAATCGTCGCCATTTTGTAGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑CCTGGGGTGGAAGCGTCATTCGCGCATTTGTTGAATGACGGGGGCT‑gg                                                                                                   >  1:295346/1‑239 (MQ=255)
                                                                                                                                     aGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCCGACGCCCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGGGTGGAAGCGTCATTCGCGCATTTGTGGGATGACGGGGGGT‑gg                                                                                                   >  1:294690/1‑239 (MQ=255)
                                                                                                                                       ttCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTcgcg                                                                                                                            <  1:338395/213‑1 (MQ=255)
                                                                                                                                            ggtCCCCGTCCTCTCCGCATGGGTCGCCTGGTTGGCGGTTATAGTGCTGCTGCTGTCGCTCATGCCCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGTGGCT‑GGTTCCCGa                                                                                            <  2:473956/237‑1 (MQ=255)
                                                                                                                                            gTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGCTAGTGCTGCGGCCGGCGGCCAATTCCCCCACCGTG‑AAAAAGTTCGGGGAGATGGGGGTCTTGTTTGCGGTGAGCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATTGGTGCGGTGCTCCTTTAGTTACCG‑GCTAGTGTGGGACCGTCATTCCCGCCTTTGGGGACTGACGGCGGGC‑GGTTgccca                                                                                            >  2:719100/1‑237 (MQ=255)
                                                                                                                                               ccccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTTCTGACGCCCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGGTTGTCGGCGGCCTGGTCATGGTGGGCGTGGCCTTTCAGTCACCGCCCGGGTGGGGAAGCG‑CCTTCGCCCCCTTTGGGAATGAAGGGGGCG‑GGGTTCCGCAAg                                                                                         >  2:621880/1‑239 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGGGATGGGGTTCTGGTTTTCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑CCTGGTGTGGAAGCGTCATTCCCGCATTTTGTGAATGACGGCGGCTGGGTTCCCGacggg                                                                                        >  1:1070276/1‑235 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGGGATGGGGTTATGGTTTGCGCTGATCAAAATCCTCGCCATTGTGTCGCTGGTTTTCGTCGGCCTGGGCATGGTGGCGATGCCACTTCAGTCACCG‑ACTGGTGGGGAGGCGCCATTTGCGCCTTTGGGGAATTAAGGCGGGC‑GGTTCCCCAGAGGt                                                                                       >  2:751407/1‑239 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGGGATGGAGTTCTGGTTTGCGATGCTCAAAATCGTCGCCATTGTGTCGCTGGTTTTCGTCGGCCTGGGCATGGTGGGGGGCGCCTTTCCTTCACCG‑GCTCGGGTGGGAAGGGATTTTGCGCCTTTGTGGAAAAACGGCGGGC‑GGTGCCCCAAAGGt                                                                                       >  2:76792/1‑239 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGGGGTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATTGTGGCGATGCCCTTTCCGTCACCC‑ACTGGTGTGGAAGCGTCATTCGCGCATTTTTGGAAAGACGGGGGCG‑GGGTCCCGAAAggg                                                                                       >  2:917220/1‑238 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGt                                                                                       >  1:502069/1‑239 (MQ=255)
                                                                                                                                                 ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGGTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGCGCTGGGCATGGTGGCGATGCACTTTTACGCACCC‑ACTGGGGGTGAAGCGTCCTTTCTGCCTGTTTGGAATGACGGCGGCT‑GGGTCCCGaaaagg                                                                                       >  2:181175/1‑236 (MQ=255)
                                                                                                                                                          ccGCCTGGGTCCCCTCGCTGGGGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGc                                                                              <  2:687327/239‑1 (MQ=255)
                                                                                                                                                                 gggccccTCGCTGGGGGTGATAGTGCTGGTGCTGACGCGCAATCTCGCCACCGTG‑AAAATTGTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTg                                                                       <  2:56209/234‑1 (MQ=255)
                                                                                                                                                                  gTCGCCTGGCTGGGGGTTATTGTGCTTCTGCTTGCGCTCAATCTCGCCACCGTG‑AAAATGGTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGGCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTg                                                                       <  1:21939/238‑1 (MQ=255)
                                                                                                                                                                    ccgcTCGCTTGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGcc                                                                     <  1:917220/235‑1 (MQ=255)
                                                                                                                                                                                      aGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑GCTGGTGTGGAAGCGCCATTCGCGCATTTTCGGAATGACGGCGGCG‑GGTTCCCGGAAGGTTTAAAGGGGTTCTTTTTCGGATCCCAGAAAGCGGgtt                                                  >  2:468806/1‑239 (MQ=255)
                                                                                                                                                                                       gtggtgCTGGTGGCGCTCAATCCTCCCACCGTG‑AAAATTTTCGGTGGGGTGGATTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGtttt                                                 <  2:468767/239‑1 (MQ=255)
                                                                                                                                                                                       gTGCTGCTGCTGACGCTTAATCCCGCCACCGTG‑AAAATTTTCGGTGGGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGtttt                                                 <  2:502069/239‑1 (MQ=255)
                                                                                                                                                                                        tgctgctgcTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGtttt                                                 >  1:373999/1‑238 (MQ=255)
                                                                                                                                                                                        tgctgctgcTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGAGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCCTTTGTGGAATGACGGCGGCT‑GGTTCCCGGAAGGTTTAAGTGGCTTCTTTGCCGGATTCCTGATAGCGGgttt                                                 >  1:603649/1‑238 (MQ=255)
                                                                                                                                                                                           tgctgctgACGCTCAATCGCGCCAACGTG‑AAAATGTTCGGTGAGATGGAGGTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCt                                             <  1:112619/239‑1 (MQ=255)
                                                                                                                                                                                            gctgctgACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCTTTCGCGCATTTGGGGAATGACGGGGGCT‑GGTTCCCGAAAGGTTTAAGGGGGTTCTTTTGCGGATTCAAGATAGCGGGTTTCGGtt                                            >  2:109245/1‑239 (MQ=255)
                                                                                                                                                                                            gctgctgACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGGGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCTTTCGCGCATTTGTGGAATGACGGCGGCC‑GGTTTCCGGCAGGTTTAAGGGGTCTCTTTGCCGGGTTCCCATTAACGGttttgcctt                                            >  2:325567/1‑234 (MQ=255)
                                                                                                                                                                                                 tgACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGt                                        <  1:588868/238‑1 (MQ=255)
                                                                                                                                                                                                     ggccAAGCTCGCCCACGTG‑AAAATGTTCGGGGAGATGGATTTTTGGTTTGCGATGATCAAAATCGGCGCCATTGTTGCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCCTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGa                                   <  2:603649/236‑1 (MQ=255)
                                                                                                                                                                                                     gccaaaTCGCCCCCCCGTG‑AAAATGTTCGGGGGGATGGGGTTTTGGTTTGGGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGa                                   <  2:94667/235‑1 (MQ=255)
                                                                                                                                                                                                     gCTCAATCTCGCCACCGTG‑AAAATTTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGa                                   <  1:381298/239‑1 (MQ=255)
                                                                                                                                                                                                          aTCTCGCCTACGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGGCAGTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAg                              <  1:55107/239‑1 (MQ=255)
                                                                                                                                                                                                             tctccacgttG‑AAAATTTTCGGTGAATTGGTGTTCTGGTTTGCGGTGATCAAAATCGGCGCCTTTGGGTCGCTGATTGTCGTCGGTCTGGTCATGGTGGCGATGCCCTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTg                           <  2:183863/230‑1 (MQ=255)
                                                                                                                                                                                                             tccccACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGAACAAAAACGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTg                           <  1:633466/236‑1 (MQ=255)
                                                                                                                                                                                                               gCCCCCGTA‑AAAATGGTCGGTTAGATGGAGTTCTGGGTTGCGATGATCAAAATCGTCGCCATTGGGTCGCTGATTGTCGTCGGCCGGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt                         <  1:156424/239‑1 (MQ=255)
                                                                                                                                                                                                               gCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAACGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAACGGTTTTCGCTTTCGTGGGGATTGGCCTgga                         >  1:497491/1‑238 (MQ=255)
                                                                                                                                                                                                               gCCACCGGT‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGAACAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt                         <  1:791993/239‑1 (MQ=255)
                                                                                                                                                                                                                tcctcGTT‑AAATTTTTCGGTTAGTTGGATTTTTGTTTTGCTAGGAACAAAGTCGTCGTCTTTGTGTGCCTGATTGGTGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGGCATGTGCGTATTTGTGGAAGGGAGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTCCTTTGCCGGATTCCAGATAGCGGTGTTCGCTTTCGTGGGGATTGAGCTggga                        <  2:55143/235‑3 (MQ=255)
                                                                                                                                                                                                                   ccGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGGTTAAGGGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGGGGGGGTTGAGGTgggaggg                     >  2:791993/1‑238 (MQ=255)
                                                                                                                                                                                                                                tctgTGATATGGAGTTCTGGTTTGCGATGATCAAAAGCGTCGGCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCTCTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCg         <  1:712609/236‑1 (MQ=255)
                                                                                                                                                                                                                                 cGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGCTAGCGGTTTTCGCTTTCTGGGGGATTGAGCTGGTAGGTACACCAGATTCgg        >  1:490158/1‑239 (MQ=255)
                                                                                                                                                                                                                                  ggTGAGGTGGATTTCTGGTTTGCGATGATCAAACTCGGCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGa       <  1:640199/239‑1 (MQ=255)
                                                                                                                                                                                                                                   gTGAGTTGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGGCGCTGATTGGCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGaa      <  2:1055163/239‑1 (MQ=255)
                                                                                                                                                                                                                                        gtGGAGTTCTGGTTTGCGGGTAACAAAATCGTCGCTATTGTGTCGCGTATTGTCGGCGGCCGGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGGATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCa  <  2:1070276/237‑1 (MQ=255)
                                                                                                                                                                                                                                        aTGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcg                                                                                                                                                                                           <  2:561491/55‑1 (MQ=255)
                                                                                                                                                                                                                                        aTGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcg                                                                                                                                                                                           >  1:561491/1‑55 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
GAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTG‑AAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCG‑ACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCT‑GGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCA  >  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: ‑