Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,539,0510AG70.5% 113.9 / 14.4 63N3N (AAT→AAC) unknownCyclic pyranopterin monophosphate synthase (EC 4.6.1.17)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (36/7);  minor base T (0/18);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 3.66e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAATTCCCACGTGCGTTTTGCCGCCTGAATGCCGGCAATGCGGGCGGTGGCGAACACATCGCCTTTGTGATGCTGACCGGACAGGATCATTTGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTATTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCGCCGAATAAACAAAGGTGCAATTTCCCTTTGGCGGAAACCCGTAAACGGTTACAG  >  1/1538805‑1539285
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             
cAATTCCCACGTGCGTTTTGCCGCCTGAATGCCGGCAATGCGGGCGGTGGCGAACACATCGCCTTTGTGATGCTGACCGGACAGGATCATTTGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGtttt                                                                                                                                                                                                                                           >  1:750914/1‑251 (MQ=255)
  aTTCCCACGTGCGTTTTGCCGCCTGAATGCCGGCAATGCGGGCGGTGGCGAACACATCGCCTTTGTGATGCTGACCGGACAGGATCATTTGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGttttt                                                                                                                                                                                                                                          >  1:62263/1‑250 (MQ=255)
                 ttGCCGCCTGAATGCCGGCAATGCGGGCGGTGGCGAACACATCGCCTTTGTGATGCTGACCGGACAGGATCATTTGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCAt                                                                                                                                                                                                                       <  1:557684/251‑1 (MQ=255)
                     cgccTGAATGCCGGCAATGCGGGCGGTGGCGAACACATCGCCTTTGTGATGCTGACCGGACAGGATCATTTGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAAc                                                                                                                                                                                                                   <  1:69813/251‑1 (MQ=255)
                       ccTGAATGCCGGCAATGCGGGCGGTGGCGAACACATCGCCTTTGTGATGCTGACCGGACAGGATCATTTGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAAcc                                                                                                                                                                                                                  <  1:376771/250‑1 (MQ=255)
                                       gcgggcggTGGCGAACACATCGCCTTTGTGATGCTGACCGGACAGGATCATTTGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAAc                                                                                                                                                                                                                   <  1:372039/233‑1 (MQ=255)
                                             ggTGGCAACCACTTCGCTTTTGTGATGCTGCCCGGACAGGATCTTTTGTAAGGTTTCGGGCGCCAAACAAACAAAAGCCTCGGCGCGGGCTTCACGGACGGTTTCGGTTTT‑GGCGGAAACGCCCACCATATTGGCTCCCCC‑ATGGA‑AA‑TTAAT‑GTGAGAAAATGTGACCATGATAATAAAAGTGGGGTTGT‑TTTTCTTGTTTTTTTTA‑AAAAT‑C‑GGGGTGCATTAACCGCCGATACTTGCCAAATGACta                                                                                                                                                                                            <  1:281306/250‑1 (MQ=255)
                                                 gCGAACACTCCGCTTTTGTGATGCTGACCGGACAGGATCTTTTGTAAGTTTTGGGGCGCCATAAAAACATAACCCTCAGCGCGGGCTTCACGGAGGTTTTCGGTTTG‑GGCGGAAACGTCCACCATTTTGGCTTCCCA‑ATGGG‑AA‑TTATT‑GTGATTAAATGTGACCATAATAATAAAGGTGCGTTTTT‑TTTTTTTGTTTTTTTTA‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACtacgta                                                                                                                                                                                        <  1:487616/250‑1 (MQ=255)
                                                            gCCTTTGTGATGCTGACCGGACAGGATCATTTGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCg                                                                                                                                                                             >  1:763609/1‑250 (MQ=255)
                                                                               gACAGGATTTTTTGTAAGGTTGCGGGCGCAATAAAAAAATAAGCCTCGGCGGGGGCTTCAGGGAGGGTTTCGGTTTT‑GGCGAAAACGCCCACCATATGGGCTCCCCC‑AGGGG‑AA‑TTAAT‑GTGAGTAAATGTGCCCATGATAATAAAAGTGCGGTTGA‑TTTTTTTGCTGTTTTAA‑AAAAT‑C‑GGCGTGCATAAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAatgat                                                                                                                                                          <  1:592206/250‑1 (MQ=255)
                                                                                   ggATCTTTTGAAAGGTTGCGGGCGCAATAAAAACAAAAGCCTCGGCGGGGGCTTCACGGACGTTTTCGGTTTT‑GGCGGAAACGTCCACCATTTTGGCTCCCCC‑ATGGA‑AA‑TTAAT‑GTGAGTAAATGTGCCCATAATAATTAAAGTGCGTTTGA‑TTTTCTTGCTTTTTTTA‑AAAAC‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTc                                                                                                                                                      <  1:68246/250‑1 (MQ=255)
                                                                                   ggATCATTTGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCa                                                                                                                                                     >  1:491158/1‑251 (MQ=255)
                                                                                    gATCATTTGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATg                                                                                                                                                                                                >  1:887906/1‑207 (MQ=255)
                                                                                       cATTTGTAGGGTTTGGGGCCCCATAAAAACATAAGCCTCGGCGGGGGCTTCACGGAGGGTTGCGGTTTT‑GGGGGAACCGTCTACCATTTTGGCTCCCCC‑ATGGA‑AA‑TTAAT‑GTGAGAAAATGTAACCATGAAAATAAAAGTGGGTTTGT‑TTTTTTTGTTTTTTTAA‑AAAAT‑C‑GGGGTGAATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCt                                                                                                                                                 <  1:438616/251‑1 (MQ=255)
                                                                                         tttGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTa                                                                                                                                               >  1:194584/1‑251 (MQ=255)
                                                                                             tagggTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGTTTT‑GGCGGAAACCTCCACAATATTGGCTCCCCC‑AGGGA‑AA‑TTAAT‑GTGAGAAAATGTGACCATAATAATAAAAGTGCGGTTGT‑TTTTCTTGTTTTTTTTA‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCt                                                                                                                                            <  1:794017/247‑1 (MQ=255)
                                                                                             aaaGGTTTGGGGCACTATAAAACCAAAAGCCTCAGCGCGGGCTACGCGAAGGGTTTGGGTTTG‑GGCGGAACCGTCCACGATTTTGG‑TTCCCCTAGGGA‑AA‑TTATG‑GTGATTAATTGTGACCATGATAATAAAAGTG‑GGTTTATTTTTCTGGCTTTTTTTA‑AAAAT‑G‑GGGGTGTATAAGCCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCt                                                                                                                                            <  1:327155/249‑1 (MQ=255)
                                                                                                             aTACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGcctgtttctt                                                                                                                           >  1:772694/1‑242 (MQ=255)
                                                                                                             aTACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAAt                                                                                                                           >  1:430650/1‑251 (MQ=255)
                                                                                                                  aaCATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTg                                                                                                                      >  1:312210/1‑251 (MQ=255)
                                                                                                                   aCATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTACAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCt                                                                                                                                            >  1:58484/1‑228 (MQ=255)
                                                                                                                   aCATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGa                                                                                                                     >  1:375268/1‑251 (MQ=255)
                                                                                                                    cATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGc                                                                                                                                                                    >  1:94239/1‑203 (MQ=255)
                                                                                                                    cATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAat                                                                                                                                                             >  1:97899/1‑210 (MQ=255)
                                                                                                                         gCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCt                                                                                                                                            >  1:718695/1‑222 (MQ=255)
                                                                                                                             cAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATg                                                                                                           >  1:400091/1‑251 (MQ=255)
                                                                                                                               gcgggggCTTCAGGGACGGTTTCGGTTTG‑GGCGAAAACGCCCACCATATGGGCTCCCCC‑ATGGG‑AA‑TAAAT‑GTGAGTAATTGTGACAATGATAATAAAGGTGCGGTTGT‑TTTTCTTGCTTTTTTTA‑AAAAT‑C‑GGCGGGCTTTAACCGCCGATACTTGCCAAATGACTACGTACGCCAc                                                                                                                                                                                  <  1:873162/174‑1 (MQ=255)
                                                                                                                                        tCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGt                                                                                                                                                        >  1:806019/1‑195 (MQ=255)
                                                                                                                                              aggttttCGGTTTT‑GGCGGAAACGTCCACAATATTGGCTCCCCC‑AGGGA‑AA‑TTAAT‑GTGAGTAAATGTGCCAATGATAATAAAAGTGCGGTTGT‑TTTTTTTGCTTTTTTTA‑AAAAT‑C‑GGCGTGAATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGt                                                                                           <  1:114468/246‑1 (MQ=255)
                                                                                                                                              acgttttCGGTTTT‑GGCGAAAACGCCCACCATATTGGCTTCCCC‑ATGGA‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGTTGTTTTTC‑AAAAT‑C‑GGCGTGCATAAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGt                                                                                           <  1:588873/246‑1 (MQ=255)
                                                                                                                                              acggTTTCGGTTTT‑GGCGGAACCGCCCACCATATTGGCTCCCCC‑ATGGA‑AA‑TTAAT‑GTGAGAAAATGTGACAATGATAATAAAAGTGCGTTTGT‑TTTTCTTGTTTTTTTTA‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGAAGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGt                                                                                           <  1:15591/250‑1 (MQ=255)
                                                                                                                                               cggTTTCGGGTTTGGGAGGAACCTTCCACCATATGGGC‑TCCCC‑AAGGG‑GA‑TTTATGGTGAGTAAATGTGCCAAGGATAATTAAAG‑GGGGTTTT‑TTTTTTTGGTTTTTTTA‑AAAATGC‑GGGTTGAATAAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGt                                                                                           <  1:363804/250‑1 (MQ=255)
                                                                                                                                                    tCGGCTTG‑GGGGGAACCGTCCACCATATTGGTTCCCCC‑AGGGG‑AT‑TTAAT‑GTGATTAATGGGGACCATGATAATGAAAGTGGGGTTGA‑TTTT‑TTTCTTGTTTTA‑AAAAA‑CGGGGGTGCATAAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTaa                                                                                    <  1:261139/251‑1 (MQ=255)
                                                                                                                                                     cgttttt‑GGCGGAACCTTCCACCTTTTTGCCTTCCCC‑A‑GGG‑AATTTAAT‑GTGATTAATTTTGACCATTATATTTAAGGTGGGTTTTT‑TTTTTTTGCT‑TTTTTCAAAAAT‑C‑GGCGTGAATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAAc                                                                                   <  1:466331/247‑1 (MQ=255)
                                                                                                                                                        gtttgggcggaacctcccaccaTTTTGCC‑TCCCC‑AAGGGAAA‑TAAAT‑GTAAGTAAATGTGCCCATGATAATTAAAGTGCGGTTGT‑TTTTCTTGTTTTTTTTA‑AAAGT‑G‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACa                                                                                  <  1:850667/234‑1 (MQ=255)
                                                                                                                                                          ttggggggaacctcccacaaTATTGGCTTCCCG‑TGGGG‑AA‑TTAAT‑GTGAGAAAATGTGCCAATGAAAATAAAAGTGGGTTTGT‑TTTTTTTGTTGTTTATA‑AAAAC‑G‑GGCGTACTTTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAAt                                                                               <  1:399756/237‑1 (MQ=255)
                                                                                                                                                             ggagaaaccgcccaccaTATTGGCTCCCCG‑AGGGG‑AA‑TTATT‑GGGAATAAGTGTGACAGTGAGAATTAAGGTG‑GGTTGATTTTTCTTGCT‑TTTTTA‑AAAATGG‑GGGGTGAATAAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGAt                                                                                                                                                                <  1:39227/156‑1 (MQ=255)
                                                                                                                                                                     cGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCa                                                                    >  1:731511/1‑251 (MQ=255)
                                                                                                                                                                      gTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTc                                                                     >  1:387996/1‑249 (MQ=255)
                                                                                                                                                                        ccaccaTATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTaa                                                                                    >  1:165674/1‑232 (MQ=255)
                                                                                                                                                                           ccaTATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAAt                                                                               >  1:299711/1‑234 (MQ=255)
                                                                                                                                                                           ccaTATTGCCTCCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATAATAATAAAAGTGCGTTTGT‑TTTTCTTGTTTTTTTTA‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGt                                                                                      <  1:201005/227‑1 (MQ=255)
                                                                                                                                                                             atatTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAAc                                                                                   >  1:527808/1‑228 (MQ=255)
                                                                                                                                                                                  ggCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCg                                                        >  1:733175/1‑250 (MQ=255)
                                                                                                                                                                                         cc‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAAt                                                                               >  1:55194/1‑220 (MQ=255)
                                                                                                                                                                                               ggaaa‑TTTAT‑GTGATTAAATGTGCCCATTAAAATTAAGGTGCGTTTTT‑TTTTTTTCTTTTTTTTC‑AAAAT‑C‑GGCGTACATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCGCCGAATAAACa                                             <  1:413606/247‑1 (MQ=255)
                                                                                                                                                                                                g‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGcc                                                                                                                                                                           >  1:743942/1‑122 (MQ=255)
                                                                                                                                                                                                g‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCaaaa                                                                                                             >  1:429311/1‑184 (MQ=255)
                                                                                                                                                                                                g‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCGCCGAATAAACAAAg                                          >  1:758091/1‑251 (MQ=255)
                                                                                                                                                                                                 aaa‑TTATT‑GTGAGTAAAGGTGACCATGAAAATAAAGGGGCGTTTGT‑TTTTCTTG‑TTTTTTTA‑AAAATGT‑GGCGGATATAAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCGCCGAATAAACaaa                                           <  1:258517/249‑1 (MQ=255)
                                                                                                                                                                                                  aa‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCt                                                                                                                                            >  1:615531/1‑152 (MQ=255)
                                                                                                                                                                                                        tt‑gtgATTAATTGTGACCATGATAATAAAAGTGCGTTTGT‑TTTTTTTGTTTTTTTTA‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAAt                                                                               <  1:362336/207‑1 (MQ=255)
                                                                                                                                                                                                        aT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCGCCGAATAAACAAAGGTGCaa                                    >  1:511733/1‑251 (MQ=255)
                                                                                                                                                                                                             gAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAAt                                                                       >  1:182124/1‑212 (MQ=255)
                                                                                                                                                                                                             gAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGc                                                                >  1:372700/1‑219 (MQ=255)
                                                                                                                                                                                                             gAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCGCCGAATAAACAAAGGTGCAAttt                                 >  1:749994/1‑250 (MQ=255)
                                                                                                                                                                                                             gAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCGCCGAATAAACAAAGGTGCAAttt                                 >  1:27230/1‑250 (MQ=255)
                                                                                                                                                                                                             aataaaaTGTGACCATTATAATAAAGGTGCGGTTGT‑TTTTTTTCCTTTTTTAC‑AAAAC‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCGCCGAATAAACAAAGGTGCAATTTc                                <  1:94136/247‑1 (MQ=255)
                                                                                                                                                                                                                               taataaAAGGGCGGTTTT‑TTTTTTTGTTTTTTTTA‑AAAAT‑C‑GGCGTGCATAAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAatgat                                                                                                                                                          <  1:319922/111‑1 (MQ=255)
                                                                                                                                                                                                                               tAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAatga                                                                                                                                                           >  1:705622/1‑110 (MQ=255)
                                                                                                                                                                                                                               tAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTAt                                                                                         >  1:486251/1‑176 (MQ=255)
                                                                                                                                                                                                                               tAATTAAAGTGCGGTTGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCGCCGAATAAACAAAGGTGCAATTTCCCTTTGGCGGAAACCCGt              >  1:487522/1‑251 (MQ=255)
                                                                                                                                                                                                                                              tGA‑TTTTCTTGCTGTTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCGCCGAATAAACAAAGGTGCAATTTCCCTTTGGCGGAAACCCGTAAACGGTTACAg  >  1:768266/1‑248 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             
CAATTCCCACGTGCGTTTTGCCGCCTGAATGCCGGCAATGCGGGCGGTGGCGAACACATCGCCTTTGTGATGCTGACCGGACAGGATCATTTGTAAGGTTTCGGGCGCCATACAAACATAAGCCTCAGCGCGGGCTTCACGGACGGTTTCGGCTTT‑GGCGGAAACGTCCACCATATTGGCTTCCCC‑ATGGG‑AA‑TTAAT‑GTGAGTAAATGTGACCATGATAATTAAAGTGCGGTTGA‑TTTTCTTGCTATTTTTC‑AAAAT‑C‑GGCGTGCATTAACCGCCGATACTTGCCAAATGACTACGTACGCCACTATCGCCGATGTGCAGATAATGATGTTCACGCTTACCTTGTAACGCAGAAAAAATTCTTGATTGCAAAATGAGTTGCTGCTTATCGTATTGTAACAAATCGCGTAATTCAATGCCTTCCTCGCCGAATAAACAAAGGTGCAATTTCCCTTTGGCGGAAACCCGTAAACGGTTACAG  >  1/1538805‑1539285

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