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	97.0%	923.9 / 17.5	303	coding (186/207 nt)	b4554	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (6/3); new base C (154/140); total (160/143)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 329 total aligned reads displayed.

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑CCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

  cgcgcCACTGCTGGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑cc                                                                                       <  2:7070429/93‑1 (MQ=255)

    cgcCACTGCTGGATTTCCAGCCAACTCCAGTAAGCCGCGCCTGCCCACTTCAACACTTCACCCGCACTCCCTGAGCGGGAAAATTGCG‑T‑cccc                                                                                     <  2:6032786/93‑2 (MQ=255)

      ccACTGCTGGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑ccccaa                                                                                   >  2:2218937/1‑90 (MQ=255)

       cACTGCTGGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑ccccaac                                                                                  >  1:5920877/1‑89 (MQ=255)

       cACTGCTGGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑ccccaac                                                                                  <  2:18422/93‑5 (MQ=255)

       cACTGCTGGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑ccccaac                                                                                  >  1:2288199/1‑89 (MQ=255)

         ctgctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑CCCCAaccc                                                                                >  1:4257073/1‑93 (MQ=255)

         ctgctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑CCCCAaccc                                                                                <  2:4069213/93‑1 (MQ=255)

          tgctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAAcccc                                                                                >  2:3005261/1‑93 (MQ=255)

           gctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCa                                                                               >  2:6384627/1‑93 (MQ=255)

           gctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCa                                                                               >  1:6514880/1‑93 (MQ=255)

            ctgGATTTCCAGCCAAATCAAGTAAGGCACGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCag                                                                              >  2:2952040/1‑92 (MQ=255)

            ctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCAc                                                                              <  1:4899381/93‑1 (MQ=255)

            ctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCAc                                                                              <  1:1711167/93‑1 (MQ=255)

             tgGATTTCCAGCCAAATCAAGTAAGCCGCGCCGGCCCACTTCAACACGTCAAACGCAATCACTGAGCGGGAAAATAGCGAC‑CCCAACCCCACg                                                                             >  1:7033378/1‑93 (MQ=255)

              ggATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGc                                                                            >  1:6598816/1‑93 (MQ=255)

               gATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGcc                                                                           <  2:2335422/93‑1 (MQ=255)

               gATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGcc                                                                           >  2:192530/1‑93 (MQ=255)

               gATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGcc                                                                           >  2:2087130/1‑93 (MQ=255)

                aTTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCa                                                                          >  2:3558465/1‑93 (MQ=255)

                   tccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAAcc                                                                       <  1:3093050/92‑1 (MQ=255)

                     cAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACGGAGCGGGAAAATAGCG‑TCCCCCACCCCACGCCAACCac                                                                     >  1:2682203/1‑92 (MQ=255)

                       gCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGca                                                                   >  1:5598402/1‑93 (MQ=255)

                       gCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGca                                                                   >  1:2823643/1‑93 (MQ=255)

                       gCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGca                                                                   >  1:2157460/1‑93 (MQ=255)

                       gCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGca                                                                   <  2:3998925/93‑1 (MQ=255)

                       gCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGca                                                                   <  1:5779567/93‑1 (MQ=255)

                          aaaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACaa                                                                >  1:6370531/1‑93 (MQ=255)

                           aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAat                                                               <  1:4591115/93‑1 (MQ=255)

                            aTCAAGTAGGCCGCGCCTCCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAata                                                              <  1:7106898/93‑1 (MQ=255)

                            aTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAata                                                              <  1:6294128/93‑1 (MQ=255)

                            aTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAata                                                              <  1:4828848/93‑1 (MQ=255)

                            aTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAata                                                              >  2:1750813/1‑93 (MQ=255)

                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGaa                                                          <  1:2370943/93‑1 (MQ=255)

                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGaa                                                          <  2:3656347/93‑1 (MQ=255)

                                 gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAAt                                                         >  2:5067196/1‑93 (MQ=255)

                                 gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAAt                                                         >  1:6458395/1‑93 (MQ=255)

                                     gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGcc                                                     >  1:4999728/1‑93 (MQ=255)

                                     gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGcc                                                     <  2:4631340/93‑1 (MQ=255)

                                      ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCa                                                    >  2:3416019/1‑93 (MQ=255)

                                      ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCa                                                    <  1:1532463/93‑1 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 >  2:161715/1‑93 (MQ=255)

                                          gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCACCCAGCACAATATGAATCGCCAGTcc                                                >  2:6707861/1‑93 (MQ=255)

                                          gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                >  2:201662/1‑93 (MQ=255)

                                             tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                                             >  2:123345/1‑93 (MQ=255)

                                              gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                                            <  1:3798069/93‑1 (MQ=255)

                                                ccACTTCAACACTTCAAACGCAAGCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTg                                          <  1:4330427/93‑1 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  2:4191740/1‑93 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  1:3702281/1‑93 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  1:6053529/1‑93 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  2:7675329/1‑93 (MQ=255)

                                                   cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCa                                                    <  1:725364/80‑1 (MQ=255)

                                                   cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCa                                                    >  2:725364/1‑80 (MQ=255)

                                                   cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                                       >  1:1095420/1‑93 (MQ=255)

                                                   cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                                       >  1:7916791/1‑93 (MQ=255)

                                                   cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                                       >  1:3344017/1‑93 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      >  2:3046500/1‑93 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      >  2:7851209/1‑93 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      >  1:2700116/1‑93 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      <  1:1879443/93‑1 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAgac                                     >  2:1260282/1‑91 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcc                                     >  1:5296813/1‑93 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcc                                     <  2:2512282/93‑1 (MQ=255)

                                                      caacacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGccc                                    <  1:4503292/93‑1 (MQ=255)

                                                         cacttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGc                                 <  2:2621669/93‑1 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                <  2:2730464/93‑1 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               <  1:7842318/93‑1 (MQ=255)

                                                            ttcaaACGCAATCATTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCCGAAGCCCAGCAAt                              <  1:4758671/93‑1 (MQ=255)

                                                            ttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              >  2:4901988/1‑93 (MQ=255)

                                                            ttcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              >  1:1454773/1‑93 (MQ=255)

                                                             tcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATa                             <  2:5933727/93‑1 (MQ=255)

                                                             tcaaACGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATa                             >  2:5162768/1‑93 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  1:3176654/93‑1 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  2:7576101/93‑1 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     >  1:5513751/1‑93 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     <  1:983976/93‑1 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     >  1:7530366/1‑93 (MQ=255)

                                                                       tCACTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg                   >  2:3230460/1‑93 (MQ=255)

                                                                          cTGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGACAccgcgc                >  1:3635331/1‑93 (MQ=255)

                                                                           tGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgcc               <  2:5552327/93‑1 (MQ=255)

                                                                           tGAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgcc               <  1:2507175/93‑1 (MQ=255)

                                                                            gAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              >  1:4562261/1‑93 (MQ=255)

                                                                            gAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              >  2:318130/1‑93 (MQ=255)

                                                                            gAGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              <  1:6239219/93‑1 (MQ=255)

                                                                             aGCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgc             <  2:5383022/93‑1 (MQ=255)

                                                                              gCGGGAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgcg            >  2:289865/1‑93 (MQ=255)

                                                                                gggAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          <  1:3400067/93‑1 (MQ=255)

                                                                                gggAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          >  1:5171359/1‑93 (MQ=255)

                                                                                 ggAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATa         >  2:5830278/1‑93 (MQ=255)

                                                                                  gAAAATAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATaa        <  2:1333564/93‑1 (MQ=255)

                                                                                     aaTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCg     <  2:7365165/93‑1 (MQ=255)

                                                                                     aaTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCg     <  1:3558466/93‑1 (MQ=255)

                                                                                     aaTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCg     <  2:6243443/93‑1 (MQ=255)

                                                                                      aTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCgg    >  2:3929005/1‑92 (MQ=255)

                                                                                      aTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt    <  2:1506205/93‑1 (MQ=255)

                                                                                      aTAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt    <  1:12005/93‑1 (MQ=255)

                                                                                       tAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   <  2:6548187/93‑1 (MQ=255)

                                                                                       tAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   <  2:6166371/93‑1 (MQ=255)

                                                                                       tAGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   <  1:3431557/93‑1 (MQ=255)

                                                                                        aGCG‑TCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTgg  <  1:1274928/93‑1 (MQ=255)

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCG‑T‑CCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

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