BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,075,220	(T)_5→6	coding (256/357 nt)	yggP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,075,215	1	.	T	94.9%	288.7 / 9.7	98	coding (261/357 nt)	yggP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (3/2); new base T (46/47); total (49/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 112 total aligned reads displayed.

CAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑C‑TT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGA  >  NC_000913/3075140‑3075285

cAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑c‑c                                                                           >  2:4783908/1‑76 (MQ=255)

cAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑c‑c                                                                           <  1:4301406/76‑1 (MQ=255)

 aaaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑C‑t                                                                         <  1:5540457/76‑1 (MQ=255)

 aaaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑C‑t                                                                         <  1:3603436/76‑1 (MQ=255)

  aaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑C‑tt                                                                        <  1:3619194/76‑1 (MQ=255)

  aaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑C‑t                                                                         >  1:3603858/1‑75 (MQ=255)

    atatTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑C‑tt‑t‑t                                                                    >  2:3225749/1‑76 (MQ=255)

    agatTCGGAGCTGTAGCCCGCTGAAGAGCCGGTTGGTAACCCCATCGCCGGTTGCTAGACAAAACGCTGG‑C‑C‑at‑t‑t                                                                    >  1:2206883/3‑76 (MQ=255)

     tngtCGTGGCTGGAGCCCGCTGAATATCCTCTTGGTAACCCCATCGCCGGTTGCAATCCAAAACGCTGG‑C‑C‑tt‑t‑tt                                                                   <  2:5686839/73‑1 (MQ=255)

     tatTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑C‑tt‑t‑tt                                                                   <  1:1084890/76‑1 (MQ=255)

       ttCGTGGCTGTAGCGCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATGCAAAACGCTGG‑C‑C‑tt‑t‑ttta                                                                 <  2:3636440/76‑3 (MQ=255)

         cGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATAGCCGGTTGCAATACAAAACGCTGG‑C‑C‑TTCT‑Ttata                                                                >  2:3739168/1‑76 (MQ=255)

            aGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATAttt                                                             >  2:930693/1‑76 (MQ=255)

            aGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATAttt                                                             <  2:5054422/76‑1 (MQ=255)

            aGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATAttt                                                             <  1:491271/76‑1 (MQ=255)

            aGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATAttt                                                             <  1:2216954/76‑1 (MQ=255)

             gCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTg                                                            <  1:1472173/76‑1 (MQ=255)

              cTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGc                                                           >  1:3939930/1‑76 (MQ=255)

               tGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGcc                                                          >  2:1531097/1‑76 (MQ=255)

               tGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGcc                                                          >  1:3518940/1‑76 (MQ=255)

                gTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATGCAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCa                                                         <  2:5375189/76‑1 (MQ=255)

                gTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCa                                                         >  2:5198832/1‑76 (MQ=255)

                gTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCa                                                         <  1:2238017/76‑1 (MQ=255)

                gTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCa                                                         >  1:1251043/1‑76 (MQ=255)

                gTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCa                                                         <  1:4972499/76‑1 (MQ=255)

                 tAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAg                                                        >  2:4941880/1‑76 (MQ=255)

                  aGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAGCGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGt                                                       <  2:2191455/76‑1 (MQ=255)

                  aGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGAC‑T‑TT‑T‑TTATATTTGCCAGt                                                       <  1:603691/76‑1 (MQ=255)

                    cccGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTga                                                     <  2:2017352/76‑1 (MQ=255)

                    cccGCTGAATATCCGCGTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTGGCCAGTga                                                     <  2:2018620/76‑1 (MQ=255)

                     ccGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTgag                                                    <  1:930693/76‑1 (MQ=255)

                      cGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTgaga                                                   >  1:5206305/1‑76 (MQ=255)

                        cTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGAtt                                                 <  1:4666554/76‑1 (MQ=255)

                         tGAATATCGGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTTTCTTTGCCCGTGAGAtct                                                >  2:5807476/1‑74 (MQ=255)

                         tGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGAttt                                                >  2:5086370/1‑76 (MQ=255)

                         tGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGAttt                                                >  2:4160073/1‑76 (MQ=255)

                           aaTATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGAttttt                                              >  2:1773089/1‑76 (MQ=255)

                            atatCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTa                                             >  1:4650732/1‑76 (MQ=255)

                            atatCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑C‑TTGT‑TTATATTTGCCAGTGAGATTTTTa                                             >  2:4377747/1‑76 (MQ=255)

                             tatCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTAc                                            >  1:2651636/1‑76 (MQ=255)

                             tatCCGCGTCGTAACCTCATCGCCGGTTGCAATACAAAACGCTGT‑C‑CTTT‑T‑TTATATTTGCCAGTGCGAttttttc                                            <  1:129556/76‑3 (MQ=255)

                              atCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑CAT‑TT‑T‑TTATATTTGCCAGTGAGATTTTTAcc                                           >  2:887872/1‑76 (MQ=255)

                               tCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTAccc                                          <  1:5739547/76‑1 (MQ=255)

                               tCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTAccc                                          <  1:183339/76‑1 (MQ=255)

                                ccGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCa                                         <  1:5305323/76‑1 (MQ=255)

                                ccGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCa                                         <  2:3791354/76‑1 (MQ=255)

                                 cGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCAc                                        <  1:5057703/76‑1 (MQ=255)

                                 cGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGACTTTTACCCAc                                        <  1:3195550/76‑1 (MQ=255)

                                  gCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACt                                       <  2:5465476/76‑1 (MQ=255)

                                   cTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACtt                                      <  2:1950977/76‑1 (MQ=255)

                                    ttGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTc                                     >  2:3940497/1‑76 (MQ=255)

                                       gTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGAc                                  <  2:3007938/76‑1 (MQ=255)

                                        tAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACa                                 <  1:971857/76‑1 (MQ=255)

                                        tAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACa                                 >  2:3576394/1‑76 (MQ=255)

                                         aaCCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACaa                                <  2:2554484/76‑1 (MQ=255)

                                         aaCCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACaa                                >  1:4834619/1‑76 (MQ=255)

                                          aCCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAAt                               <  2:3939929/76‑1 (MQ=255)

                                           ccccATCGCCGGGTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATa                              <  2:640793/76‑1 (MQ=255)

                                           ccccAGCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTCCCCACTTCGACAATa                              <  2:1883535/76‑1 (MQ=255)

                                             ccATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAc                            >  1:930363/1‑76 (MQ=255)

                                              cATCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcc                           >  2:2064767/1‑76 (MQ=255)

                                               aTCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAccc                          >  1:5231181/1‑76 (MQ=255)

                                                tCGCCGGTTGCAATCCACAACGCTGG‑C‑C‑TT‑TGTTATATTTGCCCGTGATATTTTTACCCACTTCGACAGTAGcccc                         <  1:945322/76‑1 (MQ=255)

                                                tCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcccc                         <  1:909421/76‑1 (MQ=255)

                                                tCGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcccc                         <  2:1478896/76‑1 (MQ=255)

                                                gcGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcccc                         <  2:2083791/75‑1 (MQ=255)

                                                 cGCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCg                        >  2:4868022/1‑76 (MQ=255)

                                                  gCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCgg                       >  1:61909/1‑76 (MQ=255)

                                                  gCCGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCgg                       >  2:3885827/1‑76 (MQ=255)

                                                    cGGTTGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGACATTTTTACCCACTTCGACAATAACCCCGGca                     >  2:1334287/1‑76 (MQ=255)

                                                      gTTGCGATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGcaca                   <  2:2212405/76‑1 (MQ=255)

                                                       ttGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTGTATTTGCCAGTGAGACTTTTACCCACTTCGACAATAACCCCGGCACAt                  <  2:3128164/76‑1 (MQ=255)

                                                       ttGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATa                              >  1:4092099/1‑64 (MQ=255)

                                                       ttGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATa                              <  2:4092098/64‑1 (MQ=255)

                                                       ttGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACAt                  <  2:3207923/76‑1 (MQ=255)

                                                       ttGCAATACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACAt                  <  2:1170368/76‑1 (MQ=255)

                                                           aaTACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCAt              <  2:3533883/76‑1 (MQ=255)

                                                            aTACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcc                           >  2:3137811/1‑62 (MQ=255)

                                                            aTACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcc                           <  1:3137811/62‑1 (MQ=255)

                                                            aTACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATg             <  2:812041/76‑1 (MQ=255)

                                                             tACAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGc            >  1:2649840/1‑76 (MQ=255)

                                                               cAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGccc          <  1:1066654/76‑1 (MQ=255)

                                                               cAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGccc          >  1:4519077/1‑76 (MQ=255)

                                                               cAAAACGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGccc          >  1:5671008/1‑76 (MQ=255)

                                                                aaaaCGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCg         >  2:1802260/1‑76 (MQ=255)

                                                                 aaaCGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCgg        <  1:1830024/76‑1 (MQ=255)

                                                                  aaCGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATACCCCCGGCACATTCATGCCCGGt       <  1:4377748/76‑1 (MQ=255)

                                                                  aaCGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGt       <  2:4420464/76‑1 (MQ=255)

                                                                  aaCGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGt       <  2:3518940/76‑1 (MQ=255)

                                                                  aaCGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGt       >  2:1723819/1‑76 (MQ=255)

                                                                   aCGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACCTTCATGCCCGGTa      >  2:1629259/1‑76 (MQ=255)

                                                                   aCGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTa      >  1:4603441/1‑76 (MQ=255)

                                                                   aCGCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTa      >  2:2704022/1‑76 (MQ=255)

                                                                     gCTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAAt    >  1:5586623/1‑76 (MQ=255)

                                                                      cTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATg   >  1:1539790/1‑76 (MQ=255)

                                                                      cTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATg   >  1:5150359/1‑76 (MQ=255)

                                                                      cTGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATg   >  1:2431312/1‑76 (MQ=255)

                                                                       tGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGa  >  2:4139416/1‑76 (MQ=255)

                                                                       tGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGa  >  1:1878416/1‑76 (MQ=255)

                                                                       tGG‑C‑CTTT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGa  <  1:5198833/76‑1 (MQ=255)

CAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGG‑C‑C‑TT‑T‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGA  >  NC_000913/3075140‑3075285

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