Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,151,322 A→G 15.9% A592A (GCT→GCC)  ppc ← phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,151,3220AG15.9% 145.8 / 3.7 69A592A (GCT→GCC) ppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base A (41/14);  new base G (6/6);  total (49/20)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCCGCGGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGG  >  NC_000913/4151088‑4151557
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                
gCCGCGGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑g                                                                                                                                                                                                                                            >  2:786226/1‑239 (MQ=255)
         gACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTgcggcg                                                                                                                                                                                                                                    >  2:825669/1‑238 (MQ=255)
              cGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGGCAGCAGCGGCGCATGGG‑CA‑GGTGCGCCGCcgcg                                                                                                                                                                                                                               >  2:937768/1‑238 (MQ=255)
                      aTGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCACCAGGTTGGCTTCCAGAATCGCCCCGGCATAAAGCGACAGGCTGCTGACGGTGTTTTCTGGCAGACCATATTTAAA‑GCGGATCATCCCGCCCTGTTCGGGTACGCGCAGCCCGCCTTTCAGGCCTCCTGGCGGGTGTGACCGCATCCCCCCATTGT‑CA‑GGCGcgccccccccgccacacgg                                                                                                                                                                                                                      >  1:1025937/1‑227 (MQ=255)
                        gactgacAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTTTGACAGCAGCGCCGCATGGG‑CA‑GGTGCGCCGCCGCCACCAA‑TGGa‑a                                                                                                                                                                                                                   >  2:40720/1‑238 (MQ=255)
                        gactgacAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGCTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGa‑a                                                                                                                                                                                                                   >  2:810727/1‑238 (MQ=255)
                           tgacAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGGGACAGCAGCGCCGCATGAG‑CA‑GGGGCGCCGCCGCGACCAA‑TGGA‑Accg                                                                                                                                                                                                                >  2:496399/1‑238 (MQ=255)
                             acAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGACAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTCGCGGGTGTGGCAGCAGCCCCGCCTCGG‑CA‑AGATCGCCGCCCCGGCCAA‑TGGacacacc                                                                                                                                                                                                               >  1:111395/1‑232 (MQ=255)
                                                aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTTTGACAGCAGCGCCGCATGAG‑CA‑GGGGCGCCGCCCCGACCAA‑TGGG‑AACGCCGCGGCCGGGGACCAACGTCa                                                                                                                                                                                          >  2:268179/1‑239 (MQ=255)
                                                aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGGGGGACAACGTc                                                                                                                                                                                           >  2:124247/1‑238 (MQ=255)
                                                 cGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGGG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGGACAACGTcc                                                                                                                                                                                          >  2:894230/1‑237 (MQ=255)
                                                 cGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGCG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGG‑ACCGCCCCGACCGCGGAACAACGTCa                                                                                                                                                                                          >  2:378021/1‑238 (MQ=255)
                                                             ttCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATAacc                                                                                                                                                                               >  2:70284/1‑237 (MQ=255)
                                                                 gCTCCGGCGGTGGCAGCAGGTTGGATTCCAGCATCGCCCCGGTATCAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCCGTTCGGTTACTCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAACAGCCCATCAG‑CA‑CGCGCGCCGGCCGGAACAA‑TGGG‑ACCGCCGCGCCCCTGGAACAACGGACGCTCGACACCCGCtgtt                                                                                                                                                                         >  1:575590/1‑239 (MQ=255)
                                                                     gggcggTGGCGGGAGTTTGGCTGCAAAACCCCCCCGGGTAAAAGGCGACGGGCTGGTGGCGGTTTTTTGTGGCATA‑CATTTTTAAAGGCGGATCATCCCCCCCTTTGCGTTTCCGCGGGGGCCCCTTTTCGGGTTTCGGGGGGGTTGGAAAGGCGGCGCC‑CATGGG‑GAGGGGGCCCCCCCGCCACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCg                                                                                                                                                                      <  1:236658/237‑1 (MQ=255)
                                                                     cggcggTGGCAGGGGGTTGGCTTCCAGAATCCCCCGGGTAAAACGCGACAGTCTGCGGCGGGTGTTTTCGGGCAGACCATTTTAAAG‑GCGGTTCTTCTCGCCTTGTTCGGTCACGCGCGGGCCGCCTTTCGGGCTTCCGGGGGGTTGTGACCGCAGCCCCGCAGGAG‑CG‑GGCGCCCCCCCGCGACCAA‑TGAA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCg                                                                                                                                                                      <  1:769381/238‑1 (MQ=255)
                                                                      ggcggTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGGCAGCACCGCCGCATGAG‑CA‑GGGGCGCCGCCGCGGCCAA‑TGGA‑ACCGCCGCCGCCGTGGGACAACGTCAGCCCAAAACCCCCTTTTTCGCa                                                                                                                                                                    >  2:621207/1‑239 (MQ=255)
                                                                      ggcggTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGGTTTGACCGCAGCGCCGCATGGGACA‑GG‑GCGCCGCCCCCACCCA‑TGGA‑ACCCCCCCGGCCCTGGAACCACGTCAGCCCAAAATCCCCttttttccc                                                                                                                                                                    >  1:249472/1‑235 (MQ=255)
                                                                               tggagtttGGCTGCCAAAGTCGCCCTGTTAAAAGGCGAGGGGCGGCGGCCGGTGTTTTCTGCCAGCCCATTTTTAAA‑CGGGTTCACCCCCCCTGGTTCGGTAGCGCCCGGCCCTCTTTCCGGGCTCCCTGGCGTTTGTGAGGGCGGCGCCGCATGAG‑AG‑GGTGCCCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGa                                                                                                                                                            <  1:172746/232‑1 (MQ=255)
                                                                                 cagGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCCCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCCATACCCGCCTTTTCGCCGGTTTTGATTa                                                                                                                                                         >  2:663099/1‑239 (MQ=255)
                                                                                       ggCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCa                                                                                                                                                    >  2:299960/1‑238 (MQ=255)
                                                                                                     ccccGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGatat                                                                                                                                     >  2:543359/1‑239 (MQ=255)
                                                                                                            aaaaaGCGAGAGGCTGCGGAGGGTGATCTCTGTCAGACATTTTTTAAA‑CCGGTCCATCCCGCTGTTTCCGTTCACCCGCAGGCCCCCTTTCAGGCTTCCTGGCGGTTGGGACCGCAGCGCCGCAGGAG‑AG‑GGGGCGCCCCCGCAACAAG‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGccc                                                                                                                             <  1:360996/237‑1 (MQ=255)
                                                                                                                    aCAGGCTGCTGACGGTGATTTCTTGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGCTGTGAAAGCAGCGCCGCATGAG‑CC‑AGTGCGCCGCCTCGCCCAA‑TGGA‑AACGCCGCGACCGTGGAACAACGGCAGCTCAATACCCGCCCTTTTGCCGGGTTTGCTTAAATCCTCCTTCGCCTGATCTTCGCCCCCAgaa                                                                                                                         >  1:594373/1‑234 (MQ=255)
                                                                                                                    aCAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGaa                                                                                                                        >  2:108804/1‑236 (MQ=255)
                                                                                                                     cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGTCAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCt                                                                                                                     >  2:738891/1‑238 (MQ=255)
                                                                                                                     cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCt                                                                                                                     >  2:820179/1‑238 (MQ=255)
                                                                                                                     cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCCCATGAG‑CA‑GGGGCGCCGCCCCGACCAA‑TGGA‑ACCGCCGCGGCCGGGGAACAAAGTCATCTCAATACCCGCTTTTTCTCAAGTTTTTATTAATGCATCCTTGTCCTGAGATT‑GCGCCCAGGAAgat                                                                                                                     >  1:959667/1‑236 (MQ=255)
                                                                                                                     cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAACAGCGCCGCATGAG‑CA‑GGTGCCCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCACAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCt                                                                                                                     >  2:286352/1‑238 (MQ=255)
                                                                                                                     cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGCTGTGACCGCAGGGCCGCATGCA‑CC‑AGTGCGCCCCCGCGGCCAA‑TGGA‑ACCCCCCGGACCGTGGAAACACATCACGTTCATCCCCGCTTCTTCGCACGGTTTGTTTAATGCGACCTGTCCCTGttaagccacccagccatct                                                                                                                      >  1:722626/1‑222 (MQ=255)
                                                                                                                        gctgctGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCa                                                                                                                 >  2:448614/1‑239 (MQ=255)
                                                                                                                               aCGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTc                                                                                                           >  2:828524/1‑238 (MQ=255)
                                                                                                                                  gTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGAAAGCCGCGCCGCATGAG‑CG‑GGGGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGCGGAAAAACGTCAGCTCAATACCCGCTTTTTCCCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCg                                                                                                        >  1:733390/1‑238 (MQ=255)
                                                                                                                                                  cATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATg                                                                                                                                                      >  2:1020930/1‑177 (MQ=255)
                                                                                                                                                    tatTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGa                                                                                     >  2:837849/1‑239 (MQ=255)
                                                                                                                                                         aaa‑GCGGCTCCTCTCCCCGTGTTGGTTGCCGCGCAGCCCCCCTTTCGGCTTTCGTGGCGTTTGGGGCGGGAGCGCCGAAGGGG‑CA‑GGTGCCCCCCCCCGCCCAA‑GGGA‑TCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGc                                                                                <  1:738891/239‑1 (MQ=255)
                                                                                                                                                          aa‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGcc                                                                               >  2:284041/1‑239 (MQ=255)
                                                                                                                                                           a‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGGGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTTCCATGACTCCCGCACCTTTTGCTGCGTCGGAATAGCCa                                                                              >  2:919220/1‑239 (MQ=255)
                                                                                                                                                             gCGGATCATCTCGCCCTGTTCGGTTACGCTCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATAACCGCTTTTTCGCAGGTTTTGTTTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAACTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCaa                                                                             >  1:997947/1‑239 (MQ=255)
                                                                                                                                                             gCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCCACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCCGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCaa                                                                             >  2:199068/1‑239 (MQ=255)
                                                                                                                                                               ggATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCCGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCCGGAAGCTGCCATCACTCCCCCCTCTTTTGCTGAGTCGGAATAGCCAatc                                                                           >  1:107746/1‑239 (MQ=255)
                                                                                                                                                                 atcacccccccttGTTGGGTTCCGCGCGGGCCGCCTTTCGGGCTTCCGGGCGGTGGTGACAGCGCCGCCGCGTGAG‑CG‑GGTGCGCCGCCGCGACCAA‑TGAA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatcat                                                                         <  1:284041/227‑1 (MQ=255)
                                                                                                                                                                      cccccccTTTGCGTTTCCGCGCAGCCCCCCTTCCAGGCTCCCTGGGGTTTGGCACAGCAGCGCCGCTGGAG‑CA‑GGTCCCCCCCCGCCACAAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCa                                                                    <  1:512297/235‑1 (MQ=255)
                                                                                                                                                                          cccTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCCCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGT‑ACCGCCGCGACCGTGGCACAACGTCAGCTCAATACCCGCTTTTTCGCAGGGTTTGATTAATGCTTCCTGTGCCTGATATT‑GCGCCCCGGACGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCa                                                                    >  1:606448/1‑235 (MQ=255)
                                                                                                                                                                           ccTGTTGGGTCACCGGGAGCCCGCCTTTCGGGTCCCGGGGGGGTTGTGACAGCAGCGGCCAATGAG‑CA‑GGGGCGCCGCCGCGCCCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGt                                                               <  1:825669/239‑1 (MQ=255)
                                                                                                                                                                           acTGTGCGGTTCCGCGCAGGCCGCTTTTGGGGTTCCGTGGGGGTGGGGACAGCAGCGCCGCAGGAG‑AA‑GGTGCGCCGCCGCAACCAA‑GGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGt                                                               <  2:111395/238‑1 (MQ=255)
                                                                                                                                                                              gTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGttt                                                             >  2:933111/1‑238 (MQ=255)
                                                                                                                                                                               ttCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTg                                                            >  2:215190/1‑238 (MQ=255)
                                                                                                                                                                               ttCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTg                                                            >  2:1191112/1‑238 (MQ=255)
                                                                                                                                                                               ttCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGGGCCTGATATT‑GCGCCCAGGAAGATGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTg                                                            >  2:256064/1‑238 (MQ=255)
                                                                                                                                                                                  tgtgaataggggccagccTTCCAGCCTTCCGGGGGGTTGGGACAGCGGCCCCGCGTGAG‑AG‑GGTGCCCCCCCGCGACAAA‑GGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGcc                                                          <  1:117466/222‑1 (MQ=255)
                                                                                                                                                                                       cgcgCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATc                                                   >  1:362766/1‑239 (MQ=255)
                                                                                                                                                                                         cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAg                                                 >  2:890169/1‑239 (MQ=255)
                                                                                                                                                                                         cgcAGGCCGCCTTTCAGGCTTCCTGGCGGGTGTGACAGCAGCGCCCCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑AGGA‑ACCGCCGCGACCGTGGAACAACGCCAGCTCCATACCCGCTTTTTCGCCGGTTTTTGTTAATGCATCCTGGGCCTGATATT‑GCGCCCAGCACGCTGCCATCACTCCCGCATCTTTTGCTCAGTTGGAGTAGCCCATCCTCACCCTCTGTTTGACCTGAATCa                                                  >  1:880596/1‑238 (MQ=255)
                                                                                                                                                                                          ggaGGCCGCCTTTGAGGCTTCCTGGCGGTTGAGACAGCAGCGCCGCATGAG‑GA‑GGTGCGCCTCCGCGACAAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAg                                                 <  2:1025937/236‑1 (MQ=255)
                                                                                                                                                                                              cccgccTTTCCGGCTCCCGGGCGTTTGGGCCGGCAGCGCGGAATGGG‑CA‑GGTGCGCCCCCGCGACCAA‑TGGA‑CCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCa                                             <  1:215190/237‑1 (MQ=255)
                                                                                                                                                                                              acccctgtagtgtctccgtGGCGGTTGTGCCAGCAGCGCGGTAGGGG‑CA‑GGTGCGCCCCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCa                                             <  1:837849/220‑1 (MQ=255)
                                                                                                                                                                                                  ccTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGAt                                         >  2:1057350/1‑238 (MQ=255)
                                                                                                                                                                                                      tCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCCA‑TGGA‑ACCGCCCCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGGAGGTTTTGATTAATGCATCCTGTGCCCGATATT‑GCTCCCAGGAAACTGCCATCACTCCCGCATATTTTGCTGAGTCGGAAAAGCCCATCATCACCAACTGTTTGCCCTGAATCAGGCCCCGATAcac                                     >  1:704497/1‑236 (MQ=255)
                                                                                                                                                                                                       ctggCTTCCGGGGGGTTGTGACAGAAGCGCCGCAGGAG‑AA‑GGTGCGCCGCCGCGACAAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                   <  2:798146/237‑1 (MQ=255)
                                                                                                                                                                                                       cgggCTTCGGGGCGTTGGGGACAGGCGCGCCGCATGAG‑GG‑GGGGCGCCGCCGCAACCAA‑TGAA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                   <  1:543359/237‑1 (MQ=255)
                                                                                                                                                                                                       caggccgcgggggggtggagACAGGCGCGGCGCAGGAG‑CG‑GGTGCCCCGCCGC‑ACCAAGGGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                   <  1:1057350/220‑1 (MQ=255)
                                                                                                                                                                                                             tCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAgg                                                >  2:117466/1‑220 (MQ=255)
                                                                                                                                                                                                                      tgggaagagcggcGCCCAAGGAG‑CA‑GGTGCCCCGCCGCCCCAAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGagcagc                     <  2:704497/231‑1 (MQ=255)
                                                                                                                                                                                                                       ggtgccggcagcCCCGCATAGG‑AG‑GGTGCCCCGCC‑CGACCAA‑AGGACACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCACAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCt                    <  1:806849/231‑1 (MQ=255)
                                                                                                                                                                                                                         tgACCGCAGCCCCGCATGAG‑CA‑GGTGCGCCCCCGCGGCCAA‑TGGA‑ACCGCCGGGCCCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGAGCCTGATATT‑GCGCCCAGGAAGCTGCCATCCATCGCGCATCCTTTGACGCGTCGGACTAGCCAATCATCATCCTCTGTTTGCCCTGAATCAGGCCACGATACTAGTCACTATTTAGCAGCCggg                 >  1:1010164/1‑239 (MQ=255)
                                                                                                                                                                                                                                agcGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCCTCACCATCTGTTTGCCCTGAATCAGGCCACGATCCCAGTCAATTTTGAGCAGCTTGGTCCTGAc          >  1:977495/1‑239 (MQ=255)
                                                                                                                                                                                                                                     cgcAGGGG‑CA‑GGTGCCCCGCCGCGACCAA‑TGAA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGt     <  2:107746/239‑1 (MQ=255)
                                                                                                                                                                                                                                        aTGAG‑CG‑GGTGCGCCGCCGCGTCCAA‑GGGG‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTgg  <  2:575590/239‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                
GCCGCGGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAA‑TGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGG  >  NC_000913/4151088‑4151557

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