BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,193	(C)_6→7	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	94.2%	408.1 / 17.7	120	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (4/3); new base C (56/57); total (60/60)
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.96e-01

Only 100 of 135 total aligned reads displayed.

CTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638090‑4638281

 tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cc                                                                                             <  1:2791070/100‑1 (MQ=255)

   gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccc                                                                                           <  2:2643715/100‑1 (MQ=255)

   gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccc                                                                                           <  1:2750892/100‑1 (MQ=255)

   gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccc                                                                                           >  1:2377377/1‑100 (MQ=255)

    cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑ccccc                                                                                          >  1:1786553/1‑100 (MQ=255)

     ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccc                                                                                         <  1:2183708/100‑2 (MQ=255)

     ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccc                                                                                         >  2:379172/1‑99 (MQ=255)

      tCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccct                                                                                        >  2:292925/1‑98 (MQ=255)

        gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccctgt                                                                                      >  2:277608/1‑96 (MQ=255)

        gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccctgt                                                                                      >  1:2127795/1‑96 (MQ=255)

        gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccctgt                                                                                      >  2:599279/1‑96 (MQ=255)

        gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CTCCCCCTGt                                                                                       >  2:1347944/1‑100 (MQ=255)

         aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccctgtt                                                                                     >  1:904836/1‑95 (MQ=255)

         aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccctgtt                                                                                     >  2:781745/1‑95 (MQ=255)

          cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTc                                                                                     >  1:850440/1‑100 (MQ=255)

          cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTc                                                                                     >  1:92749/1‑100 (MQ=255)

          cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccctgtt                                                                                     >  2:2901951/1‑94 (MQ=255)

           ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCt                                                                                    >  1:762423/1‑100 (MQ=255)

             cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTgg                                                                                  >  2:2951098/1‑100 (MQ=255)

             cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CCCCCTGTTCTGGa                                                                                 >  1:749495/1‑100 (MQ=255)

               ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGaa                                                                                >  2:1585257/1‑100 (MQ=255)

                  tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAAt                                                                             <  2:1038319/100‑1 (MQ=255)

                  tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CCCCCTGTTCTGGAAAATg                                                                            <  2:706894/100‑1 (MQ=255)

                   cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATg                                                                            <  2:2030032/100‑1 (MQ=255)

                   cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATg                                                                            <  2:2853358/100‑1 (MQ=255)

                    ncTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGa                                                                           >  2:3039607/2‑100 (MQ=255)

                    aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGa                                                                           >  1:2611194/1‑100 (MQ=255)

                     cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGAc                                                                          >  1:3076933/1‑100 (MQ=255)

                        cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAg                                                                       >  2:2691922/1‑100 (MQ=255)

                          caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGcc                                                                     <  2:2120129/100‑1 (MQ=255)

                          caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGcc                                                                     <  2:2611182/100‑1 (MQ=255)

                              gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:517837/1‑100 (MQ=255)

                              gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:3030322/1‑100 (MQ=255)

                               cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  2:2489675/1‑100 (MQ=255)

                               cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                <  1:2484130/100‑1 (MQ=255)

                                ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  1:961579/100‑1 (MQ=255)

                                ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  1:446650/100‑1 (MQ=255)

                                ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:620170/1‑100 (MQ=255)

                                ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:61390/1‑100 (MQ=255)

                                ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:132915/1‑100 (MQ=255)

                                 tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  2:2010644/1‑100 (MQ=255)

                                 tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  2:886779/100‑1 (MQ=255)

                                   aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            <  1:2166415/100‑1 (MQ=255)

                                    aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:1350164/100‑1 (MQ=255)

                                       tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                        <  1:1684511/100‑1 (MQ=255)

                                         ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  2:909540/1‑100 (MQ=255)

                                          tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAtt                                                     >  1:297228/1‑100 (MQ=255)

                                          tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAtt                                                     <  2:1438119/100‑1 (MQ=255)

                                          tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAtt                                                     >  2:2901357/1‑100 (MQ=255)

                                           cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  2:2671440/100‑1 (MQ=255)

                                            ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  1:789892/1‑100 (MQ=255)

                                            ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  1:2903219/100‑1 (MQ=255)

                                            ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:3148376/100‑1 (MQ=255)

                                            ccACATAGTCATCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGCAAATGACTAGCCTGTTTGGTGCAGTATTgn                                                   <  2:2201587/100‑2 (MQ=255)

                                             cacaTAGTCTTCGTATCCTGCTGCCATTGCAATGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  <  2:2352605/100‑1 (MQ=255)

                                             cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  1:3080145/1‑100 (MQ=255)

                                             cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:1557747/1‑100 (MQ=255)

                                             cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  <  1:1166202/100‑1 (MQ=255)

                                              acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  2:950806/100‑1 (MQ=255)

                                                 tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  1:2412099/100‑1 (MQ=255)

                                                  aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             <  1:1499513/100‑1 (MQ=255)

                                                   gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  2:1032521/1‑100 (MQ=255)

                                                    tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:61391/100‑1 (MQ=255)

                                                     cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  1:1030404/1‑100 (MQ=255)

                                                      ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  2:2344516/100‑1 (MQ=255)

                                                      ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  2:2962082/1‑100 (MQ=255)

                                                      ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  1:856643/1‑100 (MQ=255)

                                                        cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGTAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       <  2:2772512/100‑1 (MQ=255)

                                                         gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:951504/100‑1 (MQ=255)

                                                          tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  1:1485498/1‑100 (MQ=255)

                                                            tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  1:1976734/100‑1 (MQ=255)

                                                            tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  2:3101034/100‑1 (MQ=255)

                                                            tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CTCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  1:1347948/100‑1 (MQ=255)

                                                             ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  <  2:1894199/100‑1 (MQ=255)

                                                               tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                >  2:1589630/1‑100 (MQ=255)

                                                                  tgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             >  2:242773/1‑100 (MQ=255)

                                                                  tgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:1440207/100‑1 (MQ=255)

                                                                   gcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            <  1:2901966/100‑1 (MQ=255)

                                                                   gcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  1:2622434/1‑100 (MQ=255)

                                                                        tGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       >  1:414802/1‑100 (MQ=255)

                                                                        tGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       <  2:878831/100‑1 (MQ=255)

                                                                          cAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:2220530/1‑100 (MQ=255)

                                                                           aaaGGAGAAGACTTTGTTGAAATCC‑CCCCCTGTTCTTGAAAATGACTGGCCTGTTTGGTGCAGTATTGCTTTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:2561850/100‑1 (MQ=255)

                                                                           aaaGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  1:1692705/100‑1 (MQ=255)

                                                                           aaaGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:634169/100‑1 (MQ=255)

                                                                              ggAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  2:1276288/1‑100 (MQ=255)

                                                                              ggAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATACTGATTCGGCAGGt                 <  2:2127785/100‑1 (MQ=255)

                                                                               gagaAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                <  1:3019001/100‑1 (MQ=255)

                                                                               gagaAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                <  1:3095809/100‑1 (MQ=255)

                                                                               gagaAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                >  2:2166072/1‑100 (MQ=255)

                                                                                 gaagaCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              <  1:2901372/100‑1 (MQ=255)

                                                                                   agaCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  2:458505/1‑100 (MQ=255)

                                                                                        aTGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       <  1:2010652/100‑1 (MQ=255)

                                                                                         tGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      <  2:803511/100‑1 (MQ=255)

                                                                                          gTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:1025337/1‑100 (MQ=255)

                                                                                          gTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     <  1:277609/100‑1 (MQ=255)

                                                                                           ttGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAGGt    >  2:2610326/1‑100 (MQ=255)

                                                                                            tGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:2000362/1‑99 (MQ=255)

                                                                                            tGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   <  2:3319066/100‑1 (MQ=255)

                                                                                            tGAAAT‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  2:85107/1‑100 (MQ=255)

CTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638090‑4638281

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: `‑`