BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,088,752	A→G	V110A (GTG→GCG)	ycdR ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,088,752	0	A	G	99.2%	420.2 / ‑5.0	130	V110A (GTG→GCG)	ycdR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (74/55); total (74/56)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 132 total aligned reads displayed.

ACGGGGGCCCATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGC  >  NC_000913/1088662‑1088841

aCGGGGGCCCATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCt                                                                                         <  2:2744909/93‑1 (MQ=255)

aCGGGGGCCCATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCt                                                                                         >  1:520442/1‑93 (MQ=255)

  gggggCCCATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTAc                                                                                       <  1:3907575/93‑1 (MQ=255)

   ggggCCCATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACa                                                                                      <  2:1288859/93‑1 (MQ=255)

   ggggCCCATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACa                                                                                      >  1:368783/1‑93 (MQ=255)

   ggggCCCATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACa                                                                                      >  2:3366741/1‑93 (MQ=255)

     ggCCCATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGc                                                                                    <  1:2456563/93‑1 (MQ=255)

      gCCCATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCt                                                                                   <  1:812999/93‑1 (MQ=255)

        ccATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCttt                                                                                 <  1:1138083/93‑1 (MQ=255)

        ccATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCttt                                                                                 >  1:3691088/1‑93 (MQ=255)

         cATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCtttt                                                                                <  2:1742207/93‑1 (MQ=255)

         cATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCtttt                                                                                <  2:3474603/93‑1 (MQ=255)

          aTACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCttttt                                                                               <  2:4009000/93‑1 (MQ=255)

          aTACAGCAGGCCACTGGAAGGCCTGAAGAAGTGGGAAGACGCGTGTATAAAACCTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCttttt                                                                               <  1:2859977/93‑1 (MQ=255)

           tACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTc                                                                              >  2:4301743/1‑93 (MQ=255)

           tACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTc                                                                              >  1:564609/1‑93 (MQ=255)

            ccagcagGCCACGGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTcc                                                                             <  2:3559430/92‑1 (MQ=255)

            aCAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTcc                                                                             <  2:3559431/93‑1 (MQ=255)

            aCAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTcc                                                                             >  1:4092836/1‑93 (MQ=255)

            aCAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTcc                                                                             >  2:4750100/1‑93 (MQ=255)

            aCAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTcc                                                                             <  2:2650934/93‑1 (MQ=255)

               gcagGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGt                                                                          >  2:4424937/1‑93 (MQ=255)

                  ggCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGc                                                                       <  1:1600798/93‑1 (MQ=255)

                  ggCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGc                                                                       >  2:636646/1‑93 (MQ=255)

                   gCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCg                                                                      >  2:2260275/1‑93 (MQ=255)

                   gCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCg                                                                      >  1:1036658/1‑93 (MQ=255)

                   gCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCg                                                                      >  2:1742468/1‑93 (MQ=255)

                     cACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGt                                                                    >  2:1904576/1‑93 (MQ=255)

                     cACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGt                                                                    <  1:1374756/93‑1 (MQ=255)

                     cACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGt                                                                    >  2:1157315/1‑93 (MQ=255)

                      aCTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGtt                                                                   >  2:1888841/1‑93 (MQ=255)

                      aCTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGtt                                                                   >  2:172146/1‑93 (MQ=255)

                      aCTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGtt                                                                   <  1:4230912/93‑1 (MQ=255)

                       cTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGttt                                                                  <  1:4034474/93‑1 (MQ=255)

                        tGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGtttt                                                                 >  2:4728018/1‑93 (MQ=255)

                        tGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGtttt                                                                 >  2:2402271/1‑93 (MQ=255)

                         ggAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTtc                                                                >  2:335023/1‑93 (MQ=255)

                           aaGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTtcct                                                              >  2:3876632/1‑93 (MQ=255)

                            aGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTtcctc                                                             >  2:2693400/1‑93 (MQ=255)

                             ggCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTtcctcc                                                            <  1:1224561/93‑1 (MQ=255)

                              gCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTcctcct                                                           <  2:2563020/93‑1 (MQ=255)

                              gCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTcctcct                                                           <  1:2368935/93‑1 (MQ=255)

                              gCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTcctcct                                                           <  1:3463249/93‑1 (MQ=255)

                               ccTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCctcctc                                                          <  1:2170575/93‑1 (MQ=255)

                                 tGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGa                                                        >  1:4222202/1‑93 (MQ=255)

                                 tGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGa                                                        >  2:1425032/1‑93 (MQ=255)

                                  gaagaaTTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGAt                                                       <  2:3598552/93‑1 (MQ=255)

                                    aggaTTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGAtgt                                                     <  1:4115152/90‑1 (MQ=255)

                                      aaTTGGGAAGGCGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGCCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGc                                                   <  2:35021/93‑1 (MQ=255)

                                      aaTTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGAtgtgt                                                   >  1:1643702/1‑92 (MQ=255)

                                       aTTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCt                                                  >  1:3968993/1‑93 (MQ=255)

                                        ttGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCtt                                                 >  1:309794/1‑93 (MQ=255)

                                        ttGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCtt                                                 <  1:2342106/93‑1 (MQ=255)

                                          gggAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCa                                               <  2:2898131/93‑1 (MQ=255)

                                           ggAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCAc                                              >  2:3617272/1‑93 (MQ=255)

                                           ggAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCAc                                              <  1:1432430/93‑1 (MQ=255)

                                              aGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGaa                                           <  2:2955518/93‑1 (MQ=255)

                                              aGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGaa                                           <  1:2825015/93‑1 (MQ=255)

                                              aGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGaa                                           <  1:2644055/93‑1 (MQ=255)

                                               gACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAAt                                          >  2:1408780/1‑93 (MQ=255)

                                                 cgcgGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAAttt                                        <  1:464341/93‑1 (MQ=255)

                                                     ggTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGc                                    >  2:716004/1‑93 (MQ=255)

                                                       tataAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCaa                                  >  1:4656028/1‑93 (MQ=255)

                                                       tataAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCaa                                  >  2:2426608/1‑93 (MQ=255)

                                                        ataAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAAt                                 <  1:558322/93‑1 (MQ=255)

                                                         taAAAACTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATa                                >  2:2272476/1‑93 (MQ=255)

                                                          aaaaaCTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATAc                               >  2:2490660/1‑93 (MQ=255)

                                                          aaaaaCTCTGGTAGCCATCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATAc                               <  2:2607995/93‑1 (MQ=255)

                                                           aaaaCTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACt                              >  1:1930191/1‑93 (MQ=255)

                                                            aaaCTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTg                             >  2:3714102/1‑93 (MQ=255)

                                                            aaaCTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTg                             >  1:3164097/1‑93 (MQ=255)

                                                              aCTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGAc                           <  2:2209815/93‑1 (MQ=255)

                                                              aCTCTGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGAc                           >  1:1358980/1‑93 (MQ=255)

                                                                tctGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCg                         <  1:3022678/93‑1 (MQ=255)

                                                                  tGGTAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGt                       >  1:2669250/1‑93 (MQ=255)

                                                                     tAGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGa                    <  1:4530521/93‑1 (MQ=255)

                                                                      aGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGAt                   >  1:1927925/1‑93 (MQ=255)

                                                                      aGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGAt                   >  2:1480044/1‑93 (MQ=255)

                                                                      aGCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGCTGAt                   >  1:2731134/1‑93 (MQ=255)

                                                                       gCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATa                  >  2:756909/1‑93 (MQ=255)

                                                                       gCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATa                  >  2:140061/1‑93 (MQ=255)

                                                                       gCCGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATa                  >  1:655281/1‑93 (MQ=255)

                                                                         cGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAAc                <  1:1268738/93‑1 (MQ=255)

                                                                         cGTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAAc                <  1:502686/93‑1 (MQ=255)

                                                                          gTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAAcc               >  2:1256427/1‑93 (MQ=255)

                                                                          gTCATCAAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAAcc               >  1:2032588/1‑93 (MQ=255)

                                                                           ttatcaAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCg              >  2:2042639/3‑93 (MQ=255)

                                                                           tcatcaAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCg              >  2:3870319/1‑93 (MQ=255)

                                                                            catcaAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGt             >  2:4395083/1‑93 (MQ=255)

                                                                             atcaAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGtt            >  1:3831769/1‑93 (MQ=255)

                                                                             atcaAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGtt            >  2:4694292/1‑93 (MQ=255)

                                                                              tcaAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTtc           >  1:4355426/1‑93 (MQ=255)

                                                                              tcaAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTtc           >  1:2595315/1‑93 (MQ=255)

                                                                               caAAAGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTtct          <  1:4660462/93‑1 (MQ=255)

                                                                                aaaaGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTtctc         <  1:2598430/93‑1 (MQ=255)

                                                                                aaaaGTCAGCGCTACAGCTTTTTCCGGTAGCGGGTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTtctc         >  1:1673738/1‑93 (MQ=255)

                                                                                 aaaGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTcg        >  1:4493709/1‑93 (MQ=255)

                                                                                  aaGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTcgc       >  2:1179023/1‑93 (MQ=255)

                                                                                  aaGTCAGCGCTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTcgc       >  1:3226840/1‑93 (MQ=255)

                                                                                       aGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGc  <  1:3128540/93‑1 (MQ=255)

ACGGGGGCCCATACAGCAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGC  >  NC_000913/1088662‑1088841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 23 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`