BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,999,913	C→T	A516V (GCC→GTC)	b2866 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,999,913	0	C	T	98.7%	1057.5 / ‑2.0	304	A516V (GCC→GTC)	b2866	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/4); new base T (167/133); total (167/137)
Fisher's exact test for biased strand distribution p-value = 4.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 313 total aligned reads displayed.

GGCAGAAACCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGC  >  NC_000913/2999824‑3000002

ggCAGAAACCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGtctc                                                                                        <  1:1293385/93‑1 (MQ=255)

ggCAGAAACCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGtctc                                                                                        <  1:4970723/93‑1 (MQ=255)

ggCAGAAACCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGtctc                                                                                        <  1:3624380/93‑1 (MQ=255)

 gCAGAAGCCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCa                                                                                       >  2:7785287/1‑93 (MQ=255)

 gCAGAAACCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCa                                                                                       <  2:4050262/93‑1 (MQ=255)

 gCAGAAACCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCa                                                                                       >  2:5437517/1‑93 (MQ=255)

 gCAGAAACCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCa                                                                                       >  2:1519957/1‑93 (MQ=255)

       aCCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCag                                                                                 <  1:6686392/93‑1 (MQ=255)

       aCCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCag                                                                                 >  2:3401405/1‑93 (MQ=255)

       aCCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCag                                                                                 <  2:2284452/93‑1 (MQ=255)

           tGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCt                                                                             >  2:4664736/1‑93 (MQ=255)

            gggggTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTa                                                                            <  1:2468056/93‑1 (MQ=255)

             ggggTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTAt                                                                           <  1:578118/93‑1 (MQ=255)

              gggTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATg                                                                          <  2:1120868/93‑1 (MQ=255)

               ggTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGt                                                                         <  1:1546222/93‑1 (MQ=255)

                   ccGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGcc                                                                     >  2:1824671/1‑93 (MQ=255)

                   ccGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGcc                                                                     <  1:7586255/93‑1 (MQ=255)

                   ccGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGcc                                                                     >  2:2904952/1‑93 (MQ=255)

                   ccGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGcc                                                                     <  2:4310969/93‑1 (MQ=255)

                     ggTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCcgc                                                                   >  2:2633278/1‑93 (MQ=255)

                     ggTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCcgc                                                                   >  2:6315635/1‑93 (MQ=255)

                      gTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCcgcg                                                                  >  1:2140342/1‑93 (MQ=255)

                       tCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCgcgc                                                                 >  1:5514351/1‑93 (MQ=255)

                        cAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGcc                                                                >  1:3224835/1‑93 (MQ=255)

                         aGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCct                                                               >  2:1966106/1‑93 (MQ=255)

                         aGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCct                                                               >  1:5601636/1‑93 (MQ=255)

                            gACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCctgcg                                                            >  1:5546225/1‑93 (MQ=255)

                             aCGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCctgcgc                                                           <  1:3182863/93‑1 (MQ=255)

                              cGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCctgcgct                                                          <  2:6489884/93‑1 (MQ=255)

                               gTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCctgcgctg                                                         >  1:535937/1‑93 (MQ=255)

                                ttCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCctgcgctgc                                                        >  2:5749105/1‑93 (MQ=255)

                                ttCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCctgcgctgc                                                        <  1:6364157/93‑1 (MQ=255)

                                ttCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCctgcgctgc                                                        <  1:3442975/93‑1 (MQ=255)

                                ttCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCctgcgctgc                                                        >  1:7268896/1‑93 (MQ=255)

                                ttCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCACGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCctgcgctgc                                                        >  2:5800478/1‑93 (MQ=255)

                                gtCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTAGGTTGCCGCGCctgcgctgc                                                        <  2:5290161/92‑1 (MQ=255)

                                  cgcgTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCtgcgctgcgc                                                      <  1:6451483/93‑1 (MQ=255)

                                   gcgTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCa                                                     >  2:4093687/1‑93 (MQ=255)

                                    cgTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAg                                                    >  1:1672139/1‑93 (MQ=255)

                                     gTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGt                                                   >  2:7718118/1‑93 (MQ=255)

                                     gTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGt                                                   >  1:2744300/1‑93 (MQ=255)

                                     gTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGt                                                   >  1:1790137/1‑93 (MQ=255)

                                      ttATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTg                                                  <  2:7612684/93‑1 (MQ=255)

                                      ttATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTg                                                  >  2:6251482/1‑93 (MQ=255)

                                       tATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGc                                                 <  2:515971/93‑1 (MQ=255)

                                       tATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGc                                                 >  1:4753446/1‑93 (MQ=255)

                                       tATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGc                                                 <  2:2690727/93‑1 (MQ=255)

                                       tATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGc                                                 <  2:3175634/93‑1 (MQ=255)

                                       tATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGc                                                 <  2:3984525/93‑1 (MQ=255)

                                       tATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGc                                                 >  1:905099/1‑93 (MQ=255)

                                       tATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGc                                                 >  1:3003774/1‑93 (MQ=255)

                                           tCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCa                                             >  1:60740/1‑93 (MQ=255)

                                            cAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCAc                                            <  1:6817483/93‑1 (MQ=255)

                                              aCTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACta                                          >  1:5051372/1‑93 (MQ=255)

                                               cTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACtat                                         <  1:7894699/93‑1 (MQ=255)

                                               cTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACtat                                         <  2:6862038/93‑1 (MQ=255)

                                                tCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACtatt                                        <  2:4514090/93‑1 (MQ=255)

                                                tCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACtatt                                        >  1:5141132/1‑93 (MQ=255)

                                                tCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACtatt                                        >  2:4957543/1‑93 (MQ=255)

                                                 cAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACtatta                                       <  1:111287/93‑1 (MQ=255)

                                                  aaGATACGGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACtattat                                      <  2:7094916/93‑1 (MQ=255)

                                                  aaGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACtattat                                      >  2:5339072/1‑93 (MQ=255)

                                                  aaGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACtattat                                      >  2:2100594/1‑93 (MQ=255)

                                                  aaGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACtattat                                      >  1:7612930/1‑93 (MQ=255)

                                                    gATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTAttatta                                    >  1:5321209/1‑93 (MQ=255)

                                                      tACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTaaa                                  >  2:1253229/1‑93 (MQ=255)

                                                      tACCGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTaaa                                  <  1:7107656/93‑1 (MQ=255)

                                                         cGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAaga                               <  1:3787129/93‑1 (MQ=255)

                                                         cGACGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAaga                               >  2:4781183/1‑93 (MQ=255)

                                                           aCGTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAgaga                             >  1:6174653/1‑93 (MQ=255)

                                                             gTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGaaa                           >  2:3632378/1‑93 (MQ=255)

                                                             gTTACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGaaa                           <  1:972947/93‑1 (MQ=255)

                                                              ttACGCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGaaaa                          >  2:1614694/1‑93 (MQ=255)

                                                                 cGCCGTTCGATCCCGGCGCATTTGTCGCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAatca                       <  1:3904712/93‑1 (MQ=255)

                                                                  gCCGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAatcat                      <  1:2013202/93‑1 (MQ=255)

                                                                   ccGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAatcatc                     <  2:3825543/93‑1 (MQ=255)

                                                                   ccGTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAatcatc                     >  2:6039953/1‑93 (MQ=255)

                                                                     gTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGc                   >  2:7880748/1‑93 (MQ=255)

                                                                     gTTCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGc                   >  1:4027525/1‑93 (MQ=255)

                                                                      ttCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCt                  >  2:5650523/1‑93 (MQ=255)

                                                                      ttCGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCt                  >  1:2606066/1‑93 (MQ=255)

                                                                        cGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCa                <  2:803857/93‑1 (MQ=255)

                                                                        cGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCa                <  2:6607933/93‑1 (MQ=255)

                                                                        cGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCa                <  1:877112/93‑1 (MQ=255)

                                                                        cGATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCa                <  1:676323/93‑1 (MQ=255)

                                                                         gATCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCAc               <  2:2007315/93‑1 (MQ=255)

                                                                           tCCCGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCAcgc             >  2:3901116/1‑93 (MQ=255)

                                                                            cccGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCNCAcgcc            >  1:4291517/1‑93 (MQ=255)

                                                                             ccGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCAcgccg           >  1:547462/1‑93 (MQ=255)

                                                                             ccGACGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCAcgccg           >  2:7541041/1‑93 (MQ=255)

                                                                              cGGCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCAcgccgc          >  2:6153208/1‑93 (MQ=255)

                                                                               ggCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCa         >  2:7240882/1‑93 (MQ=255)

                                                                               ggCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCa         <  2:7415437/93‑1 (MQ=255)

                                                                               ggCGCATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCa         <  2:3247870/93‑1 (MQ=255)

                                                                                   cATTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCa     >  2:6753255/1‑93 (MQ=255)

                                                                                    aTTTGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCAt    >  2:6816535/1‑93 (MQ=255)

                                                                                     tttGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATg   >  2:5955931/1‑93 (MQ=255)

                                                                                     tttGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATg   <  2:5692480/93‑1 (MQ=255)

                                                                                     tttGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATg   >  1:5894100/1‑93 (MQ=255)

                                                                                      ttGTCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGc  <  1:4586261/93‑1 (MQ=255)

GGCAGAAACCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGC  >  NC_000913/2999824‑3000002

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