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.C95.7% 238.2 / 3.6 70coding (1286/2526 nt)ecpCECP production outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (33/34);  total (34/36)
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.80e-01

CTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCTTTTT  >  NC_000913/307903‑308192
                                                                                                                                                 |                                                                                                                                                   
cTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑c‑ccccc                                                                                                                                             >  1:673991/1‑151 (MQ=255)
 tACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑c‑cccccc                                                                                                                                            >  1:315423/1‑151 (MQ=255)
   cGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑c‑ccccccc                                                                                                                                           >  2:738078/1‑149 (MQ=255)
     tCGCCATTGATTGTAAATGCGGGTCTTTTCCTCACAAACCACCAGCGTCTTAATGCTTCCCGGTCGCGTGGCGCTGATTATGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGAATCATTTGAACGATAACGTCGATCggcacctcga                                                                                                                                             >  2:207630/1‑140 (MQ=255)
      cGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑c‑cccccccaagc                                                                                                                                       <  1:1036899/151‑6 (MQ=255)
      cGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑c‑cccccccaag                                                                                                                                        <  1:675029/150‑5 (MQ=255)
      cGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑C‑CCCCCCAAGCg                                                                                                                                       <  2:542522/151‑1 (MQ=255)
        cAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCgg                                                                                                                                      <  1:597024/151‑1 (MQ=255)
         aaTTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGc                                                                                                                                     >  2:489908/1‑151 (MQ=255)
          attgattgCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCa                                                                                                                                    <  1:1080309/151‑1 (MQ=255)
             gattgCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCg                                                                                                                                 >  1:966665/1‑151 (MQ=255)
                 gCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTc                                                                                                                               >  2:776671/1‑149 (MQ=255)
                 gCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAg                                                                                                                             <  2:1141264/151‑1 (MQ=255)
                       gCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACggg                                                                                                                        <  1:125399/150‑1 (MQ=255)
                         gggTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGttt                                                                                                                     >  1:852300/1‑151 (MQ=255)
                         gggTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGttt                                                                                                                     >  2:872452/1‑151 (MQ=255)
                          ggTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGttt                                                                                                                     >  2:926585/1‑150 (MQ=255)
                                   cTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGc                                                                                                           >  2:737511/1‑151 (MQ=255)
                                     gATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGc                                                                                                           <  1:764546/149‑1 (MQ=255)
                                      aTTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                                                                        >  1:814384/1‑151 (MQ=255)
                                         aaCCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGAtt                                                                                                     <  1:964866/151‑1 (MQ=255)
                                         aaCCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGAtt                                                                                                     >  1:1216809/1‑151 (MQ=255)
                                          aCCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGAtt                                                                                                     >  1:143517/1‑150 (MQ=255)
                                           cccACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGTGTCAGACGGGTTTCACCCACCGCCTGATTAt                                                                                                   <  1:1129003/151‑1 (MQ=255)
                                           cccACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTAt                                                                                                   >  1:285193/1‑151 (MQ=255)
                                             cacaGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCg                                                                                                 <  2:690026/151‑1 (MQ=255)
                                                aGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTAt                                                                                              >  1:2997/1‑151 (MQ=255)
                                                    aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCat                                                                                          >  1:1218951/1‑151 (MQ=255)
                                                    aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCat                                                                                          <  1:634277/151‑1 (MQ=255)
                                                    aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCat                                                                                          >  1:804812/1‑151 (MQ=255)
                                                        aaaGCGTCCCGGTAGAGTGGCGCTGATGCTGCCGATAATGCTGCGTGAGCTGTCGCTGGCCCTCATATTTTGGAGGTTAACGTTGATCT‑‑CGCCCCCCAGGCGCCAGCGCCTGACGGGCCTCACCCACCATCTGATTATCGCATCCATAGcc                                                                                      >  2:525632/1‑151 (MQ=255)
                                                        aaaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCatat                                                                                        <  1:924944/149‑1 (MQ=255)
                                                         aaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGTGGTGTCACTGGCCGGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCatat                                                                                        <  2:322984/148‑1 (MQ=255)
                                                          aGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCg                                                                                    <  2:633007/151‑1 (MQ=255)
                                                          aGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGt                                                                                   <  1:50917/151‑1 (MQ=255)
                                                            ccTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                                                                  >  1:661807/1‑151 (MQ=255)
                                                                ccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccg                                                                              >  1:623078/1‑151 (MQ=255)
                                                                     agagTGGCGCTGATGCTGCCGATGCTGCTCCATGAGTTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCa                                                                          <  2:230908/150‑1 (MQ=255)
                                                                      gagTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGc                                                                        <  1:1098682/151‑1 (MQ=255)
                                                                       agTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCt                                                                       <  2:1208196/151‑1 (MQ=255)
                                                                        gTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                                                                      <  1:692351/151‑1 (MQ=255)
                                                                         tGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCTGCCCAGCTa                                                                      >  2:1210662/1‑150 (MQ=255)
                                                                         tGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                                                                      >  2:806446/1‑150 (MQ=255)
                                                                         tGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                                                                      >  2:895281/1‑150 (MQ=255)
                                                                         tGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATAGTATCCATATCCCGTTGCCGCCCAGCTa                                                                      >  2:803458/1‑150 (MQ=255)
                                                                               tGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTAc                                                               >  1:29135/1‑151 (MQ=255)
                                                                                    ctgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCaa                                                           <  1:293599/150‑1 (MQ=255)
                                                                                      gcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa                                                        <  1:168122/151‑1 (MQ=255)
                                                                                      gcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa                                                        <  1:51397/151‑1 (MQ=255)
                                                                                      gcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa                                                        <  1:388333/151‑1 (MQ=255)
                                                                                        cGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc                                                      >  2:1142725/1‑151 (MQ=255)
                                                                                         gATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGcc                                                     >  2:500751/1‑151 (MQ=255)
                                                                                          aTGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCg                                                    <  2:990541/151‑1 (MQ=255)
                                                                                            gctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCgg                                                   <  2:1196713/150‑1 (MQ=255)
                                                                                               gctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGt                                               >  1:276676/1‑151 (MQ=255)
                                                                                                  ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTc                                              <  1:609468/149‑1 (MQ=255)
                                                                                                    aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAgg                                          >  2:1193449/1‑151 (MQ=255)
                                                                                                     tGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGc                                         >  1:319279/1‑151 (MQ=255)
                                                                                                          tGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTg                                    >  1:160889/1‑151 (MQ=255)
                                                                                                           gTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGc                                   <  1:4906/151‑1 (MQ=255)
                                                                                                           gTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGc                                   <  1:737516/151‑1 (MQ=255)
                                                                                                             cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTa                                 <  1:1015357/151‑1 (MQ=255)
                                                                                                              aCTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAg                                <  2:1208828/151‑1 (MQ=255)
                                                                                                               cTGGCCGGCATATGTTGAAGGTTAACGTTTATCGC‑C‑CCCCCCAAGCGGCATCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCGAGCTAGGCGTACTCAAGGAGCCGGAGGTCGAGGCTCCTGCTAg                                <  2:548752/150‑1 (MQ=255)
                                                                                                                       cATATTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAActctc                       <  2:799955/151‑1 (MQ=255)
                                                                                                                        atatTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAActctct                      >  2:603634/1‑151 (MQ=255)
                                                                                                                         tatTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGCGCCGGAGGTCGAGGCACCTGCTAGCCAActctct                      <  2:385416/150‑1 (MQ=255)
                                                                                                                          atTTTGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCtt                     <  1:597571/150‑1 (MQ=255)
                                                                                                                            tttGCAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCttt                    >  2:450998/1‑149 (MQ=255)
                                                                                                                                cAGGTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTgg              >  1:519076/1‑151 (MQ=255)
                                                                                                                                 aGGTTAACGTTGATCGCCC‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTgg              <  2:274443/151‑1 (MQ=255)
                                                                                                                                   gTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGTTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTGAATGAGCCGGAGGTCGAGGCTCCTGCTAGCCACCTCTCTTTAGCTGGTcg           <  2:643475/151‑1 (MQ=255)
                                                                                                                                   gTTAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTcg           >  2:1128032/1‑151 (MQ=255)
                                                                                                                                    ttAACGTTGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTcgc          >  2:1133651/1‑151 (MQ=255)
                                                                                                                                           tGATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATCATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGCCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCtttt   >  1:1192523/1‑151 (MQ=255)
                                                                                                                                            gATCGC‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCttttt  <  1:52968/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑C‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCTTTTT  >  NC_000913/307903‑308192

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