BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,208,505	G→A	R462R (CGG→CGA)	b2122 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,208,505	0	G	A	99.6%	823.5 / ‑5.2	231	R462R (CGG→CGA)	b2122	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (107/123); total (107/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

Only 100 of 238 total aligned reads displayed.

TGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGATGC  >  NC_000913/2208415‑2208594

 gCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACct                                                                                        >  2:1181790/1‑93 (MQ=255)

 gCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACct                                                                                        >  1:6539447/1‑93 (MQ=255)

 gCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACct                                                                                        <  2:324525/93‑1 (MQ=255)

  ccTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACctc                                                                                       <  2:6703797/93‑1 (MQ=255)

  ccTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACctc                                                                                       <  2:4793528/93‑1 (MQ=255)

  ccTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACctc                                                                                       >  1:2335293/1‑93 (MQ=255)

   cTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACctct                                                                                      <  2:4505864/93‑1 (MQ=255)

   cTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACctct                                                                                      >  1:916329/1‑93 (MQ=255)

         aaTGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCa                                                                                <  2:2688445/93‑1 (MQ=255)

          aTGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAg                                                                               <  1:2102172/93‑1 (MQ=255)

           tGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAgg                                                                              <  2:6582127/93‑1 (MQ=255)

           tGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAgg                                                                              >  2:1759427/1‑93 (MQ=255)

            gtggAGCGCGCCTTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGt                                                                             <  1:1162064/90‑1 (MQ=255)

            gTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGt                                                                             >  2:5189863/1‑93 (MQ=255)

            gTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGt                                                                             <  1:1072022/93‑1 (MQ=255)

            gTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGt                                                                             <  1:3920095/93‑1 (MQ=255)

               gAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGa                                                                          <  2:2301153/93‑1 (MQ=255)

                aGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGAt                                                                         >  2:3699227/1‑93 (MQ=255)

                aGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGAt                                                                         <  1:3314785/93‑1 (MQ=255)

                 gcgcgcCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATa                                                                        <  2:5493650/93‑1 (MQ=255)

                   gcgcCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAAc                                                                      >  2:5550990/1‑93 (MQ=255)

                   gcgcCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAAc                                                                      <  1:1548170/93‑1 (MQ=255)

                   gcgcCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAAc                                                                      <  2:2151069/93‑1 (MQ=255)

                    cgccgcTACTCCAGCCGGGCATTGTGGCTTGGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACt                                                                     <  1:639652/93‑1 (MQ=255)

                    cgcCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACt                                                                     <  1:3512840/93‑1 (MQ=255)

                      ctgTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACtta                                                                   <  2:2057729/91‑1 (MQ=255)

                      ccGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACtta                                                                   <  2:5864842/93‑1 (MQ=255)

                       cGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttat                                                                  <  2:4296342/93‑1 (MQ=255)

                        gTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttatt                                                                 >  2:5698444/1‑93 (MQ=255)

                        gTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttatt                                                                 <  2:5601453/93‑1 (MQ=255)

                        gTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttatt                                                                 <  2:2761369/93‑1 (MQ=255)

                        gTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttatt                                                                 <  2:3201192/93‑1 (MQ=255)

                        gTTACGCCACGCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATAGGACCTCTTAATCAGGTCGATAACttatt                                                                 <  2:309344/93‑1 (MQ=255)

                          tACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttattta                                                               <  1:5457288/93‑1 (MQ=255)

                           cccccAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttatttat                                                              <  1:4764290/90‑1 (MQ=255)

                           aCGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACtaatttat                                                              >  1:3817199/1‑93 (MQ=255)

                            cGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTtatttatt                                                             >  2:4809497/1‑93 (MQ=255)

                            cGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTtatttatt                                                             <  1:401555/93‑1 (MQ=255)

                                aGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGc                                                         <  1:32031/93‑1 (MQ=255)

                                 gCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGcc                                                        >  1:5640683/1‑93 (MQ=255)

                                 gCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGcc                                                        >  2:6461980/1‑93 (MQ=255)

                                 gCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGcc                                                        >  1:3274070/1‑93 (MQ=255)

                                 gCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCACCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGcc                                                        >  1:596401/1‑93 (MQ=255)

                                   cGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCgg                                                      <  1:460100/93‑1 (MQ=255)

                                    gggCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGt                                                     >  1:5868225/1‑93 (MQ=255)

                                     ggagttGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTc                                                    <  1:6144420/89‑1 (MQ=255)

                                     ggCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTc                                                    <  1:4849438/93‑1 (MQ=255)

                                       cATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGc                                                  >  2:4515686/1‑93 (MQ=255)

                                       cATTGTGGCTTTGCGGAAATATTTGCGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGc                                                  <  2:433921/93‑1 (MQ=255)

                                        aTTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCt                                                 >  1:6617563/1‑93 (MQ=255)

                                        aTTGTGGCTTTGCGGGAATATTTACGCCAGCGCCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCt                                                 <  1:4588225/93‑1 (MQ=255)

                                         ttGTGGCTTTGCGGGATTATTTACGACAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTg                                                <  2:200109/93‑1 (MQ=255)

                                         ttGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTg                                                >  2:1512452/1‑93 (MQ=255)

                                           gtgGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAg                                              <  2:2357671/93‑1 (MQ=255)

                                             ggCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAgtg                                            <  2:1054158/93‑1 (MQ=255)

                                             ggCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAgtg                                            >  2:1636290/1‑93 (MQ=255)

                                               cTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTa                                          >  2:1835935/1‑93 (MQ=255)

                                               cTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTa                                          >  1:422040/1‑93 (MQ=255)

                                                tttGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTAt                                         <  2:3961146/93‑1 (MQ=255)

                                                 ttGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTAtt                                        <  2:814929/93‑1 (MQ=255)

                                                 ttGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTAtt                                        <  2:5055398/93‑1 (MQ=255)

                                                 ttGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTAtt                                        <  2:6389194/93‑1 (MQ=255)

                                                  tGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTAttt                                       <  1:1852751/93‑1 (MQ=255)

                                                   gCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTc                                      <  1:2634529/93‑1 (MQ=255)

                                                    cGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCg                                     <  1:1797531/93‑1 (MQ=255)

                                                     gggAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTNTTTCGc                                    >  2:5686189/1‑93 (MQ=255)

                                                      ggAATATTTACGCCTGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCt                                   >  1:4068065/1‑93 (MQ=255)

                                                       gAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTTGCTGAGTGTATTTCGCTc                                  <  1:3130951/93‑1 (MQ=255)

                                                         atatTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCgg                                >  1:4148291/1‑93 (MQ=255)

                                                          tatTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGt                               <  1:5913636/93‑1 (MQ=255)

                                                          tatTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGt                               <  2:5177885/93‑1 (MQ=255)

                                                          tatTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGt                               >  1:1447527/1‑93 (MQ=255)

                                                            tttACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTg                             >  1:1507273/1‑93 (MQ=255)

                                                                cGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAc                         >  1:2647407/1‑93 (MQ=255)

                                                                  cccgcgcCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcag                       <  1:3460528/90‑1 (MQ=255)

                                                                   cAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcagc                      >  1:2199105/1‑93 (MQ=255)

                                                                    aGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcagca                     <  2:2492680/93‑1 (MQ=255)

                                                                    aGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcagca                     >  1:5034189/1‑93 (MQ=255)

                                                                      gctCCCCCCGCCTGCATACGTCCTCGTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACagcagc                   <  1:3006911/91‑1 (MQ=255)

                                                                        tCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACagcagcgg                 >  1:5418903/1‑91 (MQ=255)

                                                                        tCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgccg                 >  2:384166/1‑93 (MQ=255)

                                                                         ccccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgccgc                >  2:461306/1‑93 (MQ=255)

                                                                          cccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgccgcc               >  2:1367252/1‑93 (MQ=255)

                                                                          cccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgccgcc               <  2:3030883/93‑1 (MQ=255)

                                                                          cccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgccgcc               <  1:6659185/93‑1 (MQ=255)

                                                                           ccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCa              >  2:4428488/1‑93 (MQ=255)

                                                                           ccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCa              >  1:2988106/1‑93 (MQ=255)

                                                                           ccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCa              >  1:827246/1‑93 (MQ=255)

                                                                           ccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCa              <  1:4393597/93‑1 (MQ=255)

                                                                              cGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGAc           <  1:4034977/93‑1 (MQ=255)

                                                                               gCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACg          >  2:6489093/1‑93 (MQ=255)

                                                                                ccTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGc         >  2:3347960/1‑93 (MQ=255)

                                                                                ccTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGc         >  1:1602168/1‑93 (MQ=255)

                                                                                ccTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGc         >  1:5735504/1‑93 (MQ=255)

                                                                                   gCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTgg      <  2:3464125/93‑1 (MQ=255)

                                                                                     aTACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGAt    <  2:6030484/93‑1 (MQ=255)

                                                                                     aTACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGAt    >  2:2945062/1‑93 (MQ=255)

                                                                                     aTACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGAt    <  1:6269703/93‑1 (MQ=255)

                                                                                      tACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGATg   >  2:1447747/1‑93 (MQ=255)

                                                                                       aCGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGATgc  <  1:2418668/93‑1 (MQ=255)

TGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGATGC  >  NC_000913/2208415‑2208594

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