BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,892,388	A→G	E36E (GAA→GAG)	yieF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,892,388	0	A	G	96.6%	787.8 / 10.9	264	E36E (GAA→GAG)	yieF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (4/5); new base G (124/131); total (128/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

Only 100 of 272 total aligned reads displayed.

AGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAAGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCT  >  NC_000913/3892297‑3892477

 ggtggtTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGt                                                                                         >  1:2530535/1‑93 (MQ=255)

         cgttacTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGt                                                                                 <  2:1723126/93‑1 (MQ=255)

          gttacTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGtt                                                                                <  2:738353/93‑1 (MQ=255)

          gttacTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGtt                                                                                <  1:3268472/93‑1 (MQ=255)

          gttacTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGtt                                                                                >  1:717549/1‑93 (MQ=255)

           ttacTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTa                                                                               >  2:4444203/1‑93 (MQ=255)

           ttacTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTa                                                                               >  2:222212/1‑93 (MQ=255)

             acTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAACATTGCTCCGGCGAGCATGGAGGTCAATGCGTTAcc                                                                             >  1:2799067/1‑93 (MQ=255)

               tGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTAccat                                                                           >  1:4480259/1‑93 (MQ=255)

               tGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTAccat                                                                           <  1:745446/93‑1 (MQ=255)

                ggggAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTAccatc                                                                          >  2:152746/1‑93 (MQ=255)

                 gggAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTATCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTAccatcc                                                                         <  2:4621040/93‑1 (MQ=255)

                 gggAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTAccatcc                                                                         >  2:3778547/1‑93 (MQ=255)

                    aGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCAtt                                                                      <  2:5490/93‑1 (MQ=255)

                    aGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCAtt                                                                      <  1:3227641/93‑1 (MQ=255)

                    aGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCAtt                                                                      <  2:3873679/93‑1 (MQ=255)

                     gCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTg                                                                     >  2:3522982/1‑93 (MQ=255)

                       cTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGcc                                                                   >  2:1568420/1‑93 (MQ=255)

                       cTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGcc                                                                   >  1:1804967/1‑93 (MQ=255)

                       cTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAAGTCAATGCGTTACCATCCATTGcc                                                                   >  1:3054939/1‑93 (MQ=255)

                        tGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCg                                                                  >  1:2056854/1‑93 (MQ=255)

                         gcgcAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGa                                                                 >  2:3289272/1‑93 (MQ=255)

                           gcAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACa                                                               <  1:2099142/93‑1 (MQ=255)

                           gcAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACa                                                               >  1:4051561/1‑93 (MQ=255)

                           gcAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACa                                                               >  1:4632181/1‑93 (MQ=255)

                             aaaGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACAtt                                                             >  1:4035077/1‑93 (MQ=255)

                                 gCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCt                                                         <  1:3789809/93‑1 (MQ=255)

                                 gCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCt                                                         >  1:1600620/1‑93 (MQ=255)

                                  cTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCtt                                                        >  1:3007653/1‑93 (MQ=255)

                                     aTTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTa                                                     <  2:133752/93‑1 (MQ=255)

                                     aTTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTa                                                     <  2:3592252/93‑1 (MQ=255)

                                     aTTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTa                                                     <  1:30806/93‑1 (MQ=255)

                                      tttAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTAt                                                    >  1:3981690/1‑93 (MQ=255)

                                      tttAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTAt                                                    >  1:1728994/1‑93 (MQ=255)

                                      tttAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTAt                                                    >  1:996160/1‑93 (MQ=255)

                                       ttAATGGCATTGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATg                                                   <  1:2300708/93‑1 (MQ=255)

                                         aaTGGCATGGTTGCACGTCCCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGAc                                                 <  2:2353232/93‑1 (MQ=255)

                                         aaTGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGAc                                                 <  1:242323/93‑1 (MQ=255)

                                         aaTGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGAc                                                 >  1:1440760/1‑93 (MQ=255)

                                          aTGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACg                                                >  2:2255764/1‑93 (MQ=255)

                                           tGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGc                                               <  2:240279/93‑1 (MQ=255)

                                           tGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGc                                               <  1:3573287/93‑1 (MQ=255)

                                            ggCATGGTTGCACGTCCCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCt                                              <  1:51559/93‑1 (MQ=255)

                                             gCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTg                                             >  1:3490925/1‑93 (MQ=255)

                                             gCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTg                                             >  2:405137/1‑93 (MQ=255)

                                              cATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGa                                            >  2:2660957/1‑93 (MQ=255)

                                              cATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGa                                            >  1:3162420/1‑93 (MQ=255)

                                               aTGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGAc                                           <  2:4384275/93‑1 (MQ=255)

                                               aTGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGAc                                           <  2:616481/93‑1 (MQ=255)

                                               aTGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGAc                                           <  2:2300507/93‑1 (MQ=255)

                                               aTGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGAc                                           >  1:234040/1‑93 (MQ=255)

                                                tGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACg                                          >  2:2140893/1‑93 (MQ=255)

                                                  gTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTa                                        <  2:2961234/93‑1 (MQ=255)

                                                   ttGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTAc                                       <  1:4522747/93‑1 (MQ=255)

                                                    tGCACGTCCCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTAca                                      <  1:837020/93‑1 (MQ=255)

                                                    tGCACGTACCCTGCTGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTAca                                      <  2:4218309/93‑1 (MQ=255)

                                                    tGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTAca                                      <  2:2516090/93‑1 (MQ=255)

                                                    tGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTAca                                      <  1:4409218/93‑1 (MQ=255)

                                                    tGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTAca                                      <  1:333419/93‑1 (MQ=255)

                                                      cACGTACCCTGCCGAAAATTGCTCCCGCGAGCATGGAGGGCAATGCGTTACCATCGATTGCCGAGCTGCCCTTGTATGACGCTGGAGTAcagc                                    >  2:2050478/1‑93 (MQ=255)

                                                         gTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAgg                                 >  2:3441225/1‑93 (MQ=255)

                                                         gTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAgg                                 >  2:3653742/1‑93 (MQ=255)

                                                         gTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAgg                                 >  2:3151296/1‑93 (MQ=255)

                                                           aCCCTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGgaa                               <  2:3568201/93‑1 (MQ=255)

                                                             ccTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACATCAGgaaga                             >  2:304415/1‑93 (MQ=255)

                                                              cTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGgaagaa                            >  1:1001491/1‑93 (MQ=255)

                                                              cTGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGgaagaa                            >  2:1081737/1‑93 (MQ=255)

                                                               tGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGaagaag                           >  1:153175/1‑93 (MQ=255)

                                                               tGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGaagaag                           >  2:2110560/1‑93 (MQ=255)

                                                               tGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGaagaag                           <  2:262438/93‑1 (MQ=255)

                                                               tGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGaagaag                           <  2:2874425/93‑1 (MQ=255)

                                                               tGCCGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGaagaag                           >  2:104017/1‑93 (MQ=255)

                                                                gCCGAAAATTGCTCCGCCGAGCATGGATGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTCACGTACAGCAGGAAGAAgg                          <  2:4584796/93‑1 (MQ=255)

                                                                 ccGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGt                         <  2:1457047/93‑1 (MQ=255)

                                                                 ccGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGt                         <  2:2012106/93‑1 (MQ=255)

                                                                  cGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCGGGAAGAAGGtt                        <  1:1300385/93‑1 (MQ=255)

                                                                  cGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGtt                        <  1:4157113/93‑1 (MQ=255)

                                                                  cGAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGtt                        <  2:1372632/93‑1 (MQ=255)

                                                                   gAAAATTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGttt                       >  1:1044484/1‑93 (MQ=255)

                                                                    aaaaTTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGtttt                      >  1:1163602/1‑93 (MQ=255)

                                                                     aaaTTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTc                     <  2:3334379/93‑1 (MQ=255)

                                                                     aaaTTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTc                     <  2:1751130/93‑1 (MQ=255)

                                                                       aTTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCa                   <  1:4603668/93‑1 (MQ=255)

                                                                       aTTGCTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCa                   <  1:2962396/93‑1 (MQ=255)

                                                                           cTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCaa               >  2:1690297/1‑93 (MQ=255)

                                                                           cTCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCaa               <  2:2434938/93‑1 (MQ=255)

                                                                            tCCGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAAc              <  2:514736/93‑1 (MQ=255)

                                                                              cGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACgg            <  2:4307083/93‑1 (MQ=255)

                                                                              cGGCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACgg            <  2:1134073/93‑1 (MQ=255)

                                                                               ggCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGt           <  1:2281736/93‑1 (MQ=255)

                                                                                gCGAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGtt          >  2:3787154/1‑93 (MQ=255)

                                                                                  gAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGa        <  1:4424386/93‑1 (MQ=255)

                                                                                  gAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGa        >  2:1624826/1‑93 (MQ=255)

                                                                                  gAGCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGa        >  1:2700504/1‑93 (MQ=255)

                                                                                    gCATGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAg      <  1:1378332/93‑1 (MQ=255)

                                                                                       tGGAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGctc   >  1:3699842/1‑93 (MQ=255)

                                                                                       tGGAAGTCAATGAGATACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAACTTTTTCCAGCAACGGTTGAAGctc   <  2:1902872/93‑1 (MQ=255)

                                                                                        ggAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGctct  >  1:658156/1‑93 (MQ=255)

                                                                                        ggAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGctct  <  2:4421081/93‑1 (MQ=255)

                                                                                        ggAGGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGctct  <  2:286481/93‑1 (MQ=255)

AGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAAGTCAATGCGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCT  >  NC_000913/3892297‑3892477

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