BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,601,614	(A)_7→6	coding (569/726 nt)	yjjQ →	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,601,608	0	A	.	99.6%	1032.7 / ‑2.8	225	coding (563/726 nt)	yjjQ	putative regulator
Reads supporting (aligned to +/- strand): ref base A (0/1); new base . (121/103); total (121/104)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 245 total aligned reads displayed.

GAGCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTT  >  NC_000913/4601517‑4601696

  gcgcAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaa                                                                                       <  1:2065218/93‑1 (MQ=255)

    gcAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaa                                                                                     <  1:2730110/93‑1 (MQ=255)

      aaGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaaca                                                                                   >  2:16072/1‑91 (MQ=255)

       aGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaacag                                                                                  >  1:923837/1‑90 (MQ=255)

        ggTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaacagt                                                                                 >  1:738076/1‑89 (MQ=255)

            cTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAGAAACAGTGCGa                                                                            >  1:3108664/1‑93 (MQ=255)

            cTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGa                                                                            >  1:4050213/1‑93 (MQ=255)

            cTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGa                                                                            <  1:1522767/93‑1 (MQ=255)

             tGCGTCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGGAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                           <  1:527058/93‑1 (MQ=255)

             tGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                           <  1:939747/93‑1 (MQ=255)

             tGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                           >  1:2524339/1‑93 (MQ=255)

             tGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                           <  1:2820735/93‑1 (MQ=255)

             tGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                           >  2:1830334/1‑93 (MQ=255)

             tGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                           <  1:1348118/93‑1 (MQ=255)

             tGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                           >  2:3604870/1‑93 (MQ=255)

              gCGGCTATTAGTTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCTC‑AAAAAACAGTGCGATg                                                                          >  2:301714/1‑93 (MQ=255)

               cGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGc                                                                         >  2:2969172/1‑93 (MQ=255)

                ggCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCg                                                                        >  2:19029/1‑93 (MQ=255)

                 gCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGt                                                                       >  2:1502681/1‑93 (MQ=255)

                 gCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGt                                                                       >  1:2365970/1‑93 (MQ=255)

                   tATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCg                                                                     >  1:1181987/1‑93 (MQ=255)

                    aTTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGa                                                                    >  2:4393608/1‑93 (MQ=255)

                      tAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACt                                                                  <  2:1781872/93‑1 (MQ=255)

                        ggTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTgg                                                                >  1:2543554/1‑93 (MQ=255)

                           aaaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCaa                                                             >  1:2674178/1‑93 (MQ=255)

                           aaaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCaa                                                             <  2:570375/93‑1 (MQ=255)

                            aaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAAt                                                            <  1:1264109/93‑1 (MQ=255)

                             aGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAAtt                                                           >  1:3843465/1‑93 (MQ=255)

                              ggATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTc                                                          >  2:1050923/1‑93 (MQ=255)

                              ggATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTc                                                          <  2:3254255/93‑1 (MQ=255)

                                 tGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTcaca                                                       <  2:1165181/93‑1 (MQ=255)

                                  ggggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAg                                                      >  1:2394758/1‑93 (MQ=255)

                                  ggggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAg                                                      >  2:3782614/1‑93 (MQ=255)

                                     gCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCaa                                                   <  1:3543821/93‑1 (MQ=255)

                                       aTCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACg                                                 >  1:2280450/1‑93 (MQ=255)

                                        tCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGc                                                >  2:258717/1‑93 (MQ=255)

                                        tCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGc                                                >  2:1398712/1‑93 (MQ=255)

                                         cAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCt                                               <  1:3692099/93‑1 (MQ=255)

                                         cAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCt                                               <  2:2156802/93‑1 (MQ=255)

                                          aaCTAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AATAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAATGCTg                                              >  2:4200336/1‑93 (MQ=255)

                                             cAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGaaa                                           <  1:1095529/93‑1 (MQ=255)

                                               gATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATg                                         >  2:584230/1‑93 (MQ=255)

                                               gATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATg                                         >  2:2432923/1‑93 (MQ=255)

                                                 tAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTa                                       <  2:4123525/93‑1 (MQ=255)

                                                 tAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTa                                       <  1:2920938/93‑1 (MQ=255)

                                                  aGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTAt                                      <  2:451309/93‑1 (MQ=255)

                                                    cTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGc                                    >  2:1482156/1‑93 (MQ=255)

                                                     ttCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCa                                   >  2:4175723/1‑93 (MQ=255)

                                                      tCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCAt                                  >  1:3165916/1‑93 (MQ=255)

                                                       cATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATg                                 >  1:1566607/1‑93 (MQ=255)

                                                        aTTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATgg                                >  2:3671236/1‑93 (MQ=255)

                                                        aTTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATgg                                <  1:462300/93‑1 (MQ=255)

                                                           gCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATa                             >  2:1245384/1‑93 (MQ=255)

                                                           gCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATa                             <  2:1337662/93‑1 (MQ=255)

                                                             ttAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaaa                           <  2:3367198/93‑1 (MQ=255)

                                                             ttAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaaa                           >  1:3357005/1‑93 (MQ=255)

                                                             ttAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaaa                           >  1:1406344/1‑93 (MQ=255)

                                                               aaGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATa                         <  2:3536058/93‑1 (MQ=255)

                                                               aaGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATa                         <  1:54132/93‑1 (MQ=255)

                                                               aaGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATa                         <  2:3891645/93‑1 (MQ=255)

                                                               aaGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATa                         <  1:1384946/93‑1 (MQ=255)

                                                                aGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAg                        >  2:796444/1‑93 (MQ=255)

                                                                 gAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgc                       >  1:824087/1‑93 (MQ=255)

                                                                   aaaGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgc                     <  1:2217128/93‑1 (MQ=255)

                                                                    aaGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgcg                    >  2:1988115/1‑93 (MQ=255)

                                                                    aaGTAATAAAACTATCAGCACCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgcg                    >  1:4645224/1‑93 (MQ=255)

                                                                      gTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgagcgcg                  >  2:2663023/1‑92 (MQ=255)

                                                                      gTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCa                  >  2:3136977/1‑93 (MQ=255)

                                                                       taataaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAg                 <  2:3943298/93‑1 (MQ=255)

                                                                       taataaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAg                 <  2:2470184/93‑1 (MQ=255)

                                                                        tataaAACCATCAGCTCCC‑AAAAAACAGTGCGATGCGTCGCCTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAgg                <  1:1882862/92‑1 (MQ=255)

                                                                        aataaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAgg                <  2:265798/93‑1 (MQ=255)

                                                                        aataaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAgg                <  1:4607307/93‑1 (MQ=255)

                                                                        aataaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAgg                <  1:193273/93‑1 (MQ=255)

                                                                        aataaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAgg                >  1:1762601/1‑93 (MQ=255)

                                                                          taaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAgggg              >  2:3754346/1‑92 (MQ=255)

                                                                          taaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGt              <  2:162616/93‑1 (MQ=255)

                                                                          taaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGt              <  1:4127673/93‑1 (MQ=255)

                                                                          taaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGt              <  2:4656546/93‑1 (MQ=255)

                                                                            aaaCTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGc            >  1:4304675/1‑93 (MQ=255)

                                                                             aaCTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCa           >  2:2743850/1‑93 (MQ=255)

                                                                             aaCTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCa           <  1:4412771/93‑1 (MQ=255)

                                                                              aCTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCaa          <  2:843115/93‑1 (MQ=255)

                                                                              aCTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCaa          <  2:2418787/93‑1 (MQ=255)

                                                                              aCTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCaa          >  1:2058596/1‑93 (MQ=255)

                                                                               cTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAag         <  1:1235822/93‑1 (MQ=255)

                                                                               cTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAag         >  2:4427836/1‑93 (MQ=255)

                                                                                tATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAaga        >  1:4452227/1‑93 (MQ=255)

                                                                                tATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAaga        <  1:1719341/93‑1 (MQ=255)

                                                                                 aTCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAagag       <  1:702000/93‑1 (MQ=255)

                                                                                 aTCAGCGCCC‑AAAAAACAGCGCGACGCGGCGACTGGCAATTCACAGCAACGCTGAAAGGAATGCATGGAGAAATAGTCCGCAGGGCGGAagag       >  2:399598/1‑93 (MQ=255)

                                                                                  tCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAgaga      <  1:2477641/93‑1 (MQ=255)

                                                                                  tCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAgaga      >  2:234845/1‑93 (MQ=255)

                                                                                   cAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaa     <  1:1858974/93‑1 (MQ=255)

                                                                                   cAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaa     >  2:1841962/1‑93 (MQ=255)

                                                                                   cAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaa     <  1:3782157/93‑1 (MQ=255)

                                                                                     gcgcCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaact   >  2:3455213/1‑93 (MQ=255)

                                                                                     gcgcCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaact   >  1:575372/1‑93 (MQ=255)

                                                                                     gcgcCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaact   >  1:4004707/1‑93 (MQ=255)

                                                                                      cgcCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaactt  <  2:512982/93‑1 (MQ=255)

GAGCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTT  >  NC_000913/4601517‑4601696

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