BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,219,434	A→G	A211A (GCT→GCC)	arp ←	Arp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,219,434	0	A	G	99.7%	1193.9 / ‑5.2	376	A211A (GCT→GCC)	arp	Arp
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (183/192); total (184/192)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 405 total aligned reads displayed.

CACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGC  >  NC_000913/4219344‑4219523

caccaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAggcggc                                                                                         <  2:1437598/93‑1 (MQ=255)

   caGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGgcaac                                                                                      <  1:7332235/93‑1 (MQ=255)

     gTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGgcaacgc                                                                                    >  2:1188882/1‑93 (MQ=255)

     gTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAgcaacgc                                                                                    >  1:128743/1‑93 (MQ=255)

      tAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGgcaacgca                                                                                   >  2:3301921/1‑93 (MQ=255)

         tatGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGcaacgcaacg                                                                                >  1:7765390/1‑93 (MQ=255)

             gCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGcccc                                                                            <  2:4501084/93‑1 (MQ=255)

              cTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCt                                                                           <  1:5617834/93‑1 (MQ=255)

               tGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTc                                                                          >  1:1203593/1‑93 (MQ=255)

               tGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTc                                                                          >  1:6397649/1‑93 (MQ=255)

                gCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCa                                                                         >  2:6415926/1‑93 (MQ=255)

                gCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCa                                                                         >  2:6247732/1‑93 (MQ=255)

                  cATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAgg                                                                       <  2:896613/93‑1 (MQ=255)

                  cATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAgg                                                                       <  1:517044/93‑1 (MQ=255)

                   aTAACACTGATGCCATTTGCTTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCCACGCCCCTCAGGa                                                                      >  1:1009969/1‑93 (MQ=255)

                   aTAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGa                                                                      >  2:5449929/1‑93 (MQ=255)

                   aTAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGa                                                                      <  2:7579887/93‑1 (MQ=255)

                   aTAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGa                                                                      <  1:2015476/93‑1 (MQ=255)

                    tAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGAc                                                                     <  2:3062315/93‑1 (MQ=255)

                       cacTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGc                                                                  <  1:7219430/93‑1 (MQ=255)

                        acTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCt                                                                 <  1:3298045/93‑1 (MQ=255)

                        acTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCt                                                                 <  2:3820670/93‑1 (MQ=255)

                         cTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTg                                                                <  1:110141/93‑1 (MQ=255)

                         cTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTg                                                                >  1:2468295/1‑93 (MQ=255)

                           gATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGa                                                              >  1:6601862/1‑93 (MQ=255)

                           gATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGa                                                              <  1:577074/93‑1 (MQ=255)

                           gATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGa                                                              >  1:4860669/1‑93 (MQ=255)

                             tGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGaaa                                                            <  1:4269435/93‑1 (MQ=255)

                                cATTTGATTCCATCTTTTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACAGCTGGAAACCa                                                         >  1:4063245/1‑93 (MQ=255)

                                 aTTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCaa                                                        <  2:742738/93‑1 (MQ=255)

                                 aTTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCaa                                                        <  2:4005453/93‑1 (MQ=255)

                                 aTTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCaa                                                        <  2:6744577/93‑1 (MQ=255)

                                  tttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAAc                                                       <  1:6555777/93‑1 (MQ=255)

                                  tttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAAc                                                       <  1:6894504/93‑1 (MQ=255)

                                  tttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAAc                                                       >  1:3858651/1‑93 (MQ=255)

                                  tttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAAc                                                       <  1:2381591/93‑1 (MQ=255)

                                  tttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAAc                                                       <  2:7683124/93‑1 (MQ=255)

                                     gATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGcc                                                    >  2:1545704/1‑93 (MQ=255)

                                     gATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGcc                                                    >  2:5589749/1‑93 (MQ=255)

                                      aTTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCa                                                   >  2:4539908/1‑93 (MQ=255)

                                      aTTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCa                                                   <  1:1022323/93‑1 (MQ=255)

                                      aTTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCa                                                   >  1:1959540/1‑93 (MQ=255)

                                        tCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCaaa                                                 <  2:4951679/93‑1 (MQ=255)

                                        tCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCaaa                                                 <  2:6825429/93‑1 (MQ=255)

                                        tCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCaaa                                                 <  2:1931066/93‑1 (MQ=255)

                                         ccATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAg                                                <  1:3505295/93‑1 (MQ=255)

                                          cATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGt                                               >  2:5013852/1‑93 (MQ=255)

                                          cATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGt                                               >  1:2801454/1‑93 (MQ=255)

                                           aTCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTc                                              <  1:5841078/93‑1 (MQ=255)

                                            tCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCCCCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCt                                             <  2:7226308/93‑1 (MQ=255)

                                            tCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCt                                             <  2:3743313/93‑1 (MQ=255)

                                            gcTTGTCAAAATATGCTTTCATAACTAGCTCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCt                                             <  1:5402630/92‑1 (MQ=255)

                                              ttGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTgg                                           >  2:6565390/1‑93 (MQ=255)

                                                gTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGtt                                         <  1:7410188/93‑1 (MQ=255)

                                                gTCAAAATATGCTATAATATCTAGCCCGTAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCAGGAAGCCAACGCCAAAGTCTGGtt                                         <  1:7025921/93‑1 (MQ=255)

                                                 tCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGttt                                        <  1:5317487/93‑1 (MQ=255)

                                                 tCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGttt                                        >  2:281300/1‑93 (MQ=255)

                                                  cAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGtttt                                       <  1:938873/93‑1 (MQ=255)

                                                   aaaaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGttttt                                      <  2:3514045/93‑1 (MQ=255)

                                                      atatGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCa                                   >  2:7820172/1‑93 (MQ=255)

                                                       tatGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAt                                  <  2:2245224/93‑1 (MQ=255)

                                                        atGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtt                                 <  1:4675723/93‑1 (MQ=255)

                                                         tGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAttt                                <  2:1666199/93‑1 (MQ=255)

                                                           cTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAttttt                              >  1:7722382/1‑93 (MQ=255)

                                                            tttCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtttttt                             >  1:6774149/1‑93 (MQ=255)

                                                             ttCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttg                            <  1:4180435/93‑1 (MQ=255)

                                                              tCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgt                           >  1:3975635/1‑93 (MQ=255)

                                                              tCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgt                           <  1:2616859/93‑1 (MQ=255)

                                                              tCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgt                           >  1:3777721/1‑93 (MQ=255)

                                                              tCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgt                           <  2:4712473/93‑1 (MQ=255)

                                                                aTAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgttg                         <  1:5011909/93‑1 (MQ=255)

                                                                 tAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGa                        <  1:7625086/93‑1 (MQ=255)

                                                                 tAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGa                        <  1:3621384/93‑1 (MQ=255)

                                                                 tAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGa                        >  2:5756696/1‑93 (MQ=255)

                                                                  aaCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGaa                       <  2:3181284/93‑1 (MQ=255)

                                                                  aaCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGaa                       <  2:7378253/93‑1 (MQ=255)

                                                                  aaCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGaa                       <  2:6318289/93‑1 (MQ=255)

                                                                  aaCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGaa                       <  2:5313722/93‑1 (MQ=255)

                                                                   aCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAg                      >  1:1365311/1‑93 (MQ=255)

                                                                   aCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAg                      >  1:4361204/1‑93 (MQ=255)

                                                                   aCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAg                      <  2:5986696/93‑1 (MQ=255)

                                                                   aCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAg                      >  2:1570971/1‑93 (MQ=255)

                                                                    cTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGt                     >  2:4290297/1‑93 (MQ=255)

                                                                     tAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTg                    >  1:5932266/1‑93 (MQ=255)

                                                                      aGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGc                   >  2:3475772/1‑93 (MQ=255)

                                                                       gCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCg                  <  1:2098676/93‑1 (MQ=255)

                                                                       gCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCg                  <  2:2763592/93‑1 (MQ=255)

                                                                       gCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCg                  >  2:704250/1‑93 (MQ=255)

                                                                          cGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTga               <  1:859267/93‑1 (MQ=255)

                                                                            aaaTCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTgaga             >  2:2037625/1‑93 (MQ=255)

                                                                                cACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAtt         <  1:7513876/93‑1 (MQ=255)

                                                                                cACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAtt         <  1:1701431/93‑1 (MQ=255)

                                                                                 aCCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttt        >  2:2586662/1‑93 (MQ=255)

                                                                                  ccGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAtttt       <  1:3082979/93‑1 (MQ=255)

                                                                                   cGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttttt      >  2:4123619/1‑93 (MQ=255)

                                                                                    gcaggcggCAACGCAGCGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTc     <  2:2938604/93‑1 (MQ=255)

                                                                                    gcaggcggCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTc     >  1:5100576/1‑93 (MQ=255)

                                                                                    gcaggcggCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTc     >  2:4388994/1‑93 (MQ=255)

                                                                                    gcaggcggCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTc     >  1:1560677/1‑93 (MQ=255)

                                                                                       ggcggcAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGc  <  2:6097188/93‑1 (MQ=255)

CACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGC  >  NC_000913/4219344‑4219523

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