BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,051,163	Δ1 bp	coding (94/231 nt)	sfa →	Sfa

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,051,162	0	A	.	86.8%	754.2 / 124.7	219	coding (93/231 nt)	sfa	Sfa
Reads supporting (aligned to +/- strand): ref base A (15/14); new base . (87/103); total (102/117)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 230 total aligned reads displayed.

ATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAC  >  NC_000913/1051070‑1051250

aTGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTa                                                                                          <  2:4989086/93‑1 (MQ=255)

aTGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTa                                                                                          >  1:2075644/1‑93 (MQ=255)

aTGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTa                                                                                          <  1:4820357/93‑1 (MQ=255)

 tGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTag                                                                                         >  2:5476118/1‑92 (MQ=255)

 tGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTag                                                                                         <  1:3040698/93‑2 (MQ=255)

 tGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTag                                                                                         <  1:5283647/93‑2 (MQ=255)

       gcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATaa                                                                                  >  2:5182055/1‑93 (MQ=255)

       gcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATaa                                                                                  >  2:580272/1‑93 (MQ=255)

       gcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATaa                                                                                  >  1:5446407/1‑93 (MQ=255)

            aTTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTg                                                                              <  1:6450556/93‑1 (MQ=255)

            aTTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGt                                                                             <  2:4188435/93‑1 (MQ=255)

               tCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTg                                                                           >  2:4607797/1‑93 (MQ=255)

                 acaAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGatt                                                                        <  1:1606508/93‑1 (MQ=255)

                  caAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGatta                                                                       >  2:742373/1‑93 (MQ=255)

                     aataatATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTg                                                                    >  1:5646144/1‑93 (MQ=255)

                     aataatATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTg                                                                    <  2:2128626/93‑1 (MQ=255)

                       taataTCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGaa                                                                  >  2:3277635/1‑93 (MQ=255)

                        aataTCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAAt                                                                 >  1:2156628/1‑93 (MQ=255)

                         atatCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATc                                                                <  1:4583075/93‑1 (MQ=255)

                            tCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCttt                                                             <  1:7054151/93‑1 (MQ=255)

                            tCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCttt                                                             >  1:7172919/1‑93 (MQ=255)

                              aGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCg                                                           >  2:2021972/1‑93 (MQ=255)

                              aGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCg                                                           <  2:351622/93‑1 (MQ=255)

                              aGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCg                                                           >  1:4696796/1‑93 (MQ=255)

                                gCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCggg                                                         <  2:153358/93‑1 (MQ=255)

                                gCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCggg                                                         <  2:5741447/93‑1 (MQ=255)

                                   cccTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGgag                                                      >  1:5014114/1‑93 (MQ=255)

                                     cTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGAt                                                    >  2:1294172/1‑93 (MQ=255)

                                     cTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGAt                                                    <  2:3100739/93‑1 (MQ=255)

                                      tCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGAt                                                    >  1:5142019/1‑93 (MQ=255)

                                      tCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATg                                                   <  1:1055652/93‑1 (MQ=255)

                                        gTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGc                                                 <  2:4039368/93‑1 (MQ=255)

                                          ggAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCt                                                >  1:7511756/1‑93 (MQ=255)

                                           gAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTta                                              <  1:384390/93‑1 (MQ=255)

                                           gAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTta                                              >  1:1654935/1‑93 (MQ=255)

                                               cTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAAc                                          <  1:5167963/93‑1 (MQ=255)

                                               cTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAAc                                          >  1:4369912/1‑93 (MQ=255)

                                                ttttttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACa                                         >  2:728846/1‑93 (MQ=255)

                                                 tttttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACAt                                        >  2:6525283/1‑93 (MQ=255)

                                                 tttttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACAt                                        >  1:4616055/1‑93 (MQ=255)

                                                  ttttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACAtt                                       <  2:6551116/93‑1 (MQ=255)

                                                   tttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACAtt                                       <  1:1807342/93‑1 (MQ=255)

                                                   tttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACAttt                                      <  2:6329905/93‑1 (MQ=255)

                                                   tttATATCTGCCTTATCTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACAtt                                       <  1:3345616/93‑1 (MQ=255)

                                                    ttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTc                                     <  2:2449017/93‑1 (MQ=255)

                                                     tataTCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCt                                    >  1:6704226/1‑93 (MQ=255)

                                                        atCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGGGATGGCTTATAACATTTCTTAc                                 >  1:2707111/1‑93 (MQ=255)

                                                        atCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTAc                                 <  1:4016976/93‑1 (MQ=255)

                                                         tCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTAcc                                >  1:6208754/1‑93 (MQ=255)

                                                           tGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTg                              >  1:6503502/1‑93 (MQ=255)

                                                            gCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTg                              >  2:5682611/1‑93 (MQ=255)

                                                            gCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGa                             >  2:5912383/1‑93 (MQ=255)

                                                            gCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGa                             <  2:2070896/93‑1 (MQ=255)

                                                              cTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGAcc                           >  1:739870/1‑93 (MQ=255)

                                                                tATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAg                         >  2:2289129/1‑93 (MQ=255)

                                                                tATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAg                         <  2:7297177/93‑1 (MQ=255)

                                                                tATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAg                         <  2:5048972/93‑1 (MQ=255)

                                                                 aTTTTTCTTCAACGCTGTATGTATAGT‑AGCTATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTAACCTGACCAgg                        <  2:6796055/93‑1 (MQ=255)

                                                                 aTTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAgg                        <  1:1850803/93‑1 (MQ=255)

                                                                 aTTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAgg                        <  2:784080/93‑1 (MQ=255)

                                                                  tttttCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGTTTTTTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAggg                       >  2:977906/1‑93 (MQ=255)

                                                                  tttttCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAggg                       <  1:2811232/93‑1 (MQ=255)

                                                                    tttCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGt                      <  1:1089898/93‑1 (MQ=255)

                                                                     ttcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAc                    >  1:1343362/1‑93 (MQ=255)

                                                                     ttcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAc                    <  2:2912152/93‑1 (MQ=255)

                                                                     ttcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAc                    <  2:6154066/93‑1 (MQ=255)

                                                                     ttcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAc                    <  1:3561846/93‑1 (MQ=255)

                                                                      tcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAcc                   <  1:635363/93‑1 (MQ=255)

                                                                      tcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAcc                   >  1:324296/1‑93 (MQ=255)

                                                                      tcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAcc                   >  1:5878218/1‑93 (MQ=255)

                                                                       cttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCg                  >  2:700024/1‑93 (MQ=255)

                                                                       cttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCg                  <  1:572673/93‑1 (MQ=255)

                                                                        ttcAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCg                  >  2:1651722/1‑93 (MQ=255)

                                                                            aCGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGaa              >  2:5344535/1‑93 (MQ=255)

                                                                            aCGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAAc             >  1:4209671/1‑93 (MQ=255)

                                                                             cGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAAcc            >  1:2961629/1‑93 (MQ=255)

                                                                              gCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCa           >  1:1183426/1‑93 (MQ=255)

                                                                              gCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCa           >  1:6811746/1‑93 (MQ=255)

                                                                                tgtatgtaTAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAac         >  2:1014271/1‑93 (MQ=255)

                                                                                tgtatgtaTAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAac         >  2:440683/1‑93 (MQ=255)

                                                                                tgtatgtaTAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAac         <  1:1551415/93‑1 (MQ=255)

                                                                                tgtatgtaTAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAac         >  2:3427969/1‑93 (MQ=255)

                                                                                tgtatgtaTAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAac         >  2:3718259/1‑93 (MQ=255)

                                                                                 gtatgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAaca        >  1:3533536/1‑93 (MQ=255)

                                                                                 gtatgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAaca        <  1:6640551/93‑1 (MQ=255)

                                                                                  tatgtatAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAaca        >  1:1282287/1‑93 (MQ=255)

                                                                                  tatgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAacac       <  1:6310398/93‑1 (MQ=255)

                                                                                  tatgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAacac       >  2:4672769/1‑93 (MQ=255)

                                                                                   atgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACAcc      >  1:3168690/1‑93 (MQ=255)

                                                                                   atgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACAcc      <  1:2403797/93‑1 (MQ=255)

                                                                                    tgtatAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACAcc      <  2:2837124/93‑1 (MQ=255)

                                                                                    tgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCt     >  2:6880219/1‑93 (MQ=255)

                                                                                    tgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCt     >  1:5227066/1‑93 (MQ=255)

                                                                                    tgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCt     >  2:7034402/1‑93 (MQ=255)

                                                                                     gtatAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCt     >  1:2603829/1‑93 (MQ=255)

                                                                                     gtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCtt    >  1:5811357/1‑93 (MQ=255)

                                                                                     gtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCtt    >  1:3474585/1‑93 (MQ=255)

                                                                                      tataGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTa   >  1:4102687/1‑93 (MQ=255)

                                                                                      tataGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTa   <  1:418668/93‑1 (MQ=255)

                                                                                       ataGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAc  <  2:7509732/93‑1 (MQ=255)

ATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAC  >  NC_000913/1051070‑1051250

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