BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,006,483	+C	coding (183/207 nt)	b4554 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,006,480	1	.	C	94.7%	290.7 / 12.4	114	coding (186/207 nt)	b4554	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (5/1); new base C (62/46); total (67/47)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 127 total aligned reads displayed.

TTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑C‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCG  >  NC_000913/4006405‑4006550

tgcccAGCCAAAGCAAGGAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑t                                                                            <  2:654328/74‑1 (MQ=255)

 tCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑c                                                                          >  2:2066692/1‑76 (MQ=255)

  cccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACTCTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑ca                                                                         <  2:2278909/76‑2 (MQ=255)

  cccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑c‑c                                                                        <  2:3405314/76‑1 (MQ=255)

  cccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑c‑c                                                                        <  2:4306641/76‑1 (MQ=255)

   ccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑c‑c‑c                                                                      >  1:428170/1‑76 (MQ=255)

    cAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑c‑c‑cc                                                                     <  2:2360814/76‑2 (MQ=255)

    cAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑c‑c‑cc                                                                     <  1:4832196/76‑2 (MQ=255)

     aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑c‑c‑cca                                                                    <  1:1337784/76‑3 (MQ=255)

     aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑c‑c‑cca                                                                    >  1:3804448/1‑74 (MQ=255)

         aaaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑C‑C‑C‑CAaccc                                                               >  2:1582440/1‑76 (MQ=255)

          aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AAcccc                                                               <  1:3106484/76‑1 (MQ=255)

            tCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCAc                                                             <  2:3505258/76‑1 (MQ=255)

            acAAGGAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCAc                                                             <  2:4619637/75‑1 (MQ=255)

             cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACg                                                            >  2:4456655/1‑76 (MQ=255)

             cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACg                                                            <  2:1762764/76‑1 (MQ=255)

                gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCa                                                         <  1:2735372/76‑1 (MQ=255)

                 tAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCaa                                                        >  2:3789454/1‑76 (MQ=255)

                   aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAAcc                                                      <  2:3558186/76‑1 (MQ=255)

                    gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCa                                                     <  1:3866135/76‑1 (MQ=255)

                    gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCa                                                     <  1:2441443/76‑1 (MQ=255)

                     ccGCGCCTGCCCACTTCAACACTTCAAGCGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCACCCAg                                                    <  1:1020442/76‑1 (MQ=255)

                     ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCTCTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAg                                                    <  2:581122/76‑1 (MQ=255)

                     ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAg                                                    >  1:3625012/1‑76 (MQ=255)

                      cgcgCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGc                                                   >  1:2307628/1‑76 (MQ=255)

                      cgcgCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGc                                                   <  1:1529277/76‑1 (MQ=255)

                      cgcgCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGc                                                   >  2:1441266/1‑76 (MQ=255)

                        cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGcac                                                 >  2:797995/1‑76 (MQ=255)

                         gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGTGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCACcaca                                                >  2:1166474/1‑76 (MQ=255)

                          ccTGCCCTCTTCAACGCTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑C‑C‑CGAACCCCACGCCAACCAGCACaa                                               <  2:2546720/76‑1 (MQ=255)

                          ccTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACaa                                               >  1:2354702/1‑76 (MQ=255)

                           cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAat                                              >  1:968366/1‑76 (MQ=255)

                           cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAat                                              >  1:4274133/1‑76 (MQ=255)

                           cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAat                                              >  2:1423482/1‑76 (MQ=255)

                           cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAat                                              >  2:2043213/1‑76 (MQ=255)

                           cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAat                                              <  2:752938/76‑1 (MQ=255)

                            tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAata                                             >  2:2626489/1‑76 (MQ=255)

                             gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAatat                                            <  1:3226945/76‑1 (MQ=255)

                              cccACTTCAACGCTTCAAAGGCAATCGCTGAGCAGGAAAATAGCG‑TCC‑C‑C‑AACCCCACACCAGCCAGCACAATATg                                           <  1:1409546/76‑1 (MQ=255)

                              cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATg                                           >  1:2021648/1‑76 (MQ=255)

                               ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGa                                          >  1:1885553/1‑76 (MQ=255)

                               ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGa                                          >  1:1260146/1‑76 (MQ=255)

                               ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑CAC‑C‑AACCACACGCAAACCAGCACAATATGa                                          >  1:1165796/1‑76 (MQ=255)

                               ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑CCC‑C‑C‑ACCCCCACCCCACCCCCCACAATATGa                                          >  2:326556/1‑76 (MQ=255)

                               ccACTTCAACACTTAAAACGCAATAACTGAGAGGGAAAATAGCG‑T‑C‑CAC‑AACCCCACGCCAACCAGCACAATATGa                                          >  2:3038397/1‑76 (MQ=255)

                                 acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAAt                                        >  2:760696/1‑76 (MQ=255)

                                 acttcaacacttcaaACGCAATCACTAAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAAt                                        >  2:3925426/1‑76 (MQ=255)

                                   ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGAC‑C‑C‑C‑GACCCCACGCCAACCAGCACAATATGAATcc                                      >  1:2203441/1‑75 (MQ=255)

                                    tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGc                                     <  1:1630860/76‑1 (MQ=255)

                                    tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGc                                     <  1:3178305/76‑1 (MQ=255)

                                    tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGc                                     >  2:3061974/1‑76 (MQ=255)

                                     caacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGcc                                    >  2:4125468/1‑76 (MQ=255)

                                      aacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCa                                   <  1:1730599/76‑1 (MQ=255)

                                       acacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAg                                  <  1:427981/76‑1 (MQ=255)

                                         acttcaaACGCAATCACTGAGCGGGAAAATCGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACACTATGAATCGACAGTc                                >  1:678740/1‑76 (MQ=255)

                                         acttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                >  1:651376/1‑76 (MQ=255)

                                          cttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACGCCACGCCAACCAGCACAATATGAATCGCCAGTcc                               <  2:3625012/76‑1 (MQ=255)

                                           ttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCg                              >  1:3085487/1‑76 (MQ=255)

                                            tcaaACGCAATCACTGAGCGGGAAAATGTCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                             <  2:483094/76‑1 (MQ=255)

                                             caaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACACTATGAATCGCCAGACCggc                            >  2:365909/1‑75 (MQ=255)

                                             caaACGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                            <  1:4533948/76‑1 (MQ=255)

                                              aaaCGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                           <  1:4062251/76‑1 (MQ=255)

                                              aaaCGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                           >  2:3123706/1‑76 (MQ=255)

                                              aaaCGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                           <  1:3789454/76‑1 (MQ=255)

                                               agcGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                          <  2:4808746/74‑1 (MQ=255)

                                                aCGCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTg                         >  2:596832/1‑76 (MQ=255)

                                                  gCAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                       >  1:4383906/1‑76 (MQ=255)

                                                   cAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                      >  2:1009245/1‑76 (MQ=255)

                                                   cAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                      >  1:3825347/1‑76 (MQ=255)

                                                   cAATCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                      >  1:2021340/1‑76 (MQ=255)

                                                    aaTCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                     >  1:2241931/1‑76 (MQ=255)

                                                    aaTCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                     >  1:1566137/1‑76 (MQ=255)

                                                      tCACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGccc                   <  1:2770476/76‑1 (MQ=255)

                                                       cACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                  <  2:4141706/76‑1 (MQ=255)

                                                       cACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                  >  1:2902428/1‑76 (MQ=255)

                                                       cACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                  <  2:1902071/76‑1 (MQ=255)

                                                       cACTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                  >  1:2547571/1‑76 (MQ=255)

                                                       cACTGAGCGGGAAAATAGCG‑T‑CGC‑C‑AACCCCACAACAACCAGCATAATATTAACCATCAGTCCTGTCTGCAGcaca                  >  2:3245124/1‑76 (MQ=255)

                                                        aCTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAg                 >  1:991046/1‑76 (MQ=255)

                                                         cTGAGCGGGAAAATAGCG‑TCC‑C‑C‑ACCCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGc                <  1:1845408/76‑1 (MQ=255)

                                                         cTGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAga                >  1:1357451/1‑75 (MQ=255)

                                                          tGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa               <  1:4409927/76‑1 (MQ=255)

                                                          tGAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa               <  1:4456655/76‑1 (MQ=255)

                                                           gAGCGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa              <  1:2543629/76‑1 (MQ=255)

                                                              cGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAg           <  1:1360531/76‑1 (MQ=255)

                                                              cGGGAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAg           <  1:2606032/76‑1 (MQ=255)

                                                               gggAAAATGGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa          <  2:2251887/76‑1 (MQ=255)

                                                               gggAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa          >  1:3641504/1‑76 (MQ=255)

                                                               gggAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa          >  1:2809779/1‑76 (MQ=255)

                                                               gggAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa          >  2:2310662/1‑76 (MQ=255)

                                                                ggAAAATAGCG‑TCC‑C‑C‑ACCCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc         <  1:1229373/76‑1 (MQ=255)

                                                                 gAAAATAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCCGAAGCCCAGCAATAGACg        <  2:2448148/76‑1 (MQ=255)

                                                                   aaaTAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc      >  2:4678164/1‑76 (MQ=255)

                                                                    aaTAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCa     >  2:4453998/1‑76 (MQ=255)

                                                                    aaTAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCACTCTGAAACCCAGCAATACACTCCa     >  1:1444214/1‑76 (MQ=255)

                                                                    aaTAGAG‑TCC‑C‑C‑AACCCCACCCCAACCAGCACAATATGAATCGCCAGGCCGCTCGCAAGCCCAGCAAGAGACGCCa     >  2:2521817/1‑76 (MQ=255)

                                                                      tAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGCCGCcccc   >  2:3355548/1‑76 (MQ=38)

                                                                      tAGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc   >  1:3617996/1‑76 (MQ=255)

                                                                      caGCG‑TCC‑C‑C‑ACCCCCACGCCAACCAGCACAATATGACTCGCCAGCCCGGTCTGAAGCCCAGCAATAGACGCCAcc   >  2:558194/2‑76 (MQ=38)

                                                                       aGCG‑TCC‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg  <  1:2361183/76‑1 (MQ=255)

TTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑C‑C‑C‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCG  >  NC_000913/4006405‑4006550

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 10 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`