BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,074,827	A→G	Q78R (CAG→CGG)	yihZ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,074,827	0	A	G	99.7%	995.8 / NA	322	Q78R (CAG→CGG)	yihZ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (153/167); total (155/167)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 347 total aligned reads displayed.

GTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTTATATGACT  >  NC_000913/4074737‑4074916

gtgtGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGt                                                                                         >  2:1322430/1‑93 (MQ=255)

gtgtGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGt                                                                                         >  2:4015917/1‑93 (MQ=255)

 tgtgCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGtt                                                                                        >  1:6825438/1‑93 (MQ=255)

 tgtgCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCCGtt                                                                                        >  1:737263/1‑93 (MQ=255)

   tgCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTa                                                                                      >  2:6053292/1‑93 (MQ=255)

   tgCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTa                                                                                      >  1:302363/1‑93 (MQ=255)

    gCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTAc                                                                                     <  1:911189/93‑1 (MQ=255)

    gCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTAc                                                                                     <  1:5611709/93‑1 (MQ=255)

     cTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTAcc                                                                                    >  1:7753536/1‑93 (MQ=255)

      tCGTCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTAccc                                                                                   <  2:3885466/93‑1 (MQ=255)

      tCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTAccc                                                                                   <  2:2387075/93‑1 (MQ=255)

       cGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCt                                                                                  <  2:1251247/93‑1 (MQ=255)

       cGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCCGTTTACCCt                                                                                  >  2:4694027/1‑93 (MQ=255)

        ggCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTc                                                                                 <  1:1021366/93‑1 (MQ=255)

        ggCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGCTTACCCTc                                                                                 >  1:7154829/1‑93 (MQ=255)

          cTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTcgc                                                                               >  1:5001604/1‑93 (MQ=255)

           tACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTcgcc                                                                              >  1:7934018/1‑93 (MQ=255)

           tACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTcgcc                                                                              >  1:6310781/1‑93 (MQ=255)

           tACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTcgcc                                                                              >  1:5773433/1‑93 (MQ=255)

           tACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTcgcc                                                                              <  2:6423868/93‑1 (MQ=255)

           tACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGGTTCCCGGTTTACCCTcgcc                                                                              >  2:7222181/1‑93 (MQ=255)

             ccGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTcgccgc                                                                            >  1:254245/1‑93 (MQ=255)

              cGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCa                                                                           <  2:7800100/93‑1 (MQ=255)

               gCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAg                                                                          >  2:2022054/1‑93 (MQ=255)

               gCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAg                                                                          >  2:4830276/1‑93 (MQ=255)

               gCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAg                                                                          >  1:3273195/1‑93 (MQ=255)

                ctacTTTATGGATTCCGAAGGCAATATGAATCTCAGCGTGCAACAGTCGGGCGGCAGTGTGCTGGTGGTTTCCCGCTTTACCCTCGCCGCAGa                                                                         <  2:7619198/90‑1 (MQ=255)

                cATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGa                                                                         >  2:1605638/1‑93 (MQ=255)

                 aTCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGAt                                                                        >  2:2272548/1‑93 (MQ=255)

                 aTCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGAt                                                                        >  2:753320/1‑93 (MQ=255)

                  tCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACCGGCGGGCGGCAGTGTGCTGGTGGTTTCACGGTTTACCATCGCCGCAGAtc                                                                       >  2:6538943/1‑92 (MQ=255)

                  tCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTGCCGGTTTACCCTCGCCGCAGATa                                                                       >  2:4330360/1‑93 (MQ=255)

                  tCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGAtt                                                                       >  2:5929651/1‑92 (MQ=255)

                  tCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATa                                                                       >  2:2796971/1‑93 (MQ=255)

                  tCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATa                                                                       <  1:3386509/93‑1 (MQ=255)

                    tttAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATAcc                                                                     >  1:7052942/1‑93 (MQ=255)

                       aGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGaa                                                                  <  1:2362413/93‑1 (MQ=255)

                       aGCGAGGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGATGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGaa                                                                  <  2:2851030/93‑1 (MQ=255)

                        gCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTGTACCCTCGCCGCAGATCCCGGAc                                                                 >  1:726289/1‑93 (MQ=255)

                        gCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGCTTCCCGGTTTACCCTCGCCGCAGATACCGAAc                                                                 >  1:1226144/1‑93 (MQ=255)

                         cGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACg                                                                <  1:7881286/93‑1 (MQ=255)

                          gATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGGGGTTTCCCGGTTTACCCTCGCCGCAGATACCGGACgg                                                               >  1:2196336/1‑93 (MQ=255)

                            tGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACgggg                                                             >  2:821191/1‑93 (MQ=255)

                            tGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACgggg                                                             <  1:3294723/93‑1 (MQ=255)

                            tGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACgggg                                                             <  2:6785412/93‑1 (MQ=255)

                                 aaGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGCGGTTTACCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATgc                                                        >  1:3852280/1‑93 (MQ=255)

                                   ggCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATgcgc                                                      >  1:3757000/1‑93 (MQ=255)

                                     cAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGccc                                                    >  2:5235664/1‑93 (MQ=255)

                                      aaGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGCGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCa                                                   >  2:2882415/1‑93 (MQ=255)

                                       aGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCaa                                                  >  1:6163562/1‑93 (MQ=255)

                                        gATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAg                                                 <  2:1369736/93‑1 (MQ=255)

                                        gATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAg                                                 <  2:4978121/93‑1 (MQ=255)

                                        gATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAg                                                 <  2:4439865/93‑1 (MQ=255)

                                          tGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCCAGtt                                               >  2:454526/1‑93 (MQ=255)

                                          tGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGtt                                               >  1:2808796/1‑93 (MQ=255)

                                            aaTCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTtc                                             <  1:981046/93‑1 (MQ=255)

                                            aaTCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTtc                                             >  1:3474985/1‑93 (MQ=255)

                                              tctcAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTtctc                                           >  2:7265930/1‑93 (MQ=255)

                                                tcAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCa                                         >  2:7060015/1‑93 (MQ=255)

                                                 cAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGGTTCTCCaa                                        >  2:91105/1‑93 (MQ=255)

                                                  aaCGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCaaa                                       <  2:1505786/93‑1 (MQ=255)

                                                  aaCGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCaaa                                       <  1:6665278/93‑1 (MQ=255)

                                                  aaCGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCaaa                                       >  1:4572936/1‑93 (MQ=255)

                                                  aaCGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCaaa                                       <  2:6591277/93‑1 (MQ=255)

                                                    cGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAgg                                     <  2:5376689/93‑1 (MQ=255)

                                                       gCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGc                                  <  2:2553980/93‑1 (MQ=255)

                                                       gCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGc                                  <  1:473435/93‑1 (MQ=255)

                                                       gCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGc                                  <  2:2153138/93‑1 (MQ=255)

                                                         aaCAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCAt                                >  1:20056/1‑93 (MQ=255)

                                                         aaCAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCAt                                <  2:168447/93‑1 (MQ=255)

                                                          aCAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATc                               >  2:3611072/1‑93 (MQ=255)

                                                          aCAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATc                               >  2:3305283/1‑93 (MQ=255)

                                                          aCAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATc                               <  1:4036727/93‑1 (MQ=255)

                                                           cAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCa                              <  2:2088800/93‑1 (MQ=255)

                                                             ggcgggcgGCAGTGTTCTGTTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCAcc                            <  2:5138291/93‑1 (MQ=255)

                                                             ggcgggcgGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCAcg                            >  1:2284967/1‑92 (MQ=255)

                                                             ggcgggcgGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCAcc                            >  2:5971899/1‑93 (MQ=255)

                                                                gggcggCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGa                         <  2:3379776/93‑1 (MQ=255)

                                                                gggcggCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGa                         >  2:581221/1‑93 (MQ=255)

                                                                 ggcggcAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGAt                        <  2:5456551/93‑1 (MQ=255)

                                                                 ggcggcAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGAt                        >  1:7426148/1‑93 (MQ=255)

                                                                  gcggcAGTTTGCTGGTGGTTTCCCGGTTTGCCCTCGCCGCAGATACCGAACGGGGGATCGTCCCAAGTTTCTCCAAAGGTGCATCACCGGATc                       <  2:3143534/93‑1 (MQ=255)

                                                                  gcggcAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATc                       <  1:1863601/93‑1 (MQ=255)

                                                                   cggcAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATcg                      >  1:5425489/1‑93 (MQ=255)

                                                                    ggcAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATcgc                     <  1:4256455/93‑1 (MQ=255)

                                                                      cAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCgcgc                   >  2:6340237/1‑93 (MQ=255)

                                                                         tgtgCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCaga                <  1:2934346/93‑1 (MQ=255)

                                                                          gtgCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCagag               <  1:1770257/93‑1 (MQ=255)

                                                                           tgCTTGGGGTTTCCCGGTTTCCCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAgg              <  2:1463204/93‑1 (MQ=255)

                                                                            gCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGTATCGCGCAGAGGc             <  1:2036482/93‑1 (MQ=255)

                                                                            gCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGc             <  2:7859903/93‑1 (MQ=255)

                                                                            gCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGc             >  1:1474612/1‑93 (MQ=255)

                                                                             cTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCg            <  2:567562/93‑1 (MQ=255)

                                                                               ggtgatTTCCCGGTTTACCCTCGCCGCAGACACCGAGCGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGtt          <  2:4088751/93‑1 (MQ=255)

                                                                                gtggtTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTta         <  1:1297841/93‑1 (MQ=255)

                                                                                 tggtTTTCCTGTTTCCCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTtat        <  2:2643185/93‑1 (MQ=255)

                                                                                 tggtTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTtat        >  1:7593127/1‑93 (MQ=255)

                                                                                   gtTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTTatat      >  2:3546740/1‑93 (MQ=255)

                                                                                    tttCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTTATATg     <  2:1930972/93‑1 (MQ=255)

                                                                                     ttCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTTATATGa    <  1:1707053/93‑1 (MQ=255)

GTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTTATATGACT  >  NC_000913/4074737‑4074916

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