BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,134,509	A→G	E274E (GAA→GAG)	b3943 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,134,509	0	A	G	96.2%	318.5 / 1.7	106	E274E (GAA→GAG)	b3943	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (4/0); new base G (55/47); total (59/47)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

Only 100 of 114 total aligned reads displayed.

CAGTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTT  >  NC_000913/4134419‑4134597

cAGTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAggg                                                                                        <  1:4732121/93‑1 (MQ=255)

cAGTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGCGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAggg                                                                                        >  1:340273/1‑93 (MQ=255)

 aGTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGt                                                                                       <  2:3902900/93‑1 (MQ=255)

  gTATTGCCAGCTTGAGCACGCTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTa                                                                                      <  1:2262519/93‑1 (MQ=255)

  gTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTa                                                                                      <  1:3270554/93‑1 (MQ=255)

    aTTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGc                                                                                    <  2:4515569/93‑1 (MQ=255)

    aTTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGc                                                                                    <  1:2989908/93‑1 (MQ=255)

     ttGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCg                                                                                   >  1:2872227/1‑93 (MQ=255)

         cAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTat                                                                               <  1:3263650/93‑1 (MQ=255)

          aGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTatt                                                                              >  1:3869769/1‑93 (MQ=255)

          aGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTatt                                                                              <  2:3274708/93‑1 (MQ=255)

          aGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTatt                                                                              >  1:2589097/1‑93 (MQ=255)

           gtttAAGCACACTGGGCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattg                                                                             <  2:1497622/91‑1 (MQ=255)

           gCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattg                                                                             <  1:1047654/93‑1 (MQ=255)

           gCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattg                                                                             <  2:128517/93‑1 (MQ=255)

           gCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattg                                                                             >  1:2202743/1‑93 (MQ=255)

           gCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGGAGCGGTattg                                                                             >  1:4750007/1‑93 (MQ=255)

             ttAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtg                                                                           >  1:3582254/1‑93 (MQ=255)

               aaGCACACTGGCCGTTCCCGTTTACGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgct                                                                         <  2:3226005/93‑1 (MQ=255)

                aGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctg                                                                        >  1:3455377/1‑93 (MQ=255)

                  cacaCTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgat                                                                      >  2:2652583/1‑93 (MQ=255)

                   acacTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTaatgtgctgatt                                                                     >  2:2210116/1‑93 (MQ=255)

                   acacTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCAGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgatt                                                                     >  1:2600725/1‑93 (MQ=255)

                     acTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgt                                                                   >  2:399318/1‑93 (MQ=255)

                     acTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgt                                                                   >  1:4260029/1‑93 (MQ=255)

                      cTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtt                                                                  >  2:1686033/1‑92 (MQ=255)

                      cTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtg                                                                  >  1:2076613/1‑93 (MQ=255)

                      cTGGCCGTTCCCGTTTGCGGCGTACTCGTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGCGGGTAGCGGTattgtgctgattgtg                                                                  <  2:659564/93‑1 (MQ=255)

                       tGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtgc                                                                 <  1:4144912/93‑1 (MQ=255)

                       tGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtgc                                                                 >  2:2477207/1‑93 (MQ=255)

                       tGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtgc                                                                 <  2:701780/93‑1 (MQ=255)

                        ggCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtgct                                                                <  2:4128309/93‑1 (MQ=255)

                        ggCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtgcg                                                                >  2:1961737/1‑92 (MQ=255)

                        ggCCGTTCCCGTCTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtgct                                                                >  1:4383884/1‑93 (MQ=255)

                         gCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtgctg                                                               <  2:2576991/93‑1 (MQ=255)

                           cGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGc                                                             >  2:1323003/1‑93 (MQ=255)

                           cGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGc                                                             <  2:1083740/93‑1 (MQ=255)

                            gTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCa                                                            >  1:1478665/1‑93 (MQ=255)

                            gTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGAGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCa                                                            >  1:231501/1‑93 (MQ=255)

                             ttCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGCGGGTAGCGGTATTGTGCTGATTGTGCTGGCAc                                                           >  1:1482875/1‑93 (MQ=255)

                             ttCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCAc                                                           >  1:2515808/1‑93 (MQ=255)

                             ttCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCAc                                                           >  2:1871921/1‑93 (MQ=255)

                             ttCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCAc                                                           >  2:558620/1‑93 (MQ=255)

                             ttCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCAc                                                           >  1:1256908/1‑93 (MQ=255)

                             gtCCCGTCTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATGGTGCTGATTGTGCTGGCAc                                                           >  2:4222761/2‑93 (MQ=255)

                              tCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACt                                                          <  2:3336475/93‑1 (MQ=255)

                              tCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACt                                                          <  1:3628035/93‑1 (MQ=255)

                              tCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACt                                                          <  1:881987/93‑1 (MQ=255)

                               cccTTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTg                                                         <  2:3339481/93‑1 (MQ=255)

                               cccGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTg                                                         >  2:656026/1‑93 (MQ=255)

                                ccGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTgg                                                        >  2:1810032/1‑93 (MQ=255)

                                  gTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGgcg                                                      >  2:680675/1‑93 (MQ=255)

                                   tttGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGgcgc                                                     <  1:693966/93‑1 (MQ=255)

                                   tttGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGgcgc                                                     >  2:3996795/1‑93 (MQ=255)

                                       cggcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGt                                                 >  2:1840684/1‑93 (MQ=255)

                                        ggcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTg                                                <  2:3074288/93‑1 (MQ=255)

                                         gcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGa                                               <  2:3437784/93‑1 (MQ=255)

                                         gcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTgg                                               >  2:4035060/1‑92 (MQ=255)

                                           gTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGc                                             <  1:2241219/93‑1 (MQ=255)

                                            tACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGcc                                            >  1:22937/1‑93 (MQ=255)

                                              ctctTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCga                                          >  1:520141/1‑92 (MQ=255)

                                               tctTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTa                                         <  1:1037864/93‑1 (MQ=255)

                                               tctTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTa                                         >  2:1278410/1‑93 (MQ=255)

                                                ctTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTaa                                        <  1:3662698/93‑1 (MQ=255)

                                                ctTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTaa                                        <  1:1600530/93‑1 (MQ=255)

                                                  tttCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGa                                      <  1:1055742/93‑1 (MQ=255)

                                                    tCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaa                                    >  1:3405031/1‑93 (MQ=255)

                                                    tCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaa                                    >  1:4446622/1‑93 (MQ=255)

                                                      cTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaaaa                                  <  1:2723294/93‑1 (MQ=255)

                                                      cTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaaaa                                  <  2:1071877/93‑1 (MQ=255)

                                                       tggtggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaaaaa                                 >  2:4576496/1‑93 (MQ=255)

                                                        ggtggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAaag                                <  1:815901/93‑1 (MQ=255)

                                                        ggtggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAaag                                <  1:540803/93‑1 (MQ=255)

                                                         gtggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAaaga                               >  2:392968/1‑93 (MQ=255)

                                                         gtggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAaaga                               <  1:3338992/93‑1 (MQ=255)

                                                          tggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAaagaa                              >  1:3852768/1‑93 (MQ=255)

                                                          tggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAaagaa                              <  2:4095587/93‑1 (MQ=255)

                                                           ggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAaagaag                             <  1:4442386/93‑1 (MQ=255)

                                                            gctgctCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGc                            >  2:2300524/1‑93 (MQ=255)

                                                            gctgctCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGc                            >  2:2770918/1‑93 (MQ=255)

                                                              tgctCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCg                          >  2:621979/1‑93 (MQ=255)

                                                              tgctCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCg                          >  1:1351112/1‑93 (MQ=255)

                                                               gctCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGt                         >  1:4743439/1‑93 (MQ=255)

                                                               gctCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGt                         >  1:4738789/1‑93 (MQ=255)

                                                                 tCGGCGAGAATCGGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCa                       <  1:2746083/93‑1 (MQ=255)

                                                                 tCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCa                       >  1:418239/1‑93 (MQ=255)

                                                                  cGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAg                      <  2:1581272/93‑1 (MQ=255)

                                                                      gagaATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTa                  >  1:1512370/1‑93 (MQ=255)

                                                                      gagaATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTa                  >  2:4157589/1‑93 (MQ=255)

                                                                       agaATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTaa                 <  1:3888999/93‑1 (MQ=255)

                                                                       agaATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTaa                 <  2:4588278/93‑1 (MQ=255)

                                                                          aTCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTaaaaa              <  2:1447963/93‑1 (MQ=255)

                                                                              gggggCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGAt          <  2:486003/93‑1 (MQ=255)

                                                                               ggggCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATc         <  1:440251/93‑1 (MQ=255)

                                                                                  gCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGTTCTGa      >  2:4004783/1‑93 (MQ=255)

                                                                                   gcGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGaa     <  1:3149322/92‑1 (MQ=255)

                                                                                   ccGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGaa     >  1:3999742/1‑93 (MQ=255)

                                                                                    cGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAAt    >  2:1176699/1‑93 (MQ=255)

                                                                                      ttGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAAttt  >  2:2333829/1‑93 (MQ=255)

                                                                                      ttGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAAttt  >  2:3825/1‑93 (MQ=255)

CAGTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTT  >  NC_000913/4134419‑4134597

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