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	98.1%	366.6 / ‑0.6	104	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (53/49); total (53/51)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638089‑4638281

ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                                >  2:2839235/1‑100 (MQ=255)

   cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑cc                                                                                            >  2:567312/1‑100 (MQ=255)

   cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑cc                                                                                            <  1:2298084/100‑1 (MQ=255)

    gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccc                                                                                           >  2:638395/1‑100 (MQ=255)

    gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccc                                                                                           >  1:211088/1‑100 (MQ=255)

     cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑cccc                                                                                          >  1:2802581/1‑100 (MQ=255)

      ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccccc                                                                                         <  1:263379/100‑2 (MQ=255)

      ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccccc                                                                                         <  1:2294034/100‑2 (MQ=255)

       tCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccccct                                                                                        >  1:480727/1‑98 (MQ=255)

        cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccccctg                                                                                       <  1:2758934/100‑4 (MQ=255)

         gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccccctgt                                                                                      <  1:3261046/100‑5 (MQ=255)

         gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccccctgt                                                                                      >  1:2064533/1‑96 (MQ=255)

           cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTc                                                                                     >  2:87334/1‑100 (MQ=255)

           cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTc                                                                                     >  1:1644027/1‑100 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTg                                                                                   <  2:1417520/100‑1 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTgn                                                                                  <  2:2410658/100‑2 (MQ=255)

                ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATTTTGAAATCCC‑CCCCTGTTCTGGaa                                                                                >  2:636849/1‑100 (MQ=255)

                ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGaa                                                                                >  2:1127429/1‑100 (MQ=255)

                ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGaa                                                                                <  2:1581006/100‑1 (MQ=255)

                   tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAAt                                                                             <  2:371414/100‑1 (MQ=255)

                     aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGa                                                                           >  2:3198626/1‑100 (MQ=255)

                      cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGAc                                                                          <  2:2564299/100‑1 (MQ=255)

                         cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAg                                                                       <  2:2176827/100‑1 (MQ=255)

                          acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGc                                                                      <  1:680559/100‑1 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGcc                                                                     <  1:1127443/100‑1 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGcc                                                                     <  2:2463987/100‑1 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGcc                                                                     <  1:638404/100‑1 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGcc                                                                     <  1:819822/100‑1 (MQ=255)

                            aTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCt                                                                    <  1:2119409/100‑1 (MQ=255)

                             tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTg                                                                   <  1:2649068/100‑1 (MQ=255)

                             tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTg                                                                   <  2:2648427/100‑1 (MQ=255)

                               gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  2:632273/1‑100 (MQ=255)

                                cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  2:2638971/1‑100 (MQ=255)

                                cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  2:2221398/1‑100 (MQ=255)

                                 ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:2246073/1‑100 (MQ=255)

                                  tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  1:903368/1‑100 (MQ=255)

                                   cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGATTgg                                                              >  1:1172501/1‑99 (MQ=255)

                                    aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  1:1146614/1‑100 (MQ=255)

                                    aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  1:2604727/1‑100 (MQ=255)

                                     aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:2594802/100‑1 (MQ=255)

                                     aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:3220725/100‑1 (MQ=255)

                                            cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  1:1832571/100‑1 (MQ=255)

                                            cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  2:1342797/100‑1 (MQ=255)

                                            cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  2:2082494/100‑1 (MQ=255)

                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  1:223000/1‑100 (MQ=255)

                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGAATATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  2:792010/1‑100 (MQ=255)

                                               acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGcag                                                 >  1:2319434/1‑98 (MQ=255)

                                               acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 >  1:2742632/1‑100 (MQ=255)

                                                 aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  1:2976468/100‑1 (MQ=255)

                                                  tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGGCTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  2:596629/100‑1 (MQ=255)

                                                   aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             >  2:152270/1‑100 (MQ=255)

                                                    gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:1369215/100‑1 (MQ=255)

                                                    gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:2856780/100‑1 (MQ=255)

                                                    gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:1995144/100‑1 (MQ=255)

                                                     tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:464377/100‑1 (MQ=255)

                                                     tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:1647673/100‑1 (MQ=255)

                                                      nttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATCGCTGTTGTTGa                                          >  2:2262521/2‑100 (MQ=255)

                                                      cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  1:2549411/1‑100 (MQ=255)

                                                      cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  2:2329974/1‑100 (MQ=255)

                                                       ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  1:2346501/1‑100 (MQ=255)

                                                       ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  2:1257526/1‑100 (MQ=255)

                                                       ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  2:695880/100‑1 (MQ=255)

                                                        tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTCTTGCTGTTGTTGAtt                                        >  1:3212874/1‑100 (MQ=255)

                                                        tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        <  2:857788/100‑1 (MQ=255)

                                                          gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  2:3231845/100‑1 (MQ=255)

                                                          gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  2:167523/100‑1 (MQ=255)

                                                           tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     <  2:1045147/100‑1 (MQ=255)

                                                            aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    >  2:2641748/1‑100 (MQ=255)

                                                             tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  2:1850842/1‑100 (MQ=255)

                                                                   tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  2:1699836/100‑1 (MQ=255)

                                                                    gcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:200904/1‑100 (MQ=255)

                                                                     ccATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           <  1:2418462/100‑1 (MQ=255)

                                                                       aTTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         <  1:567352/100‑1 (MQ=255)

                                                                        ttGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        <  1:2018365/100‑1 (MQ=255)

                                                                         tGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       >  1:3117345/1‑100 (MQ=255)

                                                                         tGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       >  2:1086197/1‑100 (MQ=255)

                                                                         tGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       >  1:983259/1‑100 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:1923494/1‑100 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  1:304789/100‑1 (MQ=255)

                                                                           cAAAGGAAAAGACTATGTTGAAAT‑CCGCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCGGTTGTTAATTCCGATAATGCTGATTCGGc                     <  2:1532483/100‑1 (MQ=255)

                                                                            aaaGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    >  1:1614941/1‑100 (MQ=255)

                                                                               ggAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  2:1188802/1‑100 (MQ=255)

                                                                                gagaAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                <  2:3180507/100‑1 (MQ=255)

                                                                                gagaAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                >  1:3150919/1‑100 (MQ=255)

                                                                                 agaagaCTATGTTGAAGTCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTGCGATAATGCTGATTCGGCAgg                  >  1:784269/1‑97 (MQ=255)

                                                                                 agaagaCTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGa               >  2:2439966/1‑100 (MQ=255)

                                                                                 agaagaCTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGa               >  1:935525/1‑100 (MQ=255)

                                                                                   aagaCTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAtt             <  1:1602853/100‑1 (MQ=255)

                                                                                     gaCTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           <  1:1086211/100‑1 (MQ=255)

                                                                                       cTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  1:1017324/1‑100 (MQ=255)

                                                                                         aTGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       >  2:2820567/1‑100 (MQ=255)

                                                                                           gTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:251165/1‑100 (MQ=255)

                                                                                           gTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:2184708/1‑100 (MQ=255)

                                                                                            ttGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  2:500260/1‑100 (MQ=255)

                                                                                            ttGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  2:1175489/1‑100 (MQ=255)

                                                                                            ttGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:2389836/1‑100 (MQ=255)

                                                                                             tGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGGACGTg   >  1:646406/1‑100 (MQ=255)

                                                                                             tGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   <  2:2778416/100‑1 (MQ=255)

                                                                                             tGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   >  2:2314355/1‑100 (MQ=255)

                                                                                              gAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  2:986770/100‑1 (MQ=255)

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638089‑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: `‑`