BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,340,190	A→G	T93A (ACT→GCT)	yrbI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,340,190	0	A	G	99.5%	702.3 / NA	218	T93A (ACT→GCT)	yrbI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (110/107); minor base T (0/1); total (110/108)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

GTTGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCACTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTAT  >  NC_000913/3340101‑3340279

gTTGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctc                                                                                        >  2:5833695/1‑93 (MQ=255)

 ttttgCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctca                                                                                       <  1:54016/90‑1 (MQ=255)

 ttGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctca                                                                                       <  2:4221774/93‑1 (MQ=255)

  tgtgCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctcac                                                                                      >  1:4247543/1‑93 (MQ=255)

     gcgcTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTg                                                                                   <  1:5826144/93‑1 (MQ=255)

         tCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATc                                                                               >  1:2514470/1‑93 (MQ=255)

          cACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCa                                                                              <  1:1751870/93‑1 (MQ=255)

           cccTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAg                                                                             <  1:2442571/92‑1 (MQ=255)

           aCCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAg                                                                             <  2:3706671/93‑1 (MQ=255)

            ccTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgg                                                                            <  1:5539453/93‑1 (MQ=255)

             ctctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAggg                                                                           >  1:3348771/1‑93 (MQ=255)

             ctctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAggg                                                                           >  2:4044358/1‑93 (MQ=255)

              tctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgggg                                                                          >  2:4486033/1‑93 (MQ=255)

              tctGATATTGAAGTCGCTATCATTAACGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgggg                                                                          >  2:2846665/1‑93 (MQ=255)

                  atatTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGt                                                                      <  2:638301/93‑1 (MQ=255)

                  atatTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGt                                                                      >  1:3991121/1‑93 (MQ=255)

                   tatTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTc                                                                     >  2:4658751/1‑93 (MQ=255)

                   tatTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTc                                                                     <  1:1558889/93‑1 (MQ=255)

                      tGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTcaaa                                                                  <  1:6219252/93‑1 (MQ=255)

                       gAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTcaaac                                                                 <  2:534372/93‑1 (MQ=255)

                        aaGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTcaaaca                                                                <  2:2820037/93‑1 (MQ=255)

                          gTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTcaaacaaa                                                              >  2:4195081/1‑93 (MQ=255)

                           tCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCaaacaaac                                                             >  1:1863150/1‑93 (MQ=255)

                            cGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACt                                                            >  1:5978969/1‑93 (MQ=255)

                             gCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAAAAAACTg                                                           >  1:4312358/1‑93 (MQ=255)

                              cTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGa                                                          <  2:2193580/93‑1 (MQ=255)

                              cTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGa                                                          <  1:1600529/93‑1 (MQ=255)

                              cTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGa                                                          >  2:1156257/1‑93 (MQ=255)

                              cTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGa                                                          <  1:648905/93‑1 (MQ=255)

                               tATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGAt                                                         >  1:1328127/1‑93 (MQ=255)

                               tATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGAt                                                         <  2:501344/93‑1 (MQ=255)

                                aTCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAAATGATc                                                        >  2:5202843/1‑93 (MQ=255)

                                aTCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGCGGCAGACAAACAAACTGATc                                                        >  2:5748468/1‑93 (MQ=255)

                                 tCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCg                                                       >  1:1042558/1‑93 (MQ=255)

                                  cATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGc                                                      >  2:6289870/1‑93 (MQ=255)

                                  cATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGc                                                      <  2:386608/93‑1 (MQ=255)

                                    ttACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCt                                                    <  1:2781096/93‑1 (MQ=255)

                                    ttACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCt                                                    >  2:2890121/1‑93 (MQ=255)

                                    ttACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCt                                                    >  2:2991989/1‑93 (MQ=255)

                                     tACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCtt                                                   <  1:1018665/93‑1 (MQ=255)

                                     tACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCtt                                                   >  1:3229243/1‑93 (MQ=255)

                                       ccGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTa                                                 <  2:1551053/93‑1 (MQ=255)

                                        cGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAg                                                >  1:5141863/1‑93 (MQ=255)

                                        cGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAg                                                >  2:4927561/1‑93 (MQ=255)

                                          ggCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCg                                              <  1:3830251/93‑1 (MQ=255)

                                           gCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGa                                             <  1:2684638/93‑1 (MQ=255)

                                           gCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGa                                             <  2:2307388/93‑1 (MQ=255)

                                           gCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGa                                             <  2:1658002/93‑1 (MQ=255)

                                           gCGAAAGGCTAAACTGGTAGAAGATCGTGGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGa                                             <  2:2078163/93‑1 (MQ=255)

                                             gAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATc                                           >  1:920964/1‑93 (MQ=255)

                                              aaaGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATct                                          >  1:2601067/1‑93 (MQ=255)

                                              aaaGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATct                                          >  1:5214304/1‑93 (MQ=255)

                                              aaaGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATct                                          <  2:5538551/93‑1 (MQ=255)

                                              aaaGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATct                                          <  1:3084145/93‑1 (MQ=255)

                                                aGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATctgc                                        >  1:4652756/1‑93 (MQ=255)

                                                 ggCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATctgct                                       >  1:3493270/1‑93 (MQ=255)

                                                 ggCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATctgct                                       <  2:851361/93‑1 (MQ=255)

                                                  gCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATctgctg                                      >  1:4677082/1‑93 (MQ=255)

                                                  gCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCCATctgctg                                      >  1:1200738/1‑93 (MQ=255)

                                                    tAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGa                                    <  2:4905441/93‑1 (MQ=255)

                                                    tAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGa                                    >  2:4301631/1‑93 (MQ=255)

                                                      aaCTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaaa                                  >  2:4512446/1‑93 (MQ=255)

                                                      aaCTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaaa                                  >  1:2114166/1‑93 (MQ=255)

                                                      aaCTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaaa                                  >  2:2463456/1‑93 (MQ=255)

                                                       aCTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaaaa                                 >  2:2267988/1‑93 (MQ=255)

                                                        cTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaaaaa                                >  2:3394727/1‑93 (MQ=255)

                                                         tGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAAc                               <  1:8088/93‑1 (MQ=255)

                                                          ggTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACt                              <  2:2668851/93‑1 (MQ=255)

                                                           gTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTg                             <  1:3948385/93‑1 (MQ=255)

                                                           gTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTg                             <  2:5185279/93‑1 (MQ=255)

                                                            tAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTgg                            >  1:5698693/1‑93 (MQ=255)

                                                            tAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTgg                            >  1:1080234/1‑93 (MQ=255)

                                                            tAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCACCTTTAGCGATCTGCTGGAAAAACTgg                            <  1:6078393/93‑1 (MQ=255)

                                                              gaagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCg                          <  1:4549126/93‑1 (MQ=255)

                                                              gaagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCg                          >  2:1797039/1‑93 (MQ=255)

                                                               aagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGa                         <  1:2064573/93‑1 (MQ=255)

                                                               aagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGa                         <  2:480607/93‑1 (MQ=255)

                                                                agaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGAt                        <  1:4020437/93‑1 (MQ=255)

                                                                  aTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTg                      >  1:1388757/1‑93 (MQ=255)

                                                                  aTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTg                      >  1:2251824/1‑93 (MQ=255)

                                                                  aTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTg                      >  2:1994602/1‑93 (MQ=255)

                                                                   tCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGc                     >  2:2913065/1‑93 (MQ=255)

                                                                     gTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGccc                   >  1:1536648/1‑93 (MQ=255)

                                                                     gTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGccc                   >  1:4733587/1‑93 (MQ=255)

                                                                     gTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGccc                   <  1:86716/93‑1 (MQ=255)

                                                                     gTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGccc                   <  1:5749454/93‑1 (MQ=255)

                                                                         tgCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGa               >  2:5410517/1‑93 (MQ=255)

                                                                           ccACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGaaa             <  1:4093928/93‑1 (MQ=255)

                                                                             acaTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAt           <  1:532186/93‑1 (MQ=255)

                                                                             acaTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAt           <  2:3608621/93‑1 (MQ=255)

                                                                              caTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtg          >  2:3888188/1‑93 (MQ=255)

                                                                              caTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtg          <  1:3705356/93‑1 (MQ=255)

                                                                               aTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtgt         >  1:1679932/1‑93 (MQ=255)

                                                                               aTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtgt         <  1:4782624/93‑1 (MQ=255)

                                                                                   gggATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAACTGTGGCt     >  1:1424449/1‑93 (MQ=255)

                                                                                    ggATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGGTCTGCTGGAAAAACTGGCGATTGCCCCGGAGAATGTGGcct    >  2:5890961/1‑91 (MQ=255)

                                                                                      aTCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTAt  <  1:5222875/93‑1 (MQ=255)

                                                                                      aTCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTAt  <  1:4454051/93‑1 (MQ=255)

                                                                                      aTCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTAt  <  2:3675555/93‑1 (MQ=255)

                                                                                      aTCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTAt  <  2:2413173/93‑1 (MQ=255)

GTTGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCACTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTAT  >  NC_000913/3340101‑3340279

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