BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,389,803	A→G	P29P (CCT→CCC)	b1327 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,389,803	0	A	G	99.2%	770.5 / ‑3.8	248	P29P (CCT→CCC)	b1327	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (120/125); total (121/127)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 257 total aligned reads displayed.

GGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGA  >  NC_000913/1389713‑1389892

ggCATAGGTTCCGGCCGTTTGTTCTTGAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGgggg                                                                                         <  1:7252068/93‑1 (MQ=255)

  cATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGa                                                                                       >  1:5261918/1‑93 (MQ=255)

   aTAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGaa                                                                                      >  2:1163890/1‑93 (MQ=255)

   aTAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAAAACAATTGCCTGATACTGCCGGGTTTCATCCAATTTTGGGGGGaa                                                                                      >  1:1069406/1‑93 (MQ=255)

    tAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGaaa                                                                                     <  2:3459866/93‑1 (MQ=255)

    tAGGTTCCGGCCGTTTGTTCTTTAACGCACCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGaaa                                                                                     <  1:4129347/93‑1 (MQ=255)

     aGGTTCCGGTCGTTTGTTGTGTAACGTCCGCACCGGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAAt                                                                                    <  1:6351252/93‑1 (MQ=255)

     aGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAAt                                                                                    >  1:1004410/1‑93 (MQ=255)

     aGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAAt                                                                                    <  2:1269601/93‑1 (MQ=255)

     aGGTTCCGGCCGCTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAAt                                                                                    <  1:7681378/93‑1 (MQ=255)

        ttCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGa                                                                                 >  1:5985392/1‑93 (MQ=255)

         tCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATCGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGAt                                                                                <  2:7660999/93‑1 (MQ=255)

            ggCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCAc                                                                             <  2:7102645/93‑1 (MQ=255)

             gCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGTGGGAAATATATCACt                                                                            >  2:3084101/1‑93 (MQ=255)

             gCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACt                                                                            >  2:7466075/1‑93 (MQ=255)

             gCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACt                                                                            >  1:2504014/1‑93 (MQ=255)

             gCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACt                                                                            >  1:5585352/1‑93 (MQ=255)

               cGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGc                                                                          >  1:6682023/1‑93 (MQ=255)

               cGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGc                                                                          <  1:1232911/93‑1 (MQ=255)

               cGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGc                                                                          >  2:2419462/1‑93 (MQ=255)

                gTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCa                                                                         >  1:5478937/1‑93 (MQ=255)

                 tttGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAg                                                                        >  1:3829856/1‑93 (MQ=255)

                   tGTTCTTTATCGCCCCCTCCTGGGGGAGATAACACAATAGCCCGATGCAGACGGGTTTCATCAAATTTAGGAGGGAAATAGATCACTGCTGAc                                                                      <  2:6312148/93‑1 (MQ=255)

                     ttCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACaa                                                                    >  1:6785124/1‑93 (MQ=255)

                     ttCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACaa                                                                    <  1:7432773/93‑1 (MQ=255)

                        tttAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAaga                                                                 >  2:6554961/1‑93 (MQ=255)

                        tttAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAaga                                                                 <  2:2200580/93‑1 (MQ=255)

                        tttAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAaga                                                                 >  2:6364964/1‑93 (MQ=255)

                          tAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAgaga                                                               <  2:752918/93‑1 (MQ=255)

                          tAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAgaga                                                               <  2:1530011/93‑1 (MQ=255)

                          tAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAgaga                                                               >  1:5101342/1‑93 (MQ=255)

                          tAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAgaga                                                               <  2:1706897/93‑1 (MQ=255)

                          tAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAgaga                                                               >  1:4181869/1‑93 (MQ=255)

                           aaCGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGAt                                                              >  2:4600329/1‑93 (MQ=255)

                            aCGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATg                                                             <  1:6940210/93‑1 (MQ=255)

                              gCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGt                                                           >  2:2359173/1‑93 (MQ=255)

                              ccccccACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGt                                                           <  2:5542973/92‑1 (MQ=255)

                               cccccACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGtt                                                          <  2:7641195/93‑1 (MQ=255)

                                ccccACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTg                                                         >  2:6088707/1‑93 (MQ=255)

                                 cccACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTgg                                                        >  1:2113384/1‑93 (MQ=255)

                                  ccACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGa                                                       >  2:1664168/1‑93 (MQ=255)

                                   cACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGat                                                      >  1:4727612/1‑93 (MQ=255)

                                    aCCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGatt                                                     >  2:7215140/1‑93 (MQ=255)

                                     ccTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGatta                                                    <  2:7902352/93‑1 (MQ=255)

                                      cTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGattat                                                   >  1:4625137/1‑93 (MQ=255)

                                       tGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGattatt                                                  <  1:6741617/93‑1 (MQ=255)

                                       tGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGattatt                                                  >  2:2706366/1‑93 (MQ=255)

                                           tGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTg                                              <  2:7357241/93‑1 (MQ=255)

                                             agaTAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGtt                                            >  2:6518758/1‑93 (MQ=255)

                                             agaTAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGtt                                            >  2:5369197/1‑93 (MQ=255)

                                              gaTAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTa                                           <  2:4516340/93‑1 (MQ=255)

                                              aaTAACACAATGGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTa                                           <  2:4367457/92‑1 (MQ=255)

                                               aTAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAg                                          <  1:5650498/93‑1 (MQ=255)

                                               aTAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAg                                          <  2:2751021/93‑1 (MQ=255)

                                                tAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGt                                         >  1:5170066/1‑93 (MQ=255)

                                                tAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGt                                         <  1:3060314/93‑1 (MQ=255)

                                                 aaCACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTACAGACAAAGAGATGGTTGGATTATTGCTGTTAGTg                                        >  2:6872892/1‑93 (MQ=255)

                                                  agccaATTGCGTGATACTGACGGGTTTCATCAAATTTGGGGGGTAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGGTAGTGa                                       <  1:6424004/90‑1 (MQ=255)

                                                   cacaATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGaa                                      >  1:7591005/1‑93 (MQ=255)

                                                    acaATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAg                                     <  2:6999651/93‑1 (MQ=255)

                                                     caATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGc                                    >  1:6150343/1‑93 (MQ=255)

                                                      aaTTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCt                                   <  1:7736961/93‑1 (MQ=255)

                                                      aaTTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCt                                   <  1:4002454/93‑1 (MQ=255)

                                                       aTTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTg                                  >  1:821844/1‑93 (MQ=255)

                                                        ttGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGa                                 <  2:5261917/93‑1 (MQ=255)

                                                        ttGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGa                                 <  2:4394310/93‑1 (MQ=255)

                                                         tGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGAc                                >  1:1432879/1‑93 (MQ=255)

                                                         tGCATGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGAc                                <  1:1647922/93‑1 (MQ=255)

                                                           ccTGATACTGACGTGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACtt                              <  2:3263448/93‑1 (MQ=255)

                                                             tGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTtta                            >  2:3136139/1‑93 (MQ=255)

                                                             tGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTtta                            <  1:3709782/93‑1 (MQ=255)

                                                             tGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTtta                            >  1:5727877/1‑93 (MQ=255)

                                                              gATACTGACGGGTTTCATCAAATTTTGGGGGGAAATGGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttat                           >  1:2371728/1‑93 (MQ=255)

                                                                tACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttatta                         >  1:2835984/1‑93 (MQ=255)

                                                                tACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttatta                         >  1:5191679/1‑93 (MQ=255)

                                                                tACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttatta                         <  1:3322989/93‑1 (MQ=255)

                                                                  cTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTattatt                       <  1:497740/93‑1 (MQ=255)

                                                                  cTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTattatt                       >  1:996974/1‑93 (MQ=255)

                                                                   tGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtc                      <  1:4078407/93‑1 (MQ=255)

                                                                    gACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtca                     <  1:1897911/93‑1 (MQ=255)

                                                                     aCGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTAGTGCTGTTAGTGAAGCTGACTTTATTATtcat                    <  2:3412031/93‑1 (MQ=255)

                                                                      cGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatc                   <  1:4275737/93‑1 (MQ=255)

                                                                      cGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatc                   >  2:1853359/1‑93 (MQ=255)

                                                                       gggTTTCATCGAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatca                  <  1:1344160/93‑1 (MQ=255)

                                                                       gggTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatca                  >  1:2957963/1‑93 (MQ=255)

                                                                         gTTTGATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCAtt                <  2:4089125/93‑1 (MQ=255)

                                                                         gTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCAtt                <  1:7066108/93‑1 (MQ=255)

                                                                          tttCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTc               <  1:656769/93‑1 (MQ=255)

                                                                          tttCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTc               >  2:5876317/1‑93 (MQ=255)

                                                                          tttCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTc               <  1:6044437/93‑1 (MQ=255)

                                                                           ttCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTcc              <  2:3923784/93‑1 (MQ=255)

                                                                             catcaAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGt            >  2:499209/1‑93 (MQ=255)

                                                                              atcaAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGtt           <  2:5101341/93‑1 (MQ=255)

                                                                               tcaAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTc          <  2:639786/93‑1 (MQ=255)

                                                                                caAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTcc         >  1:3224842/1‑93 (MQ=255)

                                                                                  aaTTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTc       <  1:2704332/93‑1 (MQ=255)

                                                                                    ttttGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCAt     <  1:1931034/93‑1 (MQ=255)

                                                                                      ttGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATgg   >  2:3663895/1‑93 (MQ=255)

                                                                                      ttGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATgg   <  1:5805510/93‑1 (MQ=255)

                                                                                       tGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGa  >  1:3164194/1‑93 (MQ=255)

GGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGA  >  NC_000913/1389713‑1389892

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