BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,753	0	C	G	67.1%	1105.6 / inf	796	A59A (GCC→GCG)	ybcX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (146/115); new base G (222/311); total (369/427)
Fisher's exact test for biased strand distribution p-value = 2.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 854 total aligned reads displayed.

ACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTA  >  NC_000913/580661‑580845

   gcggTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCgatc                                                                                           >  1:788705/4‑93 (MQ=255)

    cggTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCgatca                                                                                          <  1:3442629/91‑1 (MQ=255)

      gTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATg                                                                                        >  2:4696032/1‑93 (MQ=255)

        tGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGaa                                                                                      <  1:9186/93‑1 (MQ=255)

        tGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGaa                                                                                      <  2:3265282/93‑1 (MQ=255)

         gAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAAt                                                                                     >  1:3799302/1‑93 (MQ=255)

         gAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAAt                                                                                     >  2:3759046/1‑93 (MQ=255)

         gAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGAGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAAt                                                                                     >  1:432036/1‑93 (MQ=255)

          aaTGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATg                                                                                    <  1:5060495/93‑1 (MQ=255)

          aaTGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATg                                                                                    <  2:4468875/93‑1 (MQ=255)

          aaTGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATg                                                                                    <  2:5973181/93‑1 (MQ=255)

             gggCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGa                                                                                 <  2:1448235/93‑1 (MQ=255)

             gggCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGa                                                                                 <  1:2340263/93‑1 (MQ=255)

             gggCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGa                                                                                 <  1:4772494/93‑1 (MQ=255)

                cGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATgg                                                                              <  1:1581088/93‑1 (MQ=255)

                   aggCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCa                                                                           <  1:4239219/91‑1 (MQ=255)

                       tAATTCCTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGGAACGa                                                                       >  1:1386599/1‑93 (MQ=255)

                       tAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGa                                                                       >  2:5443460/1‑93 (MQ=255)

                        aaTTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGat                                                                      <  1:5848200/93‑2 (MQ=255)

                        aaTTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGat                                                                      <  2:1676528/93‑2 (MQ=255)

                          ttACTACCCCCCCAAAGACGCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTa                                                                    <  1:153375/93‑1 (MQ=255)

                            aCTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACa                                                                  <  2:4651354/93‑1 (MQ=255)

                             cTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAt                                                                 >  2:4516794/1‑93 (MQ=255)

                             cTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAt                                                                 <  1:3628149/93‑1 (MQ=255)

                              tATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTAAATc                                                                >  1:1204035/1‑93 (MQ=255)

                              tATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAta                                                                <  2:1963484/93‑2 (MQ=255)

                              tATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAta                                                                <  2:4038286/93‑2 (MQ=255)

                               aTCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATcc                                                               <  2:28603/93‑1 (MQ=255)

                               aTCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAtac                                                               <  1:5949379/93‑3 (MQ=255)

                               aTCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAtac                                                               <  2:3473723/93‑3 (MQ=255)

                                tctcCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACg                                                              >  1:6098451/1‑93 (MQ=255)

                                 ctcCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGt                                                             <  2:2141978/93‑1 (MQ=255)

                                 ctcCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGt                                                             >  2:883758/1‑93 (MQ=255)

                                   ccgcaaaGAGTCGGCCTACCAGGAAGGGCGCTGCGAACCGTTGCCGTGTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGa                                                           <  1:725819/89‑1 (MQ=255)

                                     ccaaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAgg                                                         <  1:3062482/91‑1 (MQ=255)

                                     ccaaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAgg                                                         >  1:2487436/3‑93 (MQ=255)

                                     ccaaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATCGGTAACGATTACATACGTGAgg                                                         >  2:765152/3‑93 (MQ=255)

                                      gAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAggg                                                        >  1:1718258/1‑92 (MQ=255)

                                      caaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGt                                                        >  1:2873804/2‑93 (MQ=255)

                                        aaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGa                                                      >  1:5284896/1‑93 (MQ=255)

                                            aTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGGGAGGGGAAtgt                                                  >  1:2002464/1‑93 (MQ=255)

                                               cGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATgtggt                                               >  2:4166833/1‑93 (MQ=255)

                                                  ctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTaa                                            >  1:29382/2‑93 (MQ=255)

                                                   tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                           >  2:2430476/1‑93 (MQ=255)

                                                   tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                           >  2:2712256/1‑93 (MQ=255)

                                                   tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                           >  2:307677/1‑93 (MQ=255)

                                                    aCCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGt                                          >  2:262400/1‑93 (MQ=255)

                                                    aCCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGt                                          >  1:397842/1‑93 (MQ=255)

                                                     ccAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTc                                         <  1:4433234/93‑1 (MQ=255)

                                                     ccAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTc                                         <  1:690084/93‑1 (MQ=255)

                                                      cAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCt                                        >  2:5571588/1‑93 (MQ=255)

                                                      cAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATTGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCt                                        <  2:2276532/93‑1 (MQ=255)

                                                        ggAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGc                                      >  2:2992849/1‑93 (MQ=255)

                                                        ggAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGc                                      >  2:680840/1‑93 (MQ=255)

                                                         gAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcc                                     >  1:5624642/1‑93 (MQ=255)

                                                         gAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGcc                                     >  2:2510959/1‑93 (MQ=255)

                                                          aaGGGCGCTGGGTAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGccc                                    <  1:5479056/93‑1 (MQ=255)

                                                          aaGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGccc                                    <  2:4370447/93‑1 (MQ=255)

                                                          aaGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGccc                                    >  1:49985/1‑93 (MQ=255)

                                                           aGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCg                                   <  1:191396/93‑1 (MQ=255)

                                                            gggCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgt                                  <  2:5199252/93‑1 (MQ=255)

                                                            gggCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgt                                  <  2:3914320/93‑1 (MQ=255)

                                                             ggCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtg                                 >  2:1894214/1‑93 (MQ=255)

                                                              gcgcTGGGAAGCACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtgt                                <  1:3277818/93‑1 (MQ=255)

                                                              gcgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgt                                <  1:2613127/93‑1 (MQ=255)

                                                               cgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgtt                               <  2:5267187/93‑2 (MQ=255)

                                                                gcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCg                              <  2:5031466/93‑1 (MQ=255)

                                                                 cTGGGAAGCACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCgg                             <  1:4005437/93‑1 (MQ=255)

                                                                 cTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTgg                             <  2:5664371/93‑1 (MQ=255)

                                                                 cTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTgg                             <  2:4807059/93‑1 (MQ=255)

                                                                  tGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggc                            <  1:1081075/93‑2 (MQ=255)

                                                                   gggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggct                           >  1:38875/1‑91 (MQ=255)

                                                                   gggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggct                           <  1:4139693/93‑3 (MQ=255)

                                                                   gggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggct                           >  1:2164004/1‑91 (MQ=255)

                                                                    ggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTa                          >  2:4892418/1‑93 (MQ=255)

                                                                    ggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTa                          <  1:2788288/93‑1 (MQ=255)

                                                                    ggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTa                          <  2:6044473/93‑1 (MQ=255)

                                                                    ggAAACACTGCCCTTGCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTa                          <  2:535109/93‑1 (MQ=255)

                                                                     gAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTAt                         >  2:1769266/1‑93 (MQ=255)

                                                                     gAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTAt                         <  2:1335278/93‑1 (MQ=255)

                                                                     gAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTAt                         >  1:2721648/1‑93 (MQ=255)

                                                                     gAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTAt                         <  1:5785058/93‑1 (MQ=255)

                                                                       aaCACTGCCCTTTCAGCGGGCCGTCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTc                       >  2:1183188/1‑93 (MQ=255)

                                                                       aaCACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTc                       >  2:510855/1‑93 (MQ=255)

                                                                        acacTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTcc                      <  1:3190365/93‑1 (MQ=255)

                                                                          acTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaa                    >  1:4896842/1‑93 (MQ=255)

                                                                            tGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAaaaa                  <  2:1539278/93‑1 (MQ=255)

                                                                            tGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaaa                  >  2:334805/1‑93 (MQ=255)

                                                                            tGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaaa                  >  1:4668180/1‑93 (MQ=255)

                                                                             gCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAt                 <  1:2254008/93‑1 (MQ=255)

                                                                              ccgtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtg                <  1:5382829/90‑1 (MQ=255)

                                                                               cgtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgt               >  2:3190180/3‑92 (MQ=255)

                                                                                gtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgtt              <  2:1838523/92‑3 (MQ=255)

                                                                                gtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgtt              >  2:3828195/2‑91 (MQ=255)

                                                                                cTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgct              >  2:5818240/1‑93 (MQ=255)

                                                                                  ttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttgc            <  2:5858789/93‑1 (MQ=255)

                                                                                   tCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATgttgct           <  1:450599/93‑1 (MQ=255)

                                                                                     aGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTgg         >  2:3671271/1‑93 (MQ=255)

                                                                                      gCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTggg        >  2:1074748/1‑93 (MQ=255)

                                                                                      gCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTggg        <  1:5213052/93‑1 (MQ=255)

ACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTA  >  NC_000913/580661‑580845

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: `‑`