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‑‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCTTTTT  >  NC_000913/307903‑308192
                                                                                                                                                 |                                                                                                                                                   
cTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑cccccc                                                                                                                                             >  1:673991/1‑151 (MQ=255)
 tACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑ccccccc                                                                                                                                            >  1:315423/1‑151 (MQ=255)
   cGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑cccccccc                                                                                                                                           >  2:738078/1‑149 (MQ=255)
     tCGCCATTGATTGTAAATGCGGGTCTTTTCCTCACAAACCACCAGCGTCTTAATGCTTCCCGGTCGCGTGGCGCTGATTATGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGAATCATTTGAACGATAACGTCGATCggcacctcga                                                                                                                                             >  2:207630/1‑140 (MQ=255)
      cGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑ccccccccaagc                                                                                                                                       <  1:1036899/151‑6 (MQ=255)
      cGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑ccccccccaag                                                                                                                                        <  1:675029/150‑5 (MQ=255)
      cGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑CCCCCCCAAGCg                                                                                                                                       <  2:542522/151‑1 (MQ=255)
        cAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCgg                                                                                                                                      <  1:597024/151‑1 (MQ=255)
         aaTTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGc                                                                                                                                     >  2:489908/1‑151 (MQ=255)
          attgattgCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCa                                                                                                                                    <  1:1080309/151‑1 (MQ=255)
             gattgCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCg                                                                                                                                 >  1:966665/1‑151 (MQ=255)
                 gCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTc                                                                                                                               >  2:776671/1‑149 (MQ=255)
                 gCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAg                                                                                                                             <  2:1141264/151‑1 (MQ=255)
                       gCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACggg                                                                                                                        <  1:125399/150‑1 (MQ=255)
                         gggTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGttt                                                                                                                     >  1:852300/1‑151 (MQ=255)
                         gggTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGttt                                                                                                                     >  2:872452/1‑151 (MQ=255)
                          ggTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGttt                                                                                                                     >  2:926585/1‑150 (MQ=255)
                                   cTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGc                                                                                                           >  2:737511/1‑151 (MQ=255)
                                     gATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGc                                                                                                           <  1:764546/149‑1 (MQ=255)
                                      aTTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                                                                        >  1:814384/1‑151 (MQ=255)
                                         aaCCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGAtt                                                                                                     <  1:964866/151‑1 (MQ=255)
                                         aaCCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGAtt                                                                                                     >  1:1216809/1‑151 (MQ=255)
                                          aCCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGAtt                                                                                                     >  1:143517/1‑150 (MQ=255)
                                           cccACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGTGTCAGACGGGTTTCACCCACCGCCTGATTAt                                                                                                   <  1:1129003/151‑1 (MQ=255)
                                           cccACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTAt                                                                                                   >  1:285193/1‑151 (MQ=255)
                                             cacaGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCg                                                                                                 <  2:690026/151‑1 (MQ=255)
                                                aGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTAt                                                                                              >  1:2997/1‑151 (MQ=255)
                                                    aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCat                                                                                          >  1:1218951/1‑151 (MQ=255)
                                                    aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCat                                                                                          <  1:634277/151‑1 (MQ=255)
                                                    aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCat                                                                                          >  1:804812/1‑151 (MQ=255)
                                                        aaaGCGTCCCGGTAGAGTGGCGCTGATGCTGCCGATAATGCTGCGTGAGCTGTCGCTGGCCCTCATATTTTGGAGGTTAACGTTGATCT‑‑CGCCCCCCAGGCGCCAGCGCCTGACGGGCCTCACCCACCATCTGATTATCGCATCCATAGcc                                                                                      >  2:525632/1‑151 (MQ=255)
                                                        aaaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCatat                                                                                        <  1:924944/149‑1 (MQ=255)
                                                         aaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGTGGTGTCACTGGCCGGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCatat                                                                                        <  2:322984/148‑1 (MQ=255)
                                                          aGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCg                                                                                    <  2:633007/151‑1 (MQ=255)
                                                          aGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGt                                                                                   <  1:50917/151‑1 (MQ=255)
                                                            ccTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                                                                  >  1:661807/1‑151 (MQ=255)
                                                                ccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccg                                                                              >  1:623078/1‑151 (MQ=255)
                                                                     agagTGGCGCTGATGCTGCCGATGCTGCTCCATGAGTTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCa                                                                          <  2:230908/150‑1 (MQ=255)
                                                                      gagTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGc                                                                        <  1:1098682/151‑1 (MQ=255)
                                                                       agTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCt                                                                       <  2:1208196/151‑1 (MQ=255)
                                                                        gTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                                                                      <  1:692351/151‑1 (MQ=255)
                                                                         tGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCTGCCCAGCTa                                                                      >  2:1210662/1‑150 (MQ=255)
                                                                         tGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                                                                      >  2:806446/1‑150 (MQ=255)
                                                                         tGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                                                                      >  2:895281/1‑150 (MQ=255)
                                                                         tGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATAGTATCCATATCCCGTTGCCGCCCAGCTa                                                                      >  2:803458/1‑150 (MQ=255)
                                                                               tGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTAc                                                               >  1:29135/1‑151 (MQ=255)
                                                                                    ctgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCaa                                                           <  1:293599/150‑1 (MQ=255)
                                                                                      gcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa                                                        <  1:168122/151‑1 (MQ=255)
                                                                                      gcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa                                                        <  1:51397/151‑1 (MQ=255)
                                                                                      gcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa                                                        <  1:388333/151‑1 (MQ=255)
                                                                                        cGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc                                                      >  2:1142725/1‑151 (MQ=255)
                                                                                         gATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGcc                                                     >  2:500751/1‑151 (MQ=255)
                                                                                          aTGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCg                                                    <  2:990541/151‑1 (MQ=255)
                                                                                            gctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCgg                                                   <  2:1196713/150‑1 (MQ=255)
                                                                                               gctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGt                                               >  1:276676/1‑151 (MQ=255)
                                                                                                  ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTc                                              <  1:609468/149‑1 (MQ=255)
                                                                                                    aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAgg                                          >  2:1193449/1‑151 (MQ=255)
                                                                                                     tGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGc                                         >  1:319279/1‑151 (MQ=255)
                                                                                                          tGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTg                                    >  1:160889/1‑151 (MQ=255)
                                                                                                           gTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGc                                   <  1:4906/151‑1 (MQ=255)
                                                                                                           gTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGc                                   <  1:737516/151‑1 (MQ=255)
                                                                                                             cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTa                                 <  1:1015357/151‑1 (MQ=255)
                                                                                                              aCTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAg                                <  2:1208828/151‑1 (MQ=255)
                                                                                                               cTGGCCGGCATATGTTGAAGGTTAACGTTTATCGCCCCCCCCAAGCGGCATCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCGAGCTAGGCGTACTCAAGGAGCCGGAGGTCGAGGCTCCTGCTAg                                <  2:548752/150‑1 (MQ=255)
                                                                                                                       cATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAActctc                       <  2:799955/151‑1 (MQ=255)
                                                                                                                        atatTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAActctct                      >  2:603634/1‑151 (MQ=255)
                                                                                                                         tatTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGCGCCGGAGGTCGAGGCACCTGCTAGCCAActctct                      <  2:385416/150‑1 (MQ=255)
                                                                                                                          atTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCtt                     <  1:597571/150‑1 (MQ=255)
                                                                                                                            tttGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCttt                    >  2:450998/1‑149 (MQ=255)
                                                                                                                                cAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTgg              >  1:519076/1‑151 (MQ=255)
                                                                                                                                 aGGTTAACGTTGATCGCCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTgg              <  2:274443/151‑1 (MQ=255)
                                                                                                                                   gTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGTTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTGAATGAGCCGGAGGTCGAGGCTCCTGCTAGCCACCTCTCTTTAGCTGGTcg           <  2:643475/151‑1 (MQ=255)
                                                                                                                                   gTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTcg           >  2:1128032/1‑151 (MQ=255)
                                                                                                                                    ttAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTcgc          >  2:1133651/1‑151 (MQ=255)
                                                                                                                                           tGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATCATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGCCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCtttt   >  1:1192523/1‑151 (MQ=255)
                                                                                                                                            gATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCttttt  <  1:52968/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCTTTTT  >  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: