Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 308,054 (C)7→6 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,0480C.96.2% 341.5 / 8.7 80coding (1285/2526 nt)ecpCECP production outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (33/43);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGG  >  NC_000913/307949‑308143
                                                                                                   |                                                                                               
acaGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGc                                                                                                 <  2:2234541/100‑1 (MQ=255)
 caGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcc                                                                                                >  1:51108/1‑100 (MQ=255)
   gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              <  1:1360845/100‑1 (MQ=255)
   gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              <  2:941194/100‑1 (MQ=255)
   gCGAACCAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              <  2:1797340/100‑1 (MQ=255)
    cGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccc                                                                                             <  1:42894/100‑1 (MQ=255)
     gAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAGCGTTGATCGcccccc                                                                                            >  1:2935217/1‑100 (MQ=255)
        cTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccaag                                                                                         >  2:372268/1‑97 (MQ=255)
         tAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccaagc                                                                                        <  2:396085/100‑5 (MQ=255)
         tAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccaagc                                                                                        >  2:209490/1‑96 (MQ=255)
          agaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCg                                                                                      >  1:2017923/3‑100 (MQ=255)
          aaaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCg                                                                                      >  1:1540615/1‑100 (MQ=255)
            aGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCgg                                                                                     >  1:2508821/1‑100 (MQ=255)
            aGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGc                                                                                    <  2:2025451/100‑1 (MQ=255)
              ccTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAg                                                                                  <  2:1280264/100‑1 (MQ=255)
                tCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATTAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCg                                                                                >  1:1421391/1‑100 (MQ=255)
                tCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCg                                                                                <  1:2611131/100‑1 (MQ=255)
                 cccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGt                                                                               <  1:2956026/100‑1 (MQ=255)
                  ccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTc                                                                              <  1:280714/100‑1 (MQ=255)
                   cGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCa                                                                             <  2:2310743/100‑1 (MQ=255)
                    ggTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAg                                                                            <  2:1125683/100‑1 (MQ=255)
                    ggTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGAGGCCGCGTCAg                                                                            <  2:1567273/100‑1 (MQ=255)
                       agagTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACg                                                                         <  1:2545995/100‑1 (MQ=255)
                        gagTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACgg                                                                        >  2:2340066/1‑100 (MQ=255)
                         agTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACggg                                                                       <  2:2105483/100‑1 (MQ=255)
                            ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGttt                                                                    >  2:178976/1‑100 (MQ=255)
                            ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGttt                                                                    <  2:749810/100‑1 (MQ=255)
                            ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGttt                                                                    >  2:292320/1‑100 (MQ=255)
                            ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGttt                                                                    >  2:292871/1‑100 (MQ=255)
                              cgcTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTca                                                                  <  1:2004329/100‑1 (MQ=255)
                              cgcTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTca                                                                  >  2:2004317/1‑100 (MQ=255)
                                 tGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccc                                                               <  2:920196/100‑1 (MQ=255)
                                 tGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccc                                                               >  2:2396803/1‑100 (MQ=255)
                                  gATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccca                                                              >  1:729921/1‑100 (MQ=255)
                                  gATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccca                                                              >  1:730789/1‑100 (MQ=255)
                                  gATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccca                                                              >  2:789899/1‑100 (MQ=255)
                                    tgctgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTcacccacc                                                            <  1:1617228/100‑1 (MQ=255)
                                       tgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGcc                                                         <  1:2701466/100‑1 (MQ=255)
                                         ccGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                       <  1:2513771/100‑1 (MQ=255)
                                         ccGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                       <  2:2728847/100‑1 (MQ=255)
                                           gttgctgcTCTATGAGCTGTCACTTGCCAGCATAGTTTGCACGTTAACGTTGATCACCCCCTAAGCGTCTGCGTCAGACGGGTTTCACCCACCGCCTGAt                                                     <  2:2851885/98‑1 (MQ=255)
                                           gATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGAt                                                     <  1:2276143/100‑1 (MQ=255)
                                               ctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATc                                                 <  1:1644271/100‑1 (MQ=255)
                                               ctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATc                                                 >  2:475578/1‑100 (MQ=255)
                                                tgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCg                                                <  1:2178205/100‑1 (MQ=255)
                                                 gctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTAATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCg                                                <  2:1361709/99‑1 (MQ=255)
                                                  ctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTa                                              <  2:794226/100‑1 (MQ=255)
                                                   tCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTAt                                             <  1:2340080/100‑1 (MQ=255)
                                                    ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            >  2:822874/1‑100 (MQ=255)
                                                    ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            >  2:2638780/1‑100 (MQ=255)
                                                    ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            >  1:2305943/1‑100 (MQ=255)
                                                    ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            >  1:1853984/1‑100 (MQ=255)
                                                    ccATGACCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            <  2:2921959/100‑1 (MQ=255)
                                                     cATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATca                                           >  2:213384/1‑99 (MQ=255)
                                                     cATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            <  1:213386/99‑1 (MQ=255)
                                                      aTGAGCTGTCAGTGGGCAGTAGATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATcc                                           <  2:2508806/100‑1 (MQ=255)
                                                      aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCa                                          <  1:1744612/100‑1 (MQ=255)
                                                        gAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCata                                        >  1:1957388/1‑100 (MQ=255)
                                                        gAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCata                                        >  2:2479176/1‑100 (MQ=255)
                                                         aGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCatat                                       >  2:2455647/1‑100 (MQ=255)
                                                           cTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATcc                                     <  2:2526786/100‑1 (MQ=255)
                                                            tGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCatat                                       >  2:136292/1‑97 (MQ=255)
                                                            tGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATccc                                    >  1:31265/1‑100 (MQ=255)
                                                               cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                 >  1:1583278/1‑100 (MQ=255)
                                                                   ggCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccg                             >  1:1893718/1‑100 (MQ=255)
                                                                   ggCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccg                             >  2:340457/1‑100 (MQ=255)
                                                                      cAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGccc                          >  2:2108037/1‑100 (MQ=255)
                                                                         cATATTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGc                       <  2:1542143/100‑1 (MQ=255)
                                                                           tatTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                     <  1:2248510/100‑1 (MQ=255)
                                                                           tatTTTGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                     <  2:1549502/100‑1 (MQ=255)
                                                                             ttttGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAg                   <  1:965848/100‑1 (MQ=255)
                                                                             ttttGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAg                   <  2:36701/100‑1 (MQ=255)
                                                                               ttGCAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCg                 <  1:2459245/100‑1 (MQ=255)
                                                                                  cAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTAc              <  1:209492/100‑1 (MQ=255)
                                                                                  cAGGTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTAc              >  2:742438/1‑100 (MQ=255)
                                                                                    ggTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTc            <  1:91712/100‑1 (MQ=255)
                                                                                    ggTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTc            <  1:911073/100‑1 (MQ=255)
                                                                                    ggTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTc            >  1:2777741/1‑100 (MQ=255)
                                                                                     gTTAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCa           <  2:2017911/100‑1 (MQ=255)
                                                                                       tAACGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAAt         <  1:475581/100‑1 (MQ=255)
                                                                                        aaCGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATg        <  2:729918/100‑1 (MQ=255)
                                                                                        aaCGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATg        <  2:730786/100‑1 (MQ=255)
                                                                                         aCGTTGATCGCCCCCCAAGCGGCAGCGTCAGGCGGGTCTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa       <  2:2014812/100‑1 (MQ=255)
                                                                                         aCGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa       >  1:298511/1‑100 (MQ=255)
                                                                                          cGTTGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAg      >  2:128811/1‑100 (MQ=255)
                                                                                            ttGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGcc    <  2:351578/100‑1 (MQ=255)
                                                                                            ttGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGcc    <  1:777366/100‑1 (MQ=255)
                                                                                            ttGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGcc    >  1:1336787/1‑100 (MQ=255)
                                                                                            ttGATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGcc    <  2:1699944/100‑1 (MQ=255)
                                                                                              gATCGCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCgg  >  2:974612/1‑100 (MQ=255)
                                                                                                   |                                                                                               
ACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGG  >  NC_000913/307949‑308143

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