BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,298,395	A→G	L72P (CTG→CCG)	yjcP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,298,395	0	A	G	81.0%	627.2 / 161.3	315	L72P (CTG→CCG)	yjcP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (38/22); new base G (143/112); total (181/134)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 316 total aligned reads displayed.

CGGCATCTAACAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGAT  >  NC_000913/4298303‑4298486

  gCATCTAACAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGgc                                                                                           >  2:4245043/1‑93 (MQ=255)

  gCATCTAACAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGgc                                                                                           >  1:3084366/1‑93 (MQ=255)

  gCATCTAACAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGgc                                                                                           >  2:2408911/1‑93 (MQ=255)

   cATCTAACAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGgcg                                                                                          <  2:4121404/93‑1 (MQ=255)

    aTCTAACAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGgcgc                                                                                         >  1:5507227/1‑93 (MQ=255)

    aTCTAACAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGgcgc                                                                                         >  1:6076136/1‑93 (MQ=255)

     tCTAACAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCa                                                                                        >  2:4797117/1‑93 (MQ=255)

      cTAACAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCAt                                                                                       >  1:2295013/1‑93 (MQ=255)

       tAACAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATc                                                                                      >  2:1161773/1‑93 (MQ=255)

          cAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAg                                                                                   >  1:21059/1‑93 (MQ=255)

              tCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGc                                                                               <  2:3967060/93‑1 (MQ=255)

              tCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGc                                                                               >  1:2750104/1‑93 (MQ=255)

              tCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGGGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGc                                                                               >  2:5359023/1‑93 (MQ=255)

               cGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCg                                                                              <  1:1182711/93‑1 (MQ=255)

                 gCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCggg                                                                            >  2:4999454/1‑93 (MQ=255)

                   cTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGGATCCAGCTGCGGGTc                                                                          >  1:3639895/1‑93 (MQ=255)

                   cTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTc                                                                          >  2:330201/1‑93 (MQ=255)

                   cTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTc                                                                          >  1:1342181/1‑93 (MQ=255)

                    tGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCa                                                                         >  2:5085871/1‑93 (MQ=255)

                    tGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCa                                                                         >  2:4185456/1‑93 (MQ=255)

                     gggACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCATTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCAt                                                                        <  2:7385865/93‑1 (MQ=255)

                     gggACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCAt                                                                        >  2:7459648/1‑93 (MQ=255)

                     gggACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCAt                                                                        >  1:1931900/1‑93 (MQ=255)

                     gggACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCAt                                                                        <  2:5772146/93‑1 (MQ=255)

                      ggACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCAtt                                                                       <  1:2795923/93‑1 (MQ=255)

                       gACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTg                                                                      >  1:4115288/1‑93 (MQ=255)

                         cTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAg                                                                    <  1:3439784/93‑1 (MQ=255)

                          tGCGCTTTTTCTTCGCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAgt                                                                   >  2:5190992/1‑93 (MQ=255)

                          tGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAgt                                                                   <  2:2888809/93‑1 (MQ=255)

                           gcgcTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAgtt                                                                  <  1:4038734/93‑1 (MQ=255)

                              cTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAgttgtt                                                               <  1:365924/93‑1 (MQ=255)

                                ttttCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTc                                                             >  1:521805/1‑93 (MQ=255)

                                 tttCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTcc                                                            >  2:1388549/1‑93 (MQ=255)

                                  ttcttcCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTcca                                                           <  2:6033741/93‑1 (MQ=255)

                                   tcttcCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTccac                                                          >  1:2579040/1‑93 (MQ=255)

                                    cttcCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTccacc                                                         <  2:6317434/93‑1 (MQ=255)

                                      tcCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCcaccac                                                       <  1:664962/93‑1 (MQ=255)

                                      tcCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCcaccac                                                       <  1:2476700/93‑1 (MQ=255)

                                         cGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTgc                                                    <  2:3806910/93‑1 (MQ=255)

                                         cGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTgc                                                    <  2:4490925/93‑1 (MQ=255)

                                         cGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTgc                                                    <  2:1926766/93‑1 (MQ=255)

                                          gCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTgcg                                                   >  1:3912665/1‑93 (MQ=255)

                                           cAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTgcgc                                                  <  2:3587571/93‑1 (MQ=255)

                                            aGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGcc                                                 >  2:4576461/1‑93 (MQ=255)

                                             gTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCt                                                >  1:4851350/1‑93 (MQ=255)

                                             gTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCt                                                <  2:2157306/93‑1 (MQ=255)

                                               ttCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGc                                              >  1:3522063/1‑93 (MQ=255)

                                                tCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCg                                             >  1:1391731/1‑93 (MQ=255)

                                                  gCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGc                                           >  1:3133682/1‑93 (MQ=255)

                                                  gCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGc                                           >  2:72041/1‑93 (MQ=255)

                                                   cTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGcc                                          >  1:6738055/1‑93 (MQ=255)

                                                   cTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGcc                                          <  1:5684975/93‑1 (MQ=255)

                                                    ttCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCa                                         >  1:4582147/1‑93 (MQ=255)

                                                     tCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAg                                        <  2:5498847/93‑1 (MQ=255)

                                                      cGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGc                                       <  2:816298/93‑1 (MQ=255)

                                                      cGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGc                                       >  1:7038010/1‑93 (MQ=255)

                                                       ggCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGcc                                      >  1:237282/1‑93 (MQ=255)

                                                        gCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCg                                     >  1:767648/1‑93 (MQ=255)

                                                          gAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCCCCACTGCGCCTGCGGCCAGCCGGa                                   <  1:6268843/93‑1 (MQ=255)

                                                          gAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGa                                   >  1:3219595/1‑93 (MQ=255)

                                                          gAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGa                                   >  2:3724595/1‑93 (MQ=255)

                                                           aGGGTGTGTGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAg                                  <  2:4834938/93‑1 (MQ=255)

                                                              gtgtgtGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTg                               <  2:6965954/93‑1 (MQ=255)

                                                              gtgtgtGAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTg                               <  2:6389848/93‑1 (MQ=255)

                                                               tgtgtgAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGGGCCgg                              >  1:7466093/1‑93 (MQ=255)

                                                               tgtgtgAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTgg                              >  2:6547736/1‑93 (MQ=255)

                                                               tgtgtgAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTgg                              <  2:128922/93‑1 (MQ=255)

                                                               tgtgtgAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTgg                              >  2:1482183/1‑93 (MQ=255)

                                                                 tgtgAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCa                            >  2:3402157/1‑93 (MQ=255)

                                                                    gAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGa                         >  1:6410961/1‑93 (MQ=255)

                                                                    gAACCACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGa                         >  1:2189686/1‑93 (MQ=255)

                                                                     aaCCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGAt                        >  2:6805861/1‑93 (MQ=255)

                                                                        cACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGaa                     >  1:5343748/1‑93 (MQ=255)

                                                                        cACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGaa                     <  2:5908746/93‑1 (MQ=255)

                                                                        cACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGaa                     <  1:6889671/93‑1 (MQ=255)

                                                                        cACTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGaa                     >  2:5604780/1‑93 (MQ=255)

                                                                        cACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGaa                     <  1:1053249/93‑1 (MQ=255)

                                                                         aCTTAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAat                    <  1:4636082/93‑1 (MQ=255)

                                                                           ttAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAatat                  <  2:5551651/93‑1 (MQ=255)

                                                                           ttAGCGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAatat                  <  1:4178056/93‑1 (MQ=255)

                                                                              gCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATAtcgt               <  2:2267224/93‑1 (MQ=255)

                                                                              gCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATAtcgt               <  2:3085897/93‑1 (MQ=255)

                                                                              gCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATAtcgt               >  1:61505/1‑93 (MQ=255)

                                                                               cGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATAtcgtc              <  1:3758212/93‑1 (MQ=255)

                                                                               cGTCCGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATAtcgtc              <  1:5004919/93‑1 (MQ=255)

                                                                                   cGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGcc          >  2:2622218/1‑93 (MQ=255)

                                                                                   cGTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGcc          <  1:1588434/93‑1 (MQ=255)

                                                                                    gTTGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCa         <  1:5278825/93‑1 (MQ=255)

                                                                                     ttGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAg        <  1:1026534/93‑1 (MQ=255)

                                                                                     ttGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAg        >  1:654913/1‑93 (MQ=255)

                                                                                     ttGGATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAg        <  1:7302948/93‑1 (MQ=255)

                                                                                       ggATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGtt      >  1:2923498/1‑93 (MQ=255)

                                                                                       ggATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGtt      <  1:1564160/93‑1 (MQ=255)

                                                                                       ggATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGtt      <  2:1721502/93‑1 (MQ=255)

                                                                                       ggATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGtt      <  1:5863247/93‑1 (MQ=255)

                                                                                        gATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGttt     >  2:6683618/1‑93 (MQ=255)

                                                                                        gATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGttt     >  2:575474/1‑93 (MQ=255)

                                                                                        gATCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGttt     <  2:774443/93‑1 (MQ=255)

                                                                                         aTCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTg    <  1:7505439/93‑1 (MQ=255)

                                                                                         aTCGGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGCTTg    >  2:6216544/1‑93 (MQ=255)

CGGCATCTAACAAATCGGCCTGGGACTGCGCTTTTTCTTCCCGCAGTTTCGCTTCGGCGAGGGTGTGTGAACCACTTAGCGTCCGTTGGATCAGCGCATCCAGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGAT  >  NC_000913/4298303‑4298486

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