Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,909,229 A→C 18.5% G146G (GGT→GGG)  pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,909,2290AC18.5% 132.4 / 6.3 70G146G (GGT→GGG) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base A (31/24);  new base C (3/12);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGA‑CCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAACGCTCGTAGTGCCCCAGGTCCAGGTCGGTTTCA  >  NC_000913/2908992‑2909466
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
tGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATGCCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCa                                                                                                                                                                                                                                                 >  1:808096/1‑238 (MQ=255)
tGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAAAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTCCGGGGAACGGCAAGGATTCGATATCa                                                                                                                                                                                                                                                 >  1:34501/1‑238 (MQ=255)
tGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGGGCATAAACAGAGTGTGCTCCCGGCCAATTTCAACAGCCATCTGGCGAATCGCCTCGAGGAACGGCAAGAATTCGAAAacac                                                                                                                                                                                                                                                >  2:167238/1‑239 (MQ=255)
     tCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGCCAGGGGCATAAACAGAGGGTGCTCACGGCCAATTTCAACAGCCACCTGGCGAATCGCTTCGAGGAAAGGCAAGGATTCGATAACACCt                                                                                                                                                                                                                                              >  2:1005796/1‑236 (MQ=255)
          gCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGGGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGGCGGGATTCGATATCaaaaaaggtacc                                                                                                                                                                                                                                      >  2:916545/1‑228 (MQ=255)
            ccgttcgttCCCCGGAACAGCGCGATCTGAACGACAAATCGGGATGTCAGGCTGGATACCGAGGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTAccgc                                                                                                                                                                                                                                    <  1:167238/238‑1 (MQ=255)
                tcttgccccGTAACCGCGGTATCTGAACGAAAATTCAGGGTGTCAGGCTGGATCCCGAGGGGGAGCATCTTTTTTACAGGGGTCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGGGCATAAACAGAGTGTGCTCACGGCAAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGAt                                                                                                                                                                                                                                <  2:319546/232‑1 (MQ=255)
                 cgttcgCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGGGACTCGATTCAGGGAACGGCAAGGATTCGATATCCCCCACTGTAccgccgtct                                                                                                                                                                                                                               >  2:684581/1‑236 (MQ=255)
                 cgttcgCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGAtt                                                                                                                                                                                                                               >  1:732996/1‑239 (MQ=255)
                             tcgcgggttCTGAACAAAAAAGCAGGGTGGAAGGTTGGAATCCTGTGGAGATCGAGTTTTTTACAGGGTGCGGATGTGGTTTGGTTTTGTCTACCCATGGCGCTGCCTGGTACGTCGCCAGGGCCAGGTGCAAAAACGGAGTGGGCTCCCGGTCAATTTCAACAGCCATCTGGGAATTCGTTTCTAGGAACGGTAAGGATTCGATATCACCAACTGTACCGCCGATTTCTACCAGtact                                                                                                                                                                                                                   <  2:736382/231‑1 (MQ=255)
                                cgcgATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGGGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGCATCGCTTCGGGGAACGGCAAGGATTCGAAATACCCCACTGTACCCCCGATTTTTACAAggcgaacg                                                                                                                                                                                                                >  2:977924/1‑232 (MQ=255)
                                cgcgATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGAAAACACCAACTGTACCGCCCATTTCTACCAgacatac                                                                                                                                                                                                                 >  2:488668/1‑232 (MQ=255)
                                cgcgATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTGGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGACAGTGCTCCCGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGGGGAACGGCAAGGACTCGGTATCACCAACCGAACCCCCGATTTCTACCaagacaaaa                                                                                                                                                                                                                >  2:14211/1‑231 (MQ=255)
                                 gctaTCTAAACGACAAATCAGGATGTCAGGCTGGATCCCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGtact                                                                                                                                                                                                                   <  1:122370/232‑1 (MQ=255)
                                 gcgATCTGAACGACAAATCGGGCTGTCAGGCTGGATCCCGAGGGAGGCAAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGACATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGtact                                                                                                                                                                                                                   <  1:997892/235‑1 (MQ=255)
                                     tCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTCCTGGACCGCCGATTTCTTCCAGTACGACGTCCTg                                                                                                                                                                                                           >  1:872132/1‑239 (MQ=255)
                                       tGAAAGAAAAAGCGGGTTGTCGGGGTGGATGCCGATGGAGAGCCGCTTTTTTACGGGGTGCTGAGTCGGTTTGGTTTTGACTCCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGGGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGAc                                                                                                                                                                                                         <  2:808096/239‑1 (MQ=255)
                                            gACAAATCAGGTTTTCAGGGTGGATACCGAGGGAGAGCAGCTTTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCCCCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCg                                                                                                                                                                                                    <  2:139034/239‑1 (MQ=255)
                                               aaaTCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGCACGGCAAGGATTCGCTATCACCTACTGTACCGCCGATTTCTACACGTACTACGTCATGACCTTCAcc                                                                                                                                                                                                  >  1:1046863/1‑238 (MQ=255)
                                                   cAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCACCGTCAGGTGCATACACAGAGTGTGCTCACGGCCAACTTCAACAGCCATCTGGAGAAACGCTTCGAGGAAAGGCCAGGAATCGACATCCACACCTGTACCCCCGCTTTACACCAGCACAACATCCAGACCTTCCCCACCtt                                                                                                                                                                                             >  2:2972/1‑239 (MQ=255)
                                                         ggcAGGTTGGATACCGATGGAGGGCAGCTTTTTTACAGAGTTCTGAGTCGGTTTGGTTTTGAC‑TCACCGAACCCTGCCATGAACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCCCGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCa                                                                                                                                                                                       <  1:576075/236‑1 (MQ=255)
                                                           ggggggggATACCGGTGGAGGGCCGCCTTTTTACAGGGGTCTGAGTCGGTTTGGTTTTGACTTCACCAGGCGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                                     <  1:371566/233‑1 (MQ=255)
                                                           cctgCTGGATACCGTTGGAGAACCTCTTTTTGACAGGGTGCTGGGGCGGTTTGTGTTTGACTCCACCAGCCGCCGCCATGTACGGCCCCAGCGGCAGGGGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                                     <  2:830850/236‑1 (MQ=255)
                                                           cAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCCCCAGACGCTGCCATGTCCGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                                     <  1:1005796/239‑1 (MQ=255)
                                                           cAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTTTACCGCCGATTTCTACCAGTACATAGGCATACCCCTCGCCACCTTCCAGCAc                                                                                                                                                                                      >  2:129548/1‑238 (MQ=255)
                                                            aGGTTGGATACCGATGGAGGGCAGCTTTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCCCCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                                     <  1:916545/238‑1 (MQ=255)
                                                             ggttGGTTACCGGGGGAGGGCAGTTTTTTTAAAGGGGGCGGGGTGGGTTTGGTTTTGACTTCACCAGACGCTGCCAGGTACGGCACAAGCGTCAGGTGCATAAACAGGGTGTGTCCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCATGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcgcg                                                                                                                                                                                   <  2:43513/236‑1 (MQ=255)
                                                                  gATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCCTGACCTTCGCCACCTTCCCACACGCGCCCtt                                                                                                                                                                              >  2:720907/1‑239 (MQ=255)
                                                                         aTGGAGAGCAGCTTTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTAc                                                                                                                                                                                                                          <  1:146300/188‑1 (MQ=255)
                                                                         aTGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTAc                                                                                                                                                                                                                          >  2:146300/1‑188 (MQ=255)
                                                                         aTGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCtt                                                                                                                                                                                                      <  1:551561/208‑1 (MQ=255)
                                                                         aTGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCtt                                                                                                                                                                                                      >  2:551561/1‑208 (MQ=255)
                                                                            gagaGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCATAACTACGGCATGACCTTCGCCACCTTCCAGCACGCCCCCTTTGATTGCGtt                                                                                                                                                                    >  2:305218/1‑239 (MQ=255)
                                                                            gagaGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATCTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGCTTTCTACCAGTACTACGTCCTGACCTTCGCCCCCTTCCAGCACGCGCCCTTTGATTGCGtt                                                                                                                                                                    >  2:1046502/1‑239 (MQ=255)
                                                                                 cagcTCTTTTACAGGGTGCTGGGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTAATTGCGTTAGTg                                                                                                                                                                <  1:45654/238‑1 (MQ=255)
                                                                                 cagcTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGGACCGCCGCTTTCTACCAGTACTACGGCATGACCTTCCCCACCTTCCAGCACACGCCCTTTGATTGCGTTAATGa                                                                                                                                                               >  2:144191/1‑239 (MQ=255)
                                                                                  agcTCTTTTGAGGGGTGCGGGGTGGGTTTGGTTTTGACTTCCCCCGGCGCTGCAATGTCCGGCACAGGCGGCAGGTGCATAAACAGTGTGTCCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTg                                                                                                                                                                <  2:989125/237‑1 (MQ=255)
                                                                                     tctTTTACTGATTGCTGAGTCGGTTTGGGTTTGACTTCTCCAGACGGTGCCTTTTACGGCACCCGCGTCAGGTGCTTAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTTCCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGAt                                                                                                                                                                          <  1:692069/224‑1 (MQ=255)
                                                                                     tctTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCCCCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGAtgtg                                                                                                                                                           <  1:874466/239‑1 (MQ=255)
                                                                                                gTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACGGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAAGCGCCCTTTGTTTGCGTTAGTAATTTGGGGAATaa                                                                                                                                                   >  2:410210/1‑236 (MQ=255)
                                                                                                gTGCTGAGGTGGTTAGGTTTTGACTTCACCAGACTCGCCCATGTGCGCCACCAGTGTCAGGTGCATAAACGGGGTGTGATCCAGGCCAATTTCAACAGCCATTTGGCGAATCGCTTCGAGGAACGGAAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATa                                                                                                                                                    <  1:684581/235‑1 (MQ=255)
                                                                                                  gctggttggttttggtttTGACTTCCCCAGACCGGGCCATGGACGGCCCCAGGTTGAGGGGCATAAACAGGGGGTGCCCACAGACAATTTCACCCACCATCTGGCGAACCGCTTCGGGGAACGGCAAGGATTCGATATCCCCTACTGTACCGGCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTg                                                                                                                                               <  2:680180/228‑1 (MQ=255)
                                                                                                  gCTGAGTCGGTTTGGTTTTGCCTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGGGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTTTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTg                                                                                                                                               <  1:816193/238‑1 (MQ=255)
                                                                                                        tgggtttggtttTGCCTTCCCCAGACGCTGCCATGTACGGCACCAGCGTAAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGt                                                                                                                                        <  1:305218/236‑1 (MQ=255)
                                                                                                           gtttggtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTCCTGTACCGCCGGTTTCTACCAGAACTACGGCATGACCTTCCCCACCTTCCAGCCCGCCCCCTTTGATTGCGTTAGTGCTGGGCGGAAAAACACGCAC‑GGGTgg                                                                                                                                     >  2:262455/1‑237 (MQ=255)
                                                                                                           gtttggtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACAGCAAGGATTCGATATCACCTACTGTACCGCCGCTTTCTACCAGTACTACGTCCTGACCTTGCCCAACTTCCCGCACCGGCTCCTTGCTTGCGTTGGGGATTTGCGGGATAACCCGGAC‑GGGTggc                                                                                                                                    >  2:241945/1‑237 (MQ=255)
                                                                                                                                 cGCTGCCATGTACGGCCCCAGCGGCAGGTTCATAAACAGAGGGTGCCCACGGCCAATTCCAACAGCCATCTGGCGAATCTCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCgcggc                                                                                                                 <  1:129548/236‑1 (MQ=255)
                                                                                                                                          gTACGGCCCAAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATCCCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCa                                                                                                     <  1:1046502/239‑1 (MQ=255)
                                                                                                                                          gTACGGACCCAGCGTCAGGTGCAAAAACAGAGTGGGCTCACGGCCAATTTCAACAGCCCTTTGGCGAATCGCTCCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCa                                                                                                     <  1:488668/239‑1 (MQ=255)
                                                                                                                                           tACGGCCCCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAg                                                                                                    <  1:806444/239‑1 (MQ=255)
                                                                                                                                                cACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGAt                                                                                                                                                              >  2:204613/1‑177 (MQ=255)
                                                                                                                                                cACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGAt                                                                                                                                                              <  1:204613/177‑1 (MQ=255)
                                                                                                                                                cACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGt                                                                                               >  1:107421/1‑239 (MQ=255)
                                                                                                                                                cACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGCGGTTCTTTACGCAGAACg                                                                                                >  1:261109/1‑238 (MQ=255)
                                                                                                                                                 aCCAGCGTCAGGTGCATAAACAGAGGGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGt                                                                                               <  1:291825/238‑1 (MQ=255)
                                                                                                                                                     gCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCagag                                                                                          <  1:247627/239‑1 (MQ=255)
                                                                                                                                                          aGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGt                                                                                                                                        >  2:946819/1‑188 (MQ=255)
                                                                                                                                                          aGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGt                                                                                                                                        <  1:946819/188‑1 (MQ=255)
                                                                                                                                                             tGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTCCGCAGAACGTCAGAGTAGAtcg                                                                                   >  1:834733/1‑236 (MQ=255)
                                                                                                                                                               ctaaaaCGGGGTGTGCCACCGGCCAATTTCCACCGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCCACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGA‑c                                                                               <  2:842300/236‑1 (MQ=255)
                                                                                                                                                                        agTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGCGTGGAAAGGA‑CCCGTGGTGa                                                                      >  1:958580/1‑239 (MQ=255)
                                                                                                                                                                         gtgtGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGGTCTTTACGCAGAACGTCAGATTAGATACGA‑CCCGGGGTGaa                                                                     >  1:480620/1‑239 (MQ=255)
                                                                                                                                                                          tgtgCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCACGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGA‑CCCGGGGTGAAg                                                                    >  1:284982/1‑239 (MQ=255)
                                                                                                                                                                           gtgCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGCGGTAGTCACCGCGGCGTTCTTTACGCCGAACGTCCGAGAGGATACGG‑CCCGGGGGGAAgg                                                                   >  2:983770/1‑238 (MQ=255)
                                                                                                                                                                              ctaaCGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGGCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTTCGTTGGAGATGTTCGGAATAATC‑GTACAGGTTGGGCCGAGGTAGTCACCGCCGCGTTTTTTTTCCAGACGGTCAGACTAGA‑ACGATCCAGTGGGGAAggtgt                                                                >  2:850220/4‑239 (MQ=255)
                                                                                                                                                                              cTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAGGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                                     <  2:83860/124‑1 (MQ=255)
                                                                                                                                                                              cTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                                     >  1:83860/1‑124 (MQ=255)
                                                                                                                                                                                  ggcataaTTACAACCCCCGTCGGGGGACTCGCGTTGAGGAGCGGGAAGGATTCACTTTCCCCTACTGGCCCGCCGATTTCTACCAGTGCTACGTCATGACCGTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCGCCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGA‑CCCGTGGTGAAGTTGTTgcg                                                            <  2:34501/234‑1 (MQ=255)
                                                                                                                                                                                  cGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGGAGTCACCGCGGGGTTCTTTACGCAGAACGTCAGAGTAGATAAGA‑CCCGGGGTGAGATTGTTgcg                                                            >  2:624441/1‑239 (MQ=255)
                                                                                                                                                                                             aaCAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGCCAGAGCAGAGCAGA‑CCCGCGGGGACGTTTTTGcggcggggacctt                                                 >  2:46060/1‑232 (MQ=255)
                                                                                                                                                                                                          cGAATCCTTTCGAGGAACGGCAAGGATTCGATACCCCCTACTTGCCCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGA‑CCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGa                                     <  1:983770/238‑1 (MQ=255)
                                                                                                                                                                                                            aaTCGCTTCGAGGAACGGCGAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTCCGGAATAACCTTCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGA‑CCCGCGGTGACGTTGTTGCGGGGGCCCAATTTTGTagacataaaa                                   >  2:210051/1‑229 (MQ=255)
                                                                                                                                                                                                            aaTCGCTTCGAGGAACGGCAAGAATTCGAGATCCCCTACCGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGA‑CCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGaaa                                   <  1:318869/238‑1 (MQ=255)
                                                                                                                                                                                                                                            cACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGA‑CCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACAAATGAAAAGCTCGTAGTGGCCCAGGTCACGGGCGGTTTCa  >  1:374605/1‑239 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
TGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCAC‑GGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGA‑CCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAACGCTCGTAGTGCCCCAGGTCCAGGTCGGTTTCA  >  NC_000913/2908992‑2909466

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: ‑