BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,777	0	C	T	47.6%	102.5 / inf	1127	D67D (GAC→GAT)	ybcX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (330/260); new base T (248/287); total (579/548)
Fisher's exact test for biased strand distribution p-value = 1.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-05
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 1180 total aligned reads displayed.

GCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCG  >  NC_000913/580682‑580872

  tAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTa                                                                                               <  2:9816188/96‑1 (MQ=255)

   aaTTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTAc                                                                                              <  1:9978235/96‑1 (MQ=255)

   aaTTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTAc                                                                                              >  2:4116982/1‑96 (MQ=255)

     ttACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAt                                                                                            <  2:12187776/96‑1 (MQ=255)

      tACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACAtn                                                                                           <  1:1616580/96‑2 (MQ=255)

       aCTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATcc                                                                                          <  1:6975685/96‑1 (MQ=255)

          aTCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTg                                                                                       >  2:9215128/1‑96 (MQ=255)

          aTCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTg                                                                                       <  2:5298989/96‑1 (MQ=255)

          aTCTCCCCAAAGACTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGGTCATGAATGCGATGGGTAACGATTACATACGTg                                                                                       <  1:2897025/96‑1 (MQ=255)

                ccaaaGAGTCGGCCTACCAGGAAGGGCTCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGtca                                                                                 <  2:13440797/94‑3 (MQ=255)

                  aaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCaac                                                                               <  2:14015859/96‑2 (MQ=255)

                  aaaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAt                                                                               <  1:10001121/96‑1 (MQ=255)

                   aaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCaacg                                                                              <  2:6900503/96‑3 (MQ=255)

                   aaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCaacg                                                                              >  1:3175424/1‑94 (MQ=255)

                    agagTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAAcgt                                                                             >  1:8971853/1‑96 (MQ=255)

                    agagTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGGAACGATTACATACGTGAGGTCAAcgt                                                                             >  1:792168/1‑96 (MQ=255)

                      aaTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATgtgg                                                                           >  1:11573589/1‑96 (MQ=255)

                       aTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATgtggt                                                                          >  2:4552930/1‑96 (MQ=255)

                        tCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATgtggtg                                                                         >  1:9703380/1‑96 (MQ=255)

                         ccGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGa                                                                        >  1:5838167/1‑96 (MQ=255)

                          cGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGaa                                                                       <  2:240708/96‑1 (MQ=255)

                          cGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGaa                                                                       >  1:3830027/1‑96 (MQ=255)

                           gCATTCCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAg                                                                      <  2:1372691/96‑1 (MQ=255)

                            cctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGt                                                                     >  1:6300271/3‑96 (MQ=255)

                            cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGt                                                                     >  2:4290446/1‑96 (MQ=255)

                             aTACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTc                                                                    >  1:12414761/1‑96 (MQ=255)

                              tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCt                                                                   <  1:9144378/96‑1 (MQ=255)

                              tACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCt                                                                   <  1:2696224/96‑1 (MQ=255)

                              tACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCt                                                                   >  1:2097051/1‑96 (MQ=255)

                               aCCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTg                                                                  <  2:13083703/96‑1 (MQ=255)

                               aCCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTg                                                                  <  1:7386934/96‑1 (MQ=255)

                                ccAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGc                                                                 >  2:1075573/1‑96 (MQ=255)

                                  aGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGccc                                                               >  1:10251230/1‑96 (MQ=255)

                                   ggAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCg                                                              >  2:1420594/1‑96 (MQ=255)

                                    gAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgt                                                             <  2:634431/96‑1 (MQ=255)

                                    gAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgt                                                             >  1:12116568/1‑96 (MQ=255)

                                    gAAGGGCGCTGGAAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgt                                                             <  2:1523550/96‑1 (MQ=255)

                                     aaGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtg                                                            >  2:10304760/1‑96 (MQ=255)

                                      aGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgt                                                           <  2:1152873/96‑1 (MQ=255)

                                       gggCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgtn                                                          <  1:4245039/96‑2 (MQ=255)

                                        ggCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgttg                                                         <  2:4902749/96‑3 (MQ=255)

                                        ggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCg                                                         <  2:11359266/96‑1 (MQ=255)

                                         gcgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATCCGTGAGGTCAACGTTGTTAAGTCTGCCCgtgttgn                                                        <  1:5157700/96‑4 (MQ=255)

                                         ccgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTgg                                                        >  1:13105889/2‑96 (MQ=255)

                                           gcTGGGAAACGTTGCCGTTTCAACGCGCGGTCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggct                                                      >  2:11945799/1‑94 (MQ=255)

                                           gcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggct                                                      >  2:6164935/1‑94 (MQ=255)

                                           gcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGtt                                                      >  1:8331522/1‑96 (MQ=255)

                                            cagggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTa                                                     >  1:10834577/3‑96 (MQ=255)

                                            cTGNGAAACGTTGCCGTTTCAACGCGCGATCATTAATGCGATGGTTAACGCTTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggctn                                                     <  1:3656471/96‑4 (MQ=255)

                                             tGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTAt                                                    >  2:7281535/1‑96 (MQ=255)

                                             tGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTaa                                                    <  1:11334791/96‑2 (MQ=255)

                                              gggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTAtt                                                   >  2:13308199/1‑96 (MQ=255)

                                              gggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTAtt                                                   <  1:2300471/96‑1 (MQ=255)

                                              gggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTAtt                                                   <  2:8937866/96‑1 (MQ=255)

                                                gAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGGAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTca                                                 >  2:7745046/1‑95 (MQ=255)

                                                gAAACACTGACCATTCAGCTGGACATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTcc                                                 <  1:13601593/96‑1 (MQ=255)

                                                 aaaCACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGGGAAGTCTGCCCGTGTCGGTTATTCCa                                                >  2:12617670/1‑96 (MQ=255)

                                                    cgttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAaaaa                                             >  1:6673922/4‑96 (MQ=255)

                                                     gttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAt                                            <  1:2493475/94‑1 (MQ=255)

                                                     gttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAt                                            <  2:3482706/94‑1 (MQ=255)

                                                      ttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtg                                           <  2:2158163/95‑1 (MQ=255)

                                                      ttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtg                                           >  2:11142681/2‑96 (MQ=255)

                                                        gCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgtt                                         >  1:14336569/1‑94 (MQ=255)

                                                          cgtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTACTCAAAAAtgttgc                                       <  1:13894209/94‑1 (MQ=255)

                                                           gtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATgttgct                                      >  2:3829879/2‑96 (MQ=255)

                                                           cTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATgctgct                                      >  1:13350285/1‑96 (MQ=255)

                                                            tttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGttgctg                                     >  2:9463144/1‑96 (MQ=255)

                                                                 gCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtgt                                <  2:2033604/96‑1 (MQ=255)

                                                                  cGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGtt                               <  1:4530978/96‑1 (MQ=255)

                                                                    cgcgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                             >  1:12800447/5‑96 (MQ=255)

                                                                     gcgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTAt                            >  2:7501/4‑96 (MQ=255)

                                                                        atcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGcg                         >  2:8298084/1‑95 (MQ=255)

                                                                         tcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGcgt                        <  1:5421387/96‑3 (MQ=255)

                                                                         tcatGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCt                        >  1:14153623/1‑96 (MQ=255)

                                                                          catGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTa                       >  1:3463795/1‑96 (MQ=255)

                                                                            tGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTatt                     >  2:3075948/1‑94 (MQ=255)

                                                                            tGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACt                     >  1:9454001/1‑96 (MQ=255)

                                                                             gAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTAttt                    <  1:4750991/96‑1 (MQ=255)

                                                                              aaTGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTAtttt                   <  1:7208705/96‑1 (MQ=255)

                                                                              aaTGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTAtttt                   >  2:8659744/1‑96 (MQ=255)

                                                                              aaTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACttt                   <  2:5838167/96‑1 (MQ=255)

                                                                               aTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTta                  <  1:6068829/96‑1 (MQ=255)

                                                                                  cGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattc               >  2:9090287/1‑94 (MQ=255)

                                                                                  cGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATag               >  2:3254976/1‑96 (MQ=255)

                                                                                  cGATGGGCAGCGACTACATCCGTGAGGTGAAGGTGGTGAAGTCTGCCCGGGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATag               >  2:2874417/1‑96 (MQ=255)

                                                                                   gATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattca              <  2:65021/96‑4 (MQ=255)

                                                                                   gATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattca              <  1:5680284/96‑4 (MQ=255)

                                                                                    aTGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattcag             >  1:8754181/1‑92 (MQ=255)

                                                                                     tGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGc            >  2:2615519/1‑96 (MQ=255)

                                                                                       ggCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCAt          <  2:8068288/96‑1 (MQ=255)

                                                                                       ggCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCAt          >  1:189787/1‑96 (MQ=255)

                                                                                       gacAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGGGTTTATGCCTACTTTATAGAGCAt          >  1:10540186/3‑96 (MQ=255)

                                                                                          aGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAag       >  1:846799/1‑96 (MQ=255)

                                                                                           acGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagc      >  1:7680399/2‑96 (MQ=255)

                                                                                            cGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagca     <  2:10471720/96‑1 (MQ=255)

                                                                                            cGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAacca     >  1:11871784/1‑96 (MQ=255)

                                                                                               cTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAgcg  >  2:1750297/1‑96 (MQ=255)

                                                                                               cTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAgcg  >  2:5552018/1‑96 (MQ=255)

                                                                                               cTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAgcg  >  2:9699617/1‑96 (MQ=255)

                                                                                               cTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAgcg  >  1:4971093/1‑96 (MQ=255)

GCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCG  >  NC_000913/580682‑580872

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