Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 308,054 (C)7→8 coding (1279/2526 nt) ecpC ← putative fimbrial usher protein EcpC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913308,0471.C98.3% 974.0 / 6.7 233coding (1286/2526 nt)ecpCputative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (123/106);  total (125/108)
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.90e-01

Only 100 of 278 total aligned reads displayed.

CCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCG  >  NC_000913/307947‑308142
                                                                                                     |                                                                                                
ccACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCg                                                                                                   >  2:567775/1‑101 (MQ=255)
  acaGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑cc                                                                                               >  2:567780/1‑101 (MQ=255)
  acaGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACATTGATCG‑‑cc                                                                                               <  1:567180/101‑1 (MQ=255)
   caGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑ccc                                                                                              <  2:567529/101‑1 (MQ=255)
   caGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑ccc                                                                                              <  1:567387/101‑1 (MQ=255)
    aGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑cccc                                                                                             <  1:567020/101‑1 (MQ=255)
      cGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑cccccc                                                                                           >  1:567783/1‑101 (MQ=255)
       gAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTACCGTTGATCG‑‑ccccccc                                                                                          <  2:567196/101‑1 (MQ=255)
        aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑cccccccc                                                                                         <  1:567047/101‑2 (MQ=255)
        aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑cccccccc                                                                                         >  2:567787/1‑100 (MQ=255)
        aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑cccccccc                                                                                         >  2:567785/1‑100 (MQ=255)
           tAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑ccccccccaag                                                                                      >  2:567792/1‑97 (MQ=255)
           tAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑ccccccccaag                                                                                      >  1:567791/1‑97 (MQ=255)
            aaaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑ccccccccaagc                                                                                     >  2:567798/1‑96 (MQ=255)
            aaaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑ccccccccaagc                                                                                     <  1:567432/101‑6 (MQ=255)
            aaaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑ccccccccaagc                                                                                     <  2:567293/101‑6 (MQ=255)
             aaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCg                                                                                     <  2:567385/101‑1 (MQ=255)
               gCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGc                                                                                   <  1:567339/101‑1 (MQ=255)
               gCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGc                                                                                   >  1:567800/1‑101 (MQ=255)
                 cTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAg                                                                                 >  2:567803/1‑101 (MQ=255)
                  tCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCCAAGCGGCAg                                                                                 <  2:567417/101‑1 (MQ=255)
                  tCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGc                                                                                <  2:567070/101‑1 (MQ=255)
                  tCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGc                                                                                >  1:567805/1‑101 (MQ=255)
                   cccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCg                                                                               <  1:567369/101‑1 (MQ=255)
                    ccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCg                                                                               >  2:567808/1‑100 (MQ=255)
                      ggTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCa                                                                            >  1:567810/1‑101 (MQ=255)
                       gTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAg                                                                           >  1:567814/1‑101 (MQ=255)
                        tAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGa                                                                          >  2:567816/1‑101 (MQ=255)
                         agagTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGAc                                                                         <  2:567078/101‑1 (MQ=255)
                           agTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACgg                                                                       >  1:567822/1‑101 (MQ=255)
                            gTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACggg                                                                      <  1:567519/101‑1 (MQ=255)
                              ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGtt                                                                    >  1:567824/1‑101 (MQ=255)
                              ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGtt                                                                    <  2:567286/101‑1 (MQ=255)
                              ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGtt                                                                    <  1:567153/101‑1 (MQ=255)
                                cgcTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTc                                                                  <  1:567729/101‑1 (MQ=255)
                                cgcTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTc                                                                  >  2:567826/1‑101 (MQ=255)
                                  cTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTcac                                                                >  1:567829/1‑101 (MQ=255)
                                    gATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccc                                                              >  1:567832/1‑101 (MQ=255)
                                       gctgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTcacccacc                                                           <  1:567227/101‑1 (MQ=255)
                                       gctgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTcacccacc                                                           <  1:567412/101‑1 (MQ=255)
                                         tgcCGATGCTGCTCCATGAGCTGTCGCTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGc                                                         >  1:567836/1‑101 (MQ=255)
                                           ccGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCt                                                       >  2:567838/1‑101 (MQ=255)
                                            cGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                      <  1:567631/101‑1 (MQ=255)
                                                 ctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTAt                                                 <  2:567563/101‑1 (MQ=255)
                                                  tgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATc                                                >  1:567841/1‑101 (MQ=255)
                                                    ctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGt                                              >  1:567845/1‑101 (MQ=255)
                                                    ctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGt                                              >  2:567844/1‑101 (MQ=255)
                                                      ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTAt                                            >  2:567848/1‑101 (MQ=255)
                                                        aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATcc                                          >  2:567852/1‑101 (MQ=255)
                                                          gAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCat                                        >  1:567853/1‑101 (MQ=255)
                                                            gCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCatat                                      >  1:567856/1‑101 (MQ=255)
                                                             cTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATc                                     <  1:567624/101‑1 (MQ=255)
                                                              tGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATcc                                    <  2:567487/101‑1 (MQ=255)
                                                               gTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATccc                                   >  2:567861/1‑101 (MQ=255)
                                                               gTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATccc                                   <  2:567258/101‑1 (MQ=255)
                                                               gTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATccc                                   <  1:567317/101‑1 (MQ=255)
                                                                tCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCg                                  <  1:567573/101‑1 (MQ=255)
                                                                tCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCg                                  <  2:567344/101‑1 (MQ=255)
                                                                tCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCg                                  >  1:567863/1‑101 (MQ=255)
                                                                  aCTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                >  1:567865/1‑101 (MQ=255)
                                                                   cTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                <  2:567444/101‑1 (MQ=255)
                                                                   cTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTg                               <  2:567006/101‑1 (MQ=255)
                                                                   cTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTg                               <  2:567604/101‑1 (MQ=255)
                                                                    tGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGcc                                                        >  2:567870/1‑75 (MQ=255)
                                                                    tGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑CCCCCCCAAGCGGCAg                                                                                 >  1:567872/1‑49 (MQ=255)
                                                                     ggCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgcc                             >  2:567874/1‑101 (MQ=255)
                                                                      gCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccg                            >  1:567876/1‑101 (MQ=255)
                                                                         aGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccgcc                          >  1:567883/1‑101 (MQ=255)
                                                                         aGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGccc                         >  2:567881/1‑101 (MQ=255)
                                                                          gCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCa                        <  2:567392/101‑1 (MQ=255)
                                                                           cATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAg                       >  2:567886/1‑101 (MQ=255)
                                                                            atatTTTGCAGGTTAACGTTGATCGCCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCata                                       <  2:567362/85‑1 (MQ=255)
                                                                            atatTTTGCAGGTTAACGTTGATCGCCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAg                       <  2:567517/101‑1 (MQ=255)
                                                                            atatTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGc                      <  2:567161/101‑1 (MQ=255)
                                                                            atatTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGc                      >  1:567888/1‑101 (MQ=255)
                                                                             tatTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCa                        <  2:567381/98‑1 (MQ=255)
                                                                             tatTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCt                     <  1:567499/101‑1 (MQ=255)
                                                                              atTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTa                    >  1:567892/1‑101 (MQ=255)
                                                                               ttttGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTaa                   >  2:567894/1‑101 (MQ=255)
                                                                                tttGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAg                  <  1:567469/101‑1 (MQ=255)
                                                                                 ttGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGc                 >  2:567898/1‑101 (MQ=255)
                                                                                 ttGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGc                 <  1:567448/101‑1 (MQ=255)
                                                                                 ttGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGc                 <  1:567256/101‑1 (MQ=255)
                                                                                    cAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTa              >  1:567900/1‑101 (MQ=255)
                                                                                    cAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTa              >  2:567902/1‑101 (MQ=255)
                                                                                      ggTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCAAAGCGTAca            <  2:567552/101‑2 (MQ=255)
                                                                                        ttAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCa          >  2:567906/1‑101 (MQ=255)
                                                                                         tAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCaa         >  1:567908/1‑101 (MQ=255)
                                                                                         tAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCaa         <  1:567539/101‑1 (MQ=255)
                                                                                          aaCGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAGGCGTACTCAAt        <  1:567597/101‑1 (MQ=255)
                                                                                          aaCGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAAt        >  1:567912/1‑101 (MQ=255)
                                                                                           aCGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATg       >  2:567915/1‑101 (MQ=255)
                                                                                           aCGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATg       >  1:567916/1‑101 (MQ=255)
                                                                                             gTTGATCGCCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa      >  2:567921/1‑101 (MQ=255)
                                                                                             gTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAg     >  1:567918/1‑101 (MQ=255)
                                                                                             gTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAg     <  1:567484/101‑1 (MQ=255)
                                                                                             gTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAg     <  1:567296/101‑1 (MQ=255)
                                                                                              ttGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc    <  1:567456/101‑1 (MQ=255)
                                                                                               tGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgcc                             >  1:567925/1‑75 (MQ=255)
                                                                                                gATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCg  <  2:567464/101‑1 (MQ=255)
                                                                                                     |                                                                                                
CCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCG  >  NC_000913/307947‑308142

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