BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,945,896	T→C	V63A (GTA→GCA)	yifE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,945,896	0	T	C	99.7%	2112.8 / ‑4.4	658	V63A (GTA→GCA)	yifE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (342/314); total (344/314)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 683 total aligned reads displayed.

AACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACA  >  NC_000913/3945806‑3945985

  cTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGc                                                                                       >  1:4818848/1‑93 (MQ=255)

   tGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGcc                                                                                      >  2:1884102/1‑93 (MQ=255)

   tGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGcc                                                                                      >  1:5148090/1‑93 (MQ=255)

   tGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGcc                                                                                      >  1:2909841/1‑93 (MQ=255)

     cTTGAGCGTCGTGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGt                                                                                    <  1:2589569/93‑1 (MQ=255)

     cTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCgg                                                                                    >  2:303133/1‑92 (MQ=255)

     cTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGt                                                                                    >  1:1117676/1‑93 (MQ=255)

     cTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGt                                                                                    <  1:9696756/93‑1 (MQ=255)

     cTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGt                                                                                    >  2:7238776/1‑93 (MQ=255)

      ttGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTg                                                                                   >  2:318943/1‑93 (MQ=255)

       tGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTgg                                                                                  >  2:2366807/1‑93 (MQ=255)

           cGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGaa                                                                              <  1:9096067/93‑1 (MQ=255)

           cGTCAGTGTTATTCCTTCAATGAGTTTGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGaa                                                                              <  1:7701564/93‑1 (MQ=255)

            gTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAAc                                                                             <  1:2433600/93‑1 (MQ=255)

             tCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACg                                                                            >  2:4912117/1‑93 (MQ=255)

             tCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACg                                                                            >  2:3785805/1‑93 (MQ=255)

              cATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGt                                                                           >  2:5210361/1‑93 (MQ=255)

              cATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGt                                                                           >  2:4725800/1‑93 (MQ=255)

               aTGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTg                                                                          >  2:2783969/1‑93 (MQ=255)

               aTGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTg                                                                          <  2:434924/93‑1 (MQ=255)

                  gTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGc                                                                       <  1:5742142/93‑1 (MQ=255)

                   ttATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGcc                                                                      >  2:6375116/1‑93 (MQ=255)

                     aTGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAg                                                                    >  1:468923/1‑93 (MQ=255)

                     aTGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAg                                                                    >  2:7767128/1‑93 (MQ=255)

                      tGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGt                                                                   <  2:6212334/93‑1 (MQ=255)

                        ccTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGa                                                                 <  2:7029097/93‑1 (MQ=255)

                        ccTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGa                                                                 <  1:7029361/93‑1 (MQ=255)

                          ttCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACa                                                               >  2:6899941/1‑93 (MQ=255)

                            cAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGa                                                             >  2:4395008/1‑93 (MQ=255)

                             aaTGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGaa                                                            <  2:9663875/93‑1 (MQ=255)

                              aTGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAg                                                           <  2:8417485/93‑1 (MQ=255)

                              aTGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGCAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAg                                                           >  1:3865823/1‑93 (MQ=255)

                                gAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCa                                                         <  1:8896828/93‑1 (MQ=255)

                                 aGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAg                                                        <  1:6281862/93‑1 (MQ=255)

                                    tGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCTTAGCAGCATGCCGTGGCGAACGTGAGCTAGTGACAGAAGCAGAAc                                                     >  1:8981407/1‑93 (MQ=255)

                                      gATCTTGGCAAGCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcgc                                                   <  2:9008121/93‑1 (MQ=255)

                                      gATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcgg                                                   >  1:6651766/1‑92 (MQ=255)

                                       aTCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcgcg                                                  >  2:8771502/1‑93 (MQ=255)

                                       aTCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcgcg                                                  <  1:9475114/93‑1 (MQ=255)

                                       aTCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcgcg                                                  >  1:1204955/1‑93 (MQ=255)

                                        tCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgt                                                 <  1:2291841/93‑1 (MQ=255)

                                         cTTGTCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgtg                                                <  2:4696988/93‑1 (MQ=255)

                                         cTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgtg                                                <  2:1647441/93‑1 (MQ=255)

                                          ttGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgtgt                                               <  2:3454497/93‑1 (MQ=255)

                                           tGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGtgtg                                              >  1:5894177/1‑93 (MQ=255)

                                           tGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGtgtg                                              <  1:8761585/93‑1 (MQ=255)

                                            ggCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTgg                                             <  2:7563757/93‑1 (MQ=255)

                                            ggCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTgg                                             <  1:3946133/93‑1 (MQ=255)

                                              cAACCGCGAGCCGGTTCCCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAACCAGTGACAGAAGCAGAACGCGTGTGGTc                                           <  2:2189261/93‑1 (MQ=255)

                                              cAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTc                                           <  1:6403365/93‑1 (MQ=255)

                                               aaaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTcc                                          >  2:3903265/1‑93 (MQ=255)

                                                aaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCa                                         >  1:6201446/1‑93 (MQ=255)

                                                aaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCa                                         >  2:9221602/1‑93 (MQ=255)

                                                 aCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCaa                                        >  1:8208665/1‑93 (MQ=255)

                                                    cgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGta                                     <  2:6756638/93‑1 (MQ=255)

                                                    cgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGta                                     <  1:6849154/93‑1 (MQ=255)

                                                     gAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtat                                    >  2:6126665/1‑93 (MQ=255)

                                                     gAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtat                                    >  1:10034822/1‑93 (MQ=255)

                                                     gAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtat                                    >  2:4077792/1‑93 (MQ=255)

                                                     gAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtat                                    <  1:6531675/93‑1 (MQ=255)

                                                      aGCCGGTTACCGAGGAAGAGAACCTCGTCATAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtata                                   <  2:1014649/93‑1 (MQ=255)

                                                      aGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtata                                   <  1:7356990/93‑1 (MQ=255)

                                                      aGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtata                                   <  1:4456865/93‑1 (MQ=255)

                                                      aGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtata                                   >  2:3304282/1‑93 (MQ=255)

                                                       gCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTatat                                  <  1:1572743/93‑1 (MQ=255)

                                                       gCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTatat                                  >  2:5006722/1‑93 (MQ=255)

                                                        ccGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATg                                 <  1:7936401/93‑1 (MQ=255)

                                                         cGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGa                                <  2:9815087/93‑1 (MQ=255)

                                                         cGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGa                                <  1:9937852/93‑1 (MQ=255)

                                                          ggTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAc                               <  2:1743236/93‑1 (MQ=255)

                                                          ggTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAc                               >  2:6607042/1‑93 (MQ=255)

                                                           gTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcg                              >  1:4996693/1‑93 (MQ=255)

                                                           gTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcg                              >  2:5861735/1‑93 (MQ=255)

                                                            ttACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcgc                             <  1:4229569/93‑1 (MQ=255)

                                                             tACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcgcg                            >  2:1352352/1‑93 (MQ=255)

                                                              aCCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGt                           <  2:8173988/93‑1 (MQ=255)

                                                                 gAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTAtt                        <  1:8596042/93‑1 (MQ=255)

                                                                  aGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGTAGAACGCGTGCGGTTCAAGTATATGACGCGTATTa                       >  1:3498642/1‑93 (MQ=255)

                                                                  aGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTa                       <  1:3166307/93‑1 (MQ=255)

                                                                   ggAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTaa                      <  2:1460923/93‑1 (MQ=255)

                                                                    gAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAg                     <  1:1662076/93‑1 (MQ=255)

                                                                    gAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAg                     <  1:7809062/93‑1 (MQ=255)

                                                                     aaGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGc                    <  1:3094155/93‑1 (MQ=255)

                                                                      agagAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCg                   >  1:7548743/1‑93 (MQ=255)

                                                                      agagAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCg                   <  2:5490604/93‑1 (MQ=255)

                                                                       gagaAGCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGt                  <  1:6084477/93‑1 (MQ=255)

                                                                       gagaAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGt                  >  1:1337962/1‑93 (MQ=255)

                                                                        agaAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTCAGCGTc                 >  2:7937467/1‑93 (MQ=255)

                                                                          aaaCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCa               >  1:2025210/1‑93 (MQ=255)

                                                                           aaCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaa              <  1:5385189/93‑1 (MQ=255)

                                                                            aCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaaa             <  1:2164735/93‑1 (MQ=255)

                                                                                ttCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGt         >  2:5586731/1‑93 (MQ=255)

                                                                                   gTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtttt      <  2:2498415/93‑1 (MQ=255)

                                                                                   gTAGCAGCATGCCGTGGCGAACGTGAGCCAGCGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtttt      <  1:2771368/93‑1 (MQ=255)

                                                                                    tAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTc     >  1:3700790/1‑93 (MQ=255)

                                                                                    gagcagcATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTc     <  1:1811918/92‑1 (MQ=255)

                                                                                     agcagcATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTca    <  2:4232320/93‑1 (MQ=255)

                                                                                     agcagcATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTca    >  1:93520/1‑93 (MQ=255)

                                                                                     agcagcATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTca    <  2:9609379/93‑1 (MQ=255)

                                                                                      gcagcaTGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGATTTcac   >  1:5447975/1‑93 (MQ=255)

AACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACA  >  NC_000913/3945806‑3945985

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