BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,019,057	A→G	V415V (GTA→GTG)	yigR →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,019,057	0	A	G	98.2%	677.3 / 0.3	220	V415V (GTA→GTG)	yigR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/4); new base G (102/114); total (102/118)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

Only 100 of 229 total aligned reads displayed.

AATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTG  >  NC_000913/4018968‑4019145

 aTCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGAc                                                                                      <  2:2025370/93‑1 (MQ=255)

  tCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACg                                                                                     >  1:1134567/1‑93 (MQ=255)

  tCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACg                                                                                     >  1:3226099/1‑93 (MQ=255)

    tGTTTAATCCGGCGCGGCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGcc                                                                                   <  1:675803/93‑1 (MQ=255)

     gTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCa                                                                                  >  1:398710/1‑93 (MQ=255)

       ttAATCCGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGc                                                                                <  2:1186786/93‑1 (MQ=255)

         aaTACGGCGCGTCGCTTCAATATGGAAGTGCATCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCtt                                                                              <  1:1347784/93‑1 (MQ=255)

         aaTACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCtt                                                                              <  1:4590128/93‑1 (MQ=255)

         aaTACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCtt                                                                              >  2:1385737/1‑93 (MQ=255)

         aaTACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCtt                                                                              >  2:1522144/1‑93 (MQ=255)

         aaTACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCtt                                                                              >  2:2873750/1‑93 (MQ=255)

           tACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTa                                                                            <  1:212029/93‑1 (MQ=255)

           tACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTa                                                                            >  2:3198406/1‑93 (MQ=255)

           tACGGCGCGGCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTTCTCTACGTCGAAGGGGTGGGACGCCAGCTTTa                                                                            <  2:1163588/93‑1 (MQ=255)

              ggCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATcc                                                                         >  1:4075367/1‑93 (MQ=255)

                 gcgTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCa                                                                      >  2:2262386/1‑93 (MQ=255)

                  cgTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCaa                                                                     >  1:3479878/1‑93 (MQ=255)

                     cGCTCCAATATGTAATTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTc                                                                  <  1:1917613/93‑1 (MQ=255)

                        ttCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGAt                                                               >  2:1527556/1‑93 (MQ=255)

                        ttCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGAt                                                               >  1:3652437/1‑93 (MQ=255)

                         tCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGAtt                                                              <  2:1853875/93‑1 (MQ=255)

                         tCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGAtt                                                              >  1:82206/1‑93 (MQ=255)

                          aaaTATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGAttt                                                             <  2:3512814/92‑1 (MQ=255)

                           aaTATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTa                                                            <  2:1333268/93‑1 (MQ=255)

                            atatGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTAt                                                           <  1:4039936/93‑1 (MQ=255)

                             tatGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATg                                                          <  2:2996572/93‑1 (MQ=255)

                             tatGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATg                                                          <  1:2531011/93‑1 (MQ=255)

                             tatGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATg                                                          >  2:1229926/1‑93 (MQ=255)

                               tGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGa                                                        <  1:1066550/93‑1 (MQ=255)

                                  aaGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGaaaa                                                     <  1:3403399/93‑1 (MQ=255)

                                   aGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAAc                                                    <  2:1246062/93‑1 (MQ=255)

                                   aGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAAc                                                    <  2:3552328/93‑1 (MQ=255)

                                      gCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGc                                                 <  2:323246/93‑1 (MQ=255)

                                      gCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGc                                                 <  1:3171146/93‑1 (MQ=255)

                                          ccGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAg                                             >  1:3784318/1‑93 (MQ=255)

                                           cGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAACACGGCGAAGc                                            >  2:4230295/1‑93 (MQ=255)

                                           cGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGc                                            >  1:2766374/1‑93 (MQ=255)

                                            gCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGcc                                           >  2:1055368/1‑93 (MQ=255)

                                            gCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGcc                                           >  2:1529720/1‑93 (MQ=255)

                                            gCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGcc                                           >  1:906579/1‑93 (MQ=255)

                                            gCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGcc                                           >  1:4233965/1‑93 (MQ=255)

                                             cAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCt                                          >  1:2450958/1‑93 (MQ=255)

                                             cAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCt                                          <  1:3934137/93‑1 (MQ=255)

                                               aCTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGTAAACGGCGAAGCCttt                                        >  1:4707385/1‑93 (MQ=255)

                                               aCTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCttt                                        >  1:2231653/1‑93 (MQ=255)

                                               aCTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCttt                                        <  2:365532/93‑1 (MQ=255)

                                                cTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTc                                       <  2:369421/93‑1 (MQ=255)

                                                cTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTc                                       >  1:1834212/1‑93 (MQ=255)

                                                cTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTc                                       <  2:1722470/93‑1 (MQ=255)

                                                 tGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTcc                                      <  2:4401054/93‑1 (MQ=255)

                                                 tGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACAGCGAAGCCTTTcc                                      >  1:4238997/1‑93 (MQ=255)

                                                 tGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGCGGGACGACAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTcc                                      >  2:2812835/1‑93 (MQ=255)

                                                  ggTGTTACTCCAGAAACCCCTGCTCTACGTCGAAGGGGGGGGACGCAAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCt                                     <  2:3449122/93‑1 (MQ=255)

                                                  ggTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGTAGCCTTTCCt                                     <  1:1429097/93‑1 (MQ=255)

                                                    tgtTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTgg                                   <  1:676327/93‑1 (MQ=255)

                                                     gtTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGa                                  >  2:453399/1‑93 (MQ=255)

                                                      ttACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAg                                 >  1:95754/1‑93 (MQ=255)

                                                      ttACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAg                                 <  2:2172467/93‑1 (MQ=255)

                                                        aCTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTc                               <  1:1755666/93‑1 (MQ=255)

                                                         cTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCg                              <  1:3849838/93‑1 (MQ=255)

                                                         cTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCg                              <  1:3271144/93‑1 (MQ=255)

                                                          tCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGt                             <  2:718710/93‑1 (MQ=255)

                                                             aGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGa                          >  2:2587422/1‑93 (MQ=255)

                                                             aGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGa                          <  2:2925860/93‑1 (MQ=255)

                                                              gAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGAt                         >  2:1635038/1‑93 (MQ=255)

                                                              gAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGAt                         >  2:3613786/1‑93 (MQ=255)

                                                              gAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGAt                         >  2:4463695/1‑93 (MQ=255)

                                                                aaaCCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTa                       <  1:3362428/93‑1 (MQ=255)

                                                                aaaCCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTa                       >  1:4603363/1‑93 (MQ=255)

                                                                 aaCCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTaa                      <  1:2972052/93‑1 (MQ=255)

                                                                 aaCCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTaa                      >  2:3931401/1‑93 (MQ=255)

                                                                  aCCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTaaa                     <  2:341329/93‑1 (MQ=255)

                                                                  aCCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTaaa                     <  2:3886628/93‑1 (MQ=255)

                                                                   cccTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAg                    <  2:3998971/93‑1 (MQ=255)

                                                                   cccTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAg                    >  1:3813807/1‑93 (MQ=255)

                                                                   cccTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGCATTAAAg                    <  2:4652773/93‑1 (MQ=255)

                                                                     cTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGAt                  >  1:4360126/1‑93 (MQ=255)

                                                                     cTGCTCTACGTCGAAGGGGCGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGGCGAGGATTAAAGac                  >  1:2891470/1‑92 (MQ=255)

                                                                      tGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATc                 >  1:1611697/1‑93 (MQ=255)

                                                                      tGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATc                 <  1:756573/93‑1 (MQ=255)

                                                                        ctctACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGTCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAg               <  2:828263/93‑1 (MQ=255)

                                                                        ctctACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAg               <  1:2195253/93‑1 (MQ=255)

                                                                           tACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTc            >  1:507438/1‑93 (MQ=255)

                                                                           tACGTCGAAGGGGTGGGACGCAAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTc            <  1:200364/93‑1 (MQ=255)

                                                                            aCGGCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAATCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCg           <  2:1560073/93‑1 (MQ=255)

                                                                              gTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGt         >  2:2536346/1‑93 (MQ=255)

                                                                              gTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGt         <  1:2804893/93‑1 (MQ=255)

                                                                              gTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGt         >  1:420593/1‑93 (MQ=255)

                                                                               tCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTa        >  1:3502515/1‑93 (MQ=255)

                                                                               tCGAAGGGGGGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTa        >  2:4420397/1‑93 (MQ=255)

                                                                                cGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTAt       <  2:2712015/93‑1 (MQ=255)

                                                                                 gAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTAtt      >  1:4542114/1‑93 (MQ=255)

                                                                                  aaGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTc     <  1:1904182/93‑1 (MQ=255)

                                                                                  aaGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTc     <  1:4524125/93‑1 (MQ=255)

                                                                                  aaGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTc     >  2:3105687/1‑93 (MQ=255)

                                                                                   aGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTcc    >  2:4724982/1‑93 (MQ=255)

                                                                                   aGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTcc    >  2:3276692/1‑93 (MQ=255)

                                                                                    ggggTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCt   >  1:3073199/1‑93 (MQ=255)

                                                                                    ggggTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCt   <  2:524209/93‑1 (MQ=255)

                                                                                     gggTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTg  >  1:4677874/1‑93 (MQ=255)

AATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTG  >  NC_000913/4018968‑4019145

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