Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 308,054 (C)7→8 coding (1279/2526 nt) ecpC ← ECP production outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913308,0471.C96.8% 334.3 / 1.4 93coding (1286/2526 nt)ecpCECP production outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (44/46);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCTTTTT  >  NC_000913/307897‑308192
                                                                                                                                                       |                                                                                                                                                     
gCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCg                                                                                                                                                        >  1:318677/1‑151 (MQ=255)
   tCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑c‑c‑c                                                                                                                                                 >  2:1146871/1‑151 (MQ=255)
    cGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑c‑c‑cc                                                                                                                                                >  2:335957/1‑151 (MQ=255)
       tACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑c‑c‑ccccc                                                                                                                                             >  3:507495/1‑151 (MQ=255)
        aCGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑c‑c‑cccccc                                                                                                                                            >  2:1174956/1‑150 (MQ=255)
                attgattgCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCa                                                                                                                                    <  3:856453/151‑1 (MQ=255)
                 ttgattgCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGc                                                                                                                                     >  1:355282/1‑149 (MQ=255)
                    attgCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGt                                                                                                                                <  3:852690/151‑1 (MQ=255)
                    attgCAAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑C‑C‑CCCCCAAAGCGGCAGCg                                                                                                                                 >  2:499737/1‑150 (MQ=255)
                     ttgCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCg                                                                                                                                 >  1:223455/1‑149 (MQ=255)
                      tgCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTc                                                                                                                               >  2:386248/1‑150 (MQ=255)
                      tgCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTc                                                                                                                               <  2:223723/150‑1 (MQ=255)
                       gCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAATCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑C‑CTCCCCC‑AAGCGGCAGCGTc                                                                                                                               <  3:80329/149‑1 (MQ=255)
                       gCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTc                                                                                                                               >  2:1063422/1‑149 (MQ=255)
                       gCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCa                                                                                                                              <  2:245498/150‑1 (MQ=255)
                       gCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCa                                                                                                                              <  2:478103/150‑1 (MQ=255)
                        ccAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCC‑C‑CCCCC‑AAGCGGCAGCGTCAg                                                                                                                             <  3:710284/151‑1 (MQ=255)
                        ccAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTc                                                                                                                               >  2:1045265/1‑148 (MQ=255)
                         cAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGAc                                                                                                                           <  2:1053147/151‑1 (MQ=255)
                               gggTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGttt                                                                                                                     >  1:449031/1‑151 (MQ=255)
                               gggTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGttt                                                                                                                     >  1:131101/1‑151 (MQ=255)
                                ggTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTc                                                                                                                    <  3:824313/151‑1 (MQ=255)
                                 gTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTca                                                                                                                   <  4:235981/151‑1 (MQ=255)
                                   tttttCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTcac                                                                                                                  <  1:268838/150‑1 (MQ=255)
                                      ttCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTcacccac                                                                                                              >  2:997727/1‑151 (MQ=255)
                                       tCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTcacccacc                                                                                                             >  2:1265726/1‑151 (MQ=255)
                                          tGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGcc                                                                                                          >  2:164418/1‑151 (MQ=255)
                                          tGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGcc                                                                                                          >  3:320892/1‑151 (MQ=255)
                                            aTTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                                                                        <  3:386284/151‑1 (MQ=255)
                                            aTTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                                                                        >  4:199680/1‑151 (MQ=255)
                                             ttAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                                                                        >  2:1402853/1‑150 (MQ=255)
                                             ttAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                                                                        >  3:836152/1‑150 (MQ=255)
                                                 cccACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTAt                                                                                                   <  3:363594/151‑1 (MQ=255)
                                                 cccACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTAt                                                                                                   <  3:395087/151‑1 (MQ=255)
                                                   cacaGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCg                                                                                                 >  3:559232/1‑151 (MQ=255)
                                                     caGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTa                                                                                               <  3:571502/151‑1 (MQ=255)
                                                        cGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATcc                                                                                            <  2:57014/151‑1 (MQ=255)
                                                         gAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATcc                                                                                            >  1:166483/1‑150 (MQ=255)
                                                         gAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATcc                                                                                            >  2:404756/1‑150 (MQ=255)
                                                          aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCat                                                                                          >  4:44342/1‑151 (MQ=255)
                                                                   cTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTg                                                                                 >  2:1043892/1‑151 (MQ=255)
                                                                   cTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTg                                                                                 >  2:469208/1‑151 (MQ=255)
                                                                     cccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgcc                                                                               <  1:465207/151‑1 (MQ=255)
                                                                     cccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgcc                                                                               <  2:1364407/151‑1 (MQ=255)
                                                                     cccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTg                                                                                 <  2:303473/149‑1 (MQ=255)
                                                                     cccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTg                                                                                 <  1:120512/149‑1 (MQ=255)
                                                                           agagTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGccc                                                                           >  2:443888/1‑149 (MQ=255)
                                                                              gTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCCCTGGCCCGCATATTTTGCAGGTTAACGTTGATCGCCC‑C‑CCCCC‑AAGCGGCAGCGTCAGCCGGGTTTCACACACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCCGCt                                                                       <  4:45329/151‑1 (MQ=255)
                                                                              gTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCt                                                                       <  3:331990/150‑1 (MQ=255)
                                                                              gTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                                                                      <  3:1146927/151‑1 (MQ=255)
                                                                                ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTaa                                                                     <  2:155502/150‑1 (MQ=255)
                                                                                ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                                                                      >  4:249661/1‑149 (MQ=255)
                                                                                ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                                                                      >  4:134762/1‑149 (MQ=255)
                                                                                ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAg                                                                    <  2:243588/151‑1 (MQ=255)
                                                                                 gcgcTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGTATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGc                                                                   <  3:1112876/151‑1 (MQ=255)
                                                                                  cgcTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCg                                                                  >  4:252629/1‑151 (MQ=255)
                                                                                  cgcTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCg                                                                  >  1:374843/1‑151 (MQ=255)
                                                                                     tGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTAc                                                               >  3:1346823/1‑151 (MQ=255)
                                                                                      gATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑CAC‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCAACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCg                                                                  <  2:266564/147‑1 (MQ=255)
                                                                                         gctgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCaa                                                           <  3:773363/151‑1 (MQ=255)
                                                                                         gctgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCa                                                            <  3:223041/150‑1 (MQ=255)
                                                                                            gcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAAt                                                          <  2:175161/149‑1 (MQ=255)
                                                                                             ccGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATg                                                         >  3:82844/1‑149 (MQ=255)
                                                                                             ccGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAg                                                       <  3:702670/151‑1 (MQ=255)
                                                                                               gATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGcc                                                     >  3:251099/1‑151 (MQ=255)
                                                                                               gATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGcc                                                     >  2:661992/1‑151 (MQ=255)
                                                                                                    tgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAgg                                                >  3:1359351/1‑151 (MQ=255)
                                                                                                     gctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGt                                               <  2:406117/151‑1 (MQ=255)
                                                                                                      ctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTc                                              >  2:316045/1‑151 (MQ=255)
                                                                                                        ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTc                                              <  3:518794/149‑1 (MQ=255)
                                                                                                        ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCg                                             >  3:973151/1‑150 (MQ=255)
                                                                                                         cATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCg                                             >  3:903736/1‑149 (MQ=255)
                                                                                                         cATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAg                                           <  2:361777/151‑1 (MQ=255)
                                                                                                          aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAgg                                          <  1:193089/151‑1 (MQ=255)
                                                                                                           tGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGc                                         >  3:273463/1‑151 (MQ=255)
                                                                                                                 gTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTg                                    <  3:914538/150‑1 (MQ=255)
                                                                                                                   cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCt                                  <  4:31873/150‑1 (MQ=255)
                                                                                                                   cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTa                                 <  3:1045311/151‑1 (MQ=255)
                                                                                                                       ggCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCATCGTCAGACGTGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGc                               <  2:718444/149‑1 (MQ=255)
                                                                                                                       ggCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGc                               <  4:394176/149‑1 (MQ=255)
                                                                                                                        gCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCaa                            <  1:301468/151‑1 (MQ=255)
                                                                                                                         ccAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCaa                            >  4:266496/1‑150 (MQ=255)
                                                                                                                          cAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAAct                          >  1:47029/1‑151 (MQ=255)
                                                                                                                            gCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAActct                        >  2:925893/1‑151 (MQ=255)
                                                                                                                            gCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAActct                        <  4:253446/151‑1 (MQ=255)
                                                                                                                            gCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAActct                        <  4:140908/151‑1 (MQ=255)
                                                                                                                            gCATATTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAActct                        >  1:251840/1‑151 (MQ=255)
                                                                                                                              atatTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAActctct                      >  1:7395/1‑151 (MQ=255)
                                                                                                                              atatTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAActct                        <  3:499772/149‑1 (MQ=255)
                                                                                                                                atTTTGCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCttt                    >  3:145855/1‑151 (MQ=255)
                                                                                                                                     gCAGGTTAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCt                <  3:889597/150‑1 (MQ=255)
                                                                                                                                          ttAACGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTcgc          >  3:1338752/1‑151 (MQ=255)
                                                                                                                                            aaCGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGt        <  3:404787/151‑1 (MQ=255)
                                                                                                                                             aCGTTGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTc       >  2:589713/1‑151 (MQ=255)
                                                                                                                                                 tGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCtttt   >  2:472775/1‑151 (MQ=255)
                                                                                                                                                 tGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCtttt   >  3:755772/1‑151 (MQ=255)
                                                                                                                                                 tGATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCttt    <  1:75143/150‑1 (MQ=255)
                                                                                                                                                  gATCGC‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCttttt  <  3:1175012/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑C‑C‑CCCCC‑AAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCTTTTT  >  NC_000913/307897‑308192

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