BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,991,845	T→C	N62N (AAT→AAC)	b2853 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,991,845	0	T	C	99.5%	714.8 / ‑4.9	216	N62N (AAT→AAC)	b2853	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (100/115); total (100/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 220 total aligned reads displayed.

CACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTT  >  NC_000913/2991755‑2991934

cacCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAAccc                                                                                         <  1:3924834/93‑1 (MQ=255)

  ccAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTa                                                                                       >  1:1022411/1‑93 (MQ=255)

  ccAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTa                                                                                       >  2:3732743/1‑93 (MQ=255)

   cAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAt                                                                                      <  2:5484507/93‑1 (MQ=255)

     gTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtct                                                                                    >  2:2666086/1‑93 (MQ=255)

     gTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtct                                                                                    <  1:4515017/93‑1 (MQ=255)

      tGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtctc                                                                                   <  1:4223591/93‑1 (MQ=255)

      tGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtctc                                                                                   <  1:4070214/93‑1 (MQ=255)

           cTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCtt                                                                              <  1:518248/93‑1 (MQ=255)

            tGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTa                                                                             >  1:627398/1‑93 (MQ=255)

             gATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTAc                                                                            <  1:1407419/93‑1 (MQ=255)

               tACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACta                                                                          >  2:5084558/1‑93 (MQ=255)

                aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACtat                                                                         >  1:6140351/1‑93 (MQ=255)

                aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACtat                                                                         <  1:3980170/93‑1 (MQ=255)

                  aGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATaa                                                                       >  1:6014594/1‑93 (MQ=255)

                  aGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATaa                                                                       >  2:6288906/1‑93 (MQ=255)

                   gATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATaaa                                                                      >  2:3918042/1‑93 (MQ=255)

                   gATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATaaa                                                                      >  2:6176296/1‑93 (MQ=255)

                     tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAAt                                                                    >  1:3276226/1‑93 (MQ=255)

                       aaaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTg                                                                  >  1:4392923/1‑93 (MQ=255)

                         aTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAg                                                                >  1:729291/1‑93 (MQ=255)

                          tACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGt                                                               >  2:3126104/1‑93 (MQ=255)

                          tACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGt                                                               >  2:1785341/1‑93 (MQ=255)

                             aaGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGt                                                            <  1:5227326/93‑1 (MQ=255)

                             aaGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGt                                                            <  2:1536558/93‑1 (MQ=255)

                             aaGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGt                                                            >  1:2199177/1‑93 (MQ=255)

                              aGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTACAAAATTGAGTGGTa                                                           <  1:4092237/93‑1 (MQ=255)

                              aGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTa                                                           <  2:6384914/93‑1 (MQ=255)

                              aGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTa                                                           >  1:554766/1‑93 (MQ=255)

                               gCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAg                                                          <  1:3237561/93‑1 (MQ=255)

                                 aTTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGAc                                                        <  2:2227445/93‑1 (MQ=255)

                                 aTTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGGAGAc                                                        >  2:1366185/1‑93 (MQ=255)

                                  tttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACg                                                       >  2:5866529/1‑93 (MQ=255)

                                   ttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGa                                                      >  1:5294562/1‑93 (MQ=255)

                                    tataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGAt                                                     <  2:4768041/93‑1 (MQ=255)

                                     ataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGAtt                                                    >  2:4461399/1‑93 (MQ=255)

                                         tttCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGAcc                                                <  1:1522747/93‑1 (MQ=255)

                                          ttCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCa                                               >  2:121639/1‑93 (MQ=255)

                                          ttCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCa                                               >  1:2694047/1‑93 (MQ=255)

                                           tCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaa                                              <  2:802418/93‑1 (MQ=255)

                                             aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaat                                            >  2:2716701/1‑93 (MQ=255)

                                             aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaat                                            >  1:4027012/1‑93 (MQ=255)

                                             aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaat                                            <  1:1704924/93‑1 (MQ=255)

                                              aCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaata                                           >  2:4753219/1‑93 (MQ=255)

                                              aCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaata                                           >  2:5292684/1‑93 (MQ=255)

                                                 tggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaataaat                                        >  1:500983/1‑93 (MQ=255)

                                                 tggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaataaat                                        <  1:4985305/93‑1 (MQ=255)

                                                  ggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATc                                       >  2:4244566/1‑93 (MQ=255)

                                                   gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCa                                      <  2:3045460/93‑1 (MQ=255)

                                                    atggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCAt                                     >  1:6597456/1‑93 (MQ=255)

                                                    atggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCAt                                     <  1:4547550/93‑1 (MQ=255)

                                                     tggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATa                                    <  2:5002831/93‑1 (MQ=255)

                                                     tggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATa                                    <  2:6098532/93‑1 (MQ=255)

                                                      ggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATaa                                   <  1:4726694/93‑1 (MQ=255)

                                                      ggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATaa                                   <  2:1299610/93‑1 (MQ=255)

                                                         gcCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACga                                <  2:6463157/93‑1 (MQ=255)

                                                         gcCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACga                                <  2:3635648/93‑1 (MQ=255)

                                                         gcCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACga                                >  2:3366547/1‑93 (MQ=255)

                                                           cTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACgaga                              <  2:5800130/93‑1 (MQ=255)

                                                            ttttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGAt                             >  2:4506148/1‑93 (MQ=255)

                                                            ttttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGAt                             >  1:5031889/1‑93 (MQ=255)

                                                            ttttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGAt                             >  1:5958953/1‑93 (MQ=255)

                                                            ttttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGAt                             >  1:943210/1‑93 (MQ=255)

                                                              ttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATAc                           <  2:745884/93‑1 (MQ=255)

                                                              ttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATAc                           <  1:4957435/93‑1 (MQ=255)

                                                               tCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACt                          <  1:604800/93‑1 (MQ=255)

                                                               tCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACt                          <  1:278610/93‑1 (MQ=255)

                                                                cAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtt                         <  1:1237838/93‑1 (MQ=255)

                                                                cAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtt                         <  1:5711311/93‑1 (MQ=255)

                                                                cAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtt                         <  1:1869577/93‑1 (MQ=255)

                                                                cAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtt                         >  1:6354086/1‑93 (MQ=255)

                                                                 aaaaTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACttt                        >  2:679994/1‑93 (MQ=255)

                                                                   aaTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTta                      >  2:5578472/1‑93 (MQ=255)

                                                                    aTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtat                     <  2:198628/93‑1 (MQ=255)

                                                                    aTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtat                     >  1:135292/1‑93 (MQ=255)

                                                                     tGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGAAGACGAGTGACCAAATAAATCATAACGAGATACTATtata                    >  1:6722822/1‑93 (MQ=255)

                                                                      gATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataa                   >  1:851159/1‑93 (MQ=255)

                                                                       atatGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                  <  1:6405546/93‑1 (MQ=255)

                                                                        tatGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataata                 <  2:1117131/93‑1 (MQ=255)

                                                                         atGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                <  1:3574541/93‑1 (MQ=255)

                                                                         atGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                <  2:6574291/93‑1 (MQ=255)

                                                                           gTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa              >  1:6390534/1‑93 (MQ=255)

                                                                           gTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa              <  1:5466590/93‑1 (MQ=255)

                                                                           gTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa              <  1:2133321/93‑1 (MQ=255)

                                                                           gTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa              >  1:1762517/1‑93 (MQ=255)

                                                                            tATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAt             <  1:2923809/93‑1 (MQ=255)

                                                                            tATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAt             <  1:4294728/93‑1 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc            <  2:3807048/93‑1 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc            <  2:5701709/93‑1 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc            <  1:5870171/93‑1 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc            <  2:1285882/93‑1 (MQ=255)

                                                                               tCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtctc          <  2:1855328/93‑1 (MQ=255)

                                                                                 agaagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt        >  1:3349077/1‑93 (MQ=255)

                                                                                  gaagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt       >  1:3474536/1‑93 (MQ=255)

                                                                                  gaagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt       >  1:2328854/1‑93 (MQ=255)

                                                                                  gaagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt       >  2:3283186/1‑93 (MQ=255)

                                                                                   aagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTc      >  1:2988952/1‑93 (MQ=255)

                                                                                     gaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCt    <  2:2369776/93‑1 (MQ=255)

                                                                                     gaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCt    >  2:6598668/1‑93 (MQ=255)

                                                                                      aTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCtt   <  2:4651716/93‑1 (MQ=255)

CACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTT  >  NC_000913/2991755‑2991934

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