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	.	98.9%	1278.1 / 3.4	278	coding (563/726 nt)	yjjQ	putative regulator
Reads supporting (aligned to +/- strand): ref base A (2/1); new base . (161/114); total (163/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

Only 100 of 312 total aligned reads displayed.

AGAGCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTT  >  NC_000913/4601516‑4601696

agagCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCa                                                                                          >  1:1776912/1‑93 (MQ=255)

   gcgcAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaa                                                                                       >  1:6633440/1‑93 (MQ=255)

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

    cgcAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaa                                                                                      <  2:600707/93‑1 (MQ=255)

     gcAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaa                                                                                     >  2:4143446/1‑93 (MQ=255)

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

     gcAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaa                                                                                     >  2:4768220/1‑93 (MQ=255)

      cAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaac                                                                                    >  2:4912607/1‑92 (MQ=255)

      cAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaac                                                                                    <  2:4529797/93‑2 (MQ=255)

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

       aaGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaaca                                                                                   <  1:1272891/93‑3 (MQ=255)

        aGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaacag                                                                                  <  1:2422992/93‑4 (MQ=255)

         ggTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaacagt                                                                                 <  2:6507667/93‑5 (MQ=255)

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

          gTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTg                                                                               >  1:3368750/1‑93 (MQ=255)

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

                cGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGc                                                                         <  1:2940161/93‑1 (MQ=255)

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

                      ttAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGAc                                                                   >  1:3236265/1‑93 (MQ=255)

                        aGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTg                                                                 >  2:3351252/1‑93 (MQ=255)

                        aGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTg                                                                 <  1:1983658/93‑1 (MQ=255)

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

                           tAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCa                                                              <  1:338494/93‑1 (MQ=255)

                             aaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAAt                                                            >  2:1610958/1‑93 (MQ=255)

                             aaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAAt                                                            >  2:5905478/1‑93 (MQ=255)

                             aaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAAt                                                            >  1:101596/1‑93 (MQ=255)

                              aGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAAtt                                                           <  2:2291259/93‑1 (MQ=255)

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

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

                                gATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTca                                                         >  1:3943435/1‑93 (MQ=255)

                                gATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTca                                                         >  2:5201820/1‑93 (MQ=255)

                                 aTGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTcac                                                        <  2:3821704/93‑1 (MQ=255)

                                  tGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑ACAAAACAGTGCGATGCGTCGACTGGCAATTcaca                                                       >  1:6065952/1‑93 (MQ=255)

                                  tGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTcaca                                                       >  2:4082798/1‑93 (MQ=255)

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

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

                                       cATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAAc                                                  <  2:2667543/93‑1 (MQ=255)

                                         tCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGc                                                <  1:2316656/93‑1 (MQ=255)

                                         tCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGc                                                <  2:4668210/93‑1 (MQ=255)

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

                                         tCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGc                                                <  2:2364509/93‑1 (MQ=255)

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

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

                                           aaCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTg                                              >  2:1782898/1‑93 (MQ=255)

                                            aCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGa                                             >  1:5598964/1‑93 (MQ=255)

                                            aCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGa                                             <  2:1080693/93‑1 (MQ=255)

                                            aCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGa                                             >  1:4692354/1‑93 (MQ=255)

                                            aCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGa                                             >  2:5837762/1‑93 (MQ=255)

                                            aCCAGATAGCTTCATTGCTTAAGAAAAGCAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGa                                             >  2:1448924/1‑93 (MQ=255)

                                             ccAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGaa                                            >  2:199051/1‑93 (MQ=255)

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

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

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

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

                                                 aTAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGt                                        >  1:3706462/1‑93 (MQ=255)

                                                  tAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTa                                       >  2:1237849/1‑93 (MQ=255)

                                                   aGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTAt                                      >  2:5974334/1‑93 (MQ=255)

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

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

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

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

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

                                                        cATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATg                                 <  2:5426732/93‑1 (MQ=255)

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

                                                          ttGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGa                               <  1:4630103/93‑1 (MQ=255)

                                                             cTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaa                            >  2:802033/1‑93 (MQ=255)

                                                             cTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaa                            >  2:2258115/1‑93 (MQ=255)

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

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

                                                               tAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAAt                          <  1:6463780/93‑1 (MQ=255)

                                                                aaGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATa                         >  1:1725929/1‑93 (MQ=255)

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

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

                                                                    aaaGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgc                     >  1:1094208/1‑93 (MQ=255)

                                                                    aaaGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgc                     >  1:4713892/1‑93 (MQ=255)

                                                                    aaaGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgc                     >  2:2831827/1‑93 (MQ=255)

                                                                    aaaGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgc                     >  1:6358010/1‑93 (MQ=255)

                                                                    aaaGTAATAAAACTATCAGCGCCC‑AAAAAACAGCGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgc                     >  1:5647711/1‑93 (MQ=255)

                                                                      aGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgcgc                   >  2:2438797/1‑93 (MQ=255)

                                                                       gTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgcgcg                  >  2:5460243/1‑92 (MQ=255)

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

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

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

                                                                          ataaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAggg               <  2:4842382/93‑1 (MQ=255)

                                                                          ataaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAggg               >  2:1518336/1‑93 (MQ=255)

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

                                                                           taaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAgggc              <  1:5067421/93‑2 (MQ=255)

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

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

                                                                            aaaaCTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTg             >  1:4097175/1‑93 (MQ=255)

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

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

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

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

                                                                                  aTCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAagag       >  1:856194/1‑93 (MQ=255)

                                                                                  aTCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAagag       >  1:4271364/1‑93 (MQ=255)

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

                                                                                     aGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaac    >  2:6698858/1‑93 (MQ=255)

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

                                                                                      agcgccAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaac    <  2:3459718/89‑1 (MQ=255)

AGAGCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTT  >  NC_000913/4601516‑4601696

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