BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,305,969	Δ1 bp	coding (271/1461 nt)	cls ←	cardiolipin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,305,969	0	T	.	98.7%	761.5 / 2.3	154	coding (271/1461 nt)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand): ref base T (2/0); new base . (92/60); total (94/60)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 155 total aligned reads displayed.

AGCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCAGC  >  REL606/1305898‑1306038

aGCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCt                                                                  >  2:1700534/1‑76 (MQ=255)

 gCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCtt                                                                 <  2:2981498/76‑1 (MQ=255)

 gCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCtt                                                                 <  2:5140174/76‑1 (MQ=255)

  cGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTa                                                                >  1:5816392/1‑76 (MQ=255)

   gATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTAc                                                               >  2:3263248/1‑76 (MQ=255)

   gATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTAc                                                               <  1:4825954/76‑1 (MQ=255)

   gATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTAc                                                               <  2:2852609/76‑1 (MQ=255)

    aTCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACa                                                              >  2:60555/1‑76 (MQ=255)

     tCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAg                                                             >  2:3868457/1‑76 (MQ=255)

       cccTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGc                                                           <  1:386057/76‑1 (MQ=255)

       cccTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGc                                                           >  1:3802775/1‑76 (MQ=255)

        ccTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCt                                                          <  1:401246/76‑1 (MQ=255)

        ccTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCt                                                          >  1:3259380/1‑76 (MQ=255)

         cTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCtt                                                         <  1:690289/76‑1 (MQ=255)

             cgacgCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAg                                                     <  2:2894571/76‑1 (MQ=255)

              gacgCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAgg                                                    <  2:4078854/76‑1 (MQ=255)

              gacgCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAgg                                                    <  1:6306749/76‑1 (MQ=255)

              gacgCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGATATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAgg                                                    <  1:1418929/76‑1 (MQ=255)

                 gCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCg                                                 >  1:611474/1‑76 (MQ=255)

                 gCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCg                                                 >  1:1204695/1‑76 (MQ=255)

                    cGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTa                                              >  2:4133096/1‑76 (MQ=255)

                    cGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTa                                              >  2:2181596/1‑76 (MQ=255)

                    cGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTa                                              <  1:2079635/76‑1 (MQ=255)

                    cGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTa                                              <  1:167813/76‑1 (MQ=255)

                     gCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTaa                                             <  2:6157756/76‑1 (MQ=255)

                      cAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTaa                                             >  2:4488183/1‑76 (MQ=255)

                        aaGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGcc                                          >  1:1184587/1‑76 (MQ=255)

                         aGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCa                                         >  2:3312516/1‑76 (MQ=255)

                         aGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCa                                         >  2:4852848/1‑76 (MQ=255)

                         aGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCa                                         >  2:2732771/1‑76 (MQ=255)

                          gCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCAt                                        <  1:4703262/76‑1 (MQ=255)

                          gCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCAt                                        <  1:4904222/76‑1 (MQ=255)

                          gCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCAt                                        >  2:5394928/1‑76 (MQ=255)

                          gCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCAt                                        <  2:4698908/76‑1 (MQ=255)

                           cTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCAtt                                       <  2:4887898/76‑1 (MQ=255)

                           cTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCAtt                                       >  1:2704850/1‑76 (MQ=255)

                           cTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCAtt                                       >  1:6609069/1‑76 (MQ=255)

                            ttGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCAttt                                      >  2:4841539/1‑76 (MQ=255)

                            ttGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCAttt                                      >  2:5141119/1‑76 (MQ=255)

                            ttGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCAttt                                      >  1:5040865/1‑76 (MQ=255)

                             tGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCAtttt                                     >  1:6407414/1‑76 (MQ=255)

                              gAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATTCTTACAGGCTTTAAGGTCGTTAAGCCATTTTg                                    <  1:5956495/76‑1 (MQ=255)

                              gAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTg                                    <  1:5762723/76‑1 (MQ=255)

                               aataatGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGc                                   >  2:3145172/1‑76 (MQ=255)

                               aataatGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGc                                   >  2:1463198/1‑76 (MQ=255)

                                ataatGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCg                                  >  2:3399016/1‑76 (MQ=255)

                                 taatGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCgg                                 >  2:3724002/1‑76 (MQ=255)

                                  aatGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGt                                >  2:6529845/1‑76 (MQ=255)

                                   atGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTAAAGCCATTTTGCGGTg                               >  1:4148931/1‑76 (MQ=255)

                                    tGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTgg                              >  1:6262477/1‑76 (MQ=255)

                                    tGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTgg                              >  2:4222127/1‑76 (MQ=255)

                                    tGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTgg                              >  2:2048560/1‑76 (MQ=255)

                                    tGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTgg                              >  1:1053212/1‑76 (MQ=255)

                                     ggCGCAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGa                             <  2:1910030/76‑1 (MQ=255)

                                      gcgcAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGaa                            <  1:260640/76‑1 (MQ=255)

                                       cgcAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAg                           >  2:763343/1‑76 (MQ=255)

                                         cAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGc                         >  2:4697842/1‑76 (MQ=255)

                                         cAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGc                         >  1:1995197/1‑76 (MQ=255)

                                         cAGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGc                         >  1:4519066/1‑76 (MQ=255)

                                          aGCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGcc                        >  1:3517234/1‑76 (MQ=255)

                                          aGCGACACTGCTATTTTCATCGGCGAAG‑GATGCTTACAGGCTTTAAGGGCGTTAAGCCATTTTGCGGTGGAAGGcc                        <  1:2679549/76‑1 (MQ=255)

                                           gCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCca                       <  1:2990024/76‑1 (MQ=255)

                                           gCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCca                       <  2:2773504/76‑1 (MQ=255)

                                           gCGACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCca                       <  2:3868009/76‑1 (MQ=255)

                                             gACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCcaca                     >  2:5356337/1‑76 (MQ=255)

                                             gACACTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCcaca                     <  2:6288228/76‑1 (MQ=255)

                                              acacTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACAt                    <  1:1056313/76‑1 (MQ=255)

                                              acacTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACAt                    <  2:1865826/76‑1 (MQ=255)

                                               cacTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATc                   >  1:6413720/1‑76 (MQ=255)

                                                acTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCg                  <  2:2599649/76‑1 (MQ=255)

                                                acTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCg                  <  1:1405807/76‑1 (MQ=255)

                                                acTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCg                  >  1:4403194/1‑76 (MQ=255)

                                                acTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCg                  <  1:1690799/76‑1 (MQ=255)

                                                 cTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGc                 <  2:2429192/76‑1 (MQ=255)

                                                 cTGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGc                 <  2:2004331/76‑1 (MQ=255)

                                                  tGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGct                >  1:1496456/1‑76 (MQ=255)

                                                  tGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGct                <  2:5496005/76‑1 (MQ=255)

                                                  tGCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGct                <  2:5966149/76‑1 (MQ=255)

                                                   gCTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGctc               >  1:2982572/1‑76 (MQ=255)

                                                    cTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGctct              >  1:805601/1‑76 (MQ=255)

                                                    cTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGctct              >  1:5060134/1‑76 (MQ=255)

                                                    cTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGctct              >  2:6153193/1‑76 (MQ=255)

                                                    cTATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGctct              <  1:2404811/76‑1 (MQ=255)

                                                     tATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTg             >  2:4504421/1‑76 (MQ=255)

                                                     tATTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTg             <  2:4531277/76‑1 (MQ=255)

                                                      aTTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTgg            <  1:1803108/76‑1 (MQ=255)

                                                      aTTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTgg            >  1:5327001/1‑76 (MQ=255)

                                                      aTTTTCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTgg            >  2:6416197/1‑76 (MQ=255)

                                                       ttttCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGgc           >  1:714457/1‑76 (MQ=255)

                                                       ttttCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGgc           >  2:718925/1‑76 (MQ=255)

                                                        tttCTTCGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGgcg          >  2:2363651/1‑76 (MQ=255)

                                                          tcttcGGCGAAGA‑ATGCTTACAGGCTTTAAGGCCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGgcgcg        >  2:2937725/1‑76 (MQ=255)

                                                           cttcGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGgcgcgc       >  1:3594993/1‑76 (MQ=255)

                                                            ttcGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCt      >  1:6195124/1‑76 (MQ=255)

                                                             tcGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTc     <  1:5251487/76‑1 (MQ=255)

                                                              cGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCa    <  2:4921633/76‑1 (MQ=255)

                                                              cGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCa    <  2:2179859/76‑1 (MQ=255)

                                                              cGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCa    >  2:4334226/1‑76 (MQ=255)

                                                              cGGCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCa    >  1:146118/1‑76 (MQ=255)

                                                               ggCGAAGA‑ATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCAg   >  2:1292625/1‑76 (MQ=255)

AGCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCAGC  >  REL606/1305898‑1306038

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