BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	177,386	C→T	R80H (CGC→CAC)	yadS ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	177,386	0	C	T	99.6%	954.8 / NA	265	R80H (CGC→CAC)	yadS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (140/124); minor base G (1/0); total (141/124)
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.99e-01

Only 100 of 290 total aligned reads displayed.

AAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGC  >  NC_000913/177296‑177475

aaGGCTTTATTCACACCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGgtgc                                                                                         >  1:2435079/1‑93 (MQ=255)

 aGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCa                                                                                        >  1:297657/1‑93 (MQ=255)

  ggCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCAc                                                                                       >  2:4156433/1‑93 (MQ=255)

  ggCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCAc                                                                                       <  1:3702598/93‑1 (MQ=255)

   gCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCAcc                                                                                      <  1:1338640/93‑1 (MQ=255)

     tttATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAg                                                                                    >  2:6174391/1‑93 (MQ=255)

     tttATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAg                                                                                    >  1:189083/1‑93 (MQ=255)

     tttATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAg                                                                                    >  1:1967924/1‑93 (MQ=255)

       tATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCa                                                                                  >  2:2851952/1‑93 (MQ=255)

         ttCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTTGTGCACCAGCACg                                                                                >  2:3802055/1‑93 (MQ=255)

         ttCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACg                                                                                <  2:4118/93‑1 (MQ=255)

         ttCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACg                                                                                >  1:2919431/1‑93 (MQ=255)

           cACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGAt                                                                              <  1:2180189/93‑1 (MQ=255)

           cACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGAt                                                                              <  1:4377373/93‑1 (MQ=255)

           cACGCCAATACAGACAACCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGAt                                                                              <  1:3624490/93‑1 (MQ=255)

             cGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATgg                                                                            <  1:5635535/93‑1 (MQ=255)

              gCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGt                                                                           >  1:2318451/1‑93 (MQ=255)

              gCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGt                                                                           >  1:1572000/1‑93 (MQ=255)

               ccAATGCCGACAAACACCGCCAGACCAACGGGGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTc                                                                          >  2:4451741/1‑93 (MQ=255)

               ccAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTc                                                                          <  2:552380/93‑1 (MQ=255)

                cAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCa                                                                         >  1:5737065/1‑93 (MQ=255)

                cAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCa                                                                         <  2:1217355/93‑1 (MQ=255)

                cAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCa                                                                         >  2:2082871/1‑93 (MQ=255)

                 aaTGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTccg                                                                        >  1:6553139/1‑91 (MQ=255)

                  aTGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGc                                                                       >  2:5056247/1‑93 (MQ=255)

                     ccGACAAACACCACCAGACCAACGTCGTCCAGCACCGGCATCATCCATTTTGGTAAGCGTCTTGACTGGTGCACCAGCACGATGGTCAGCATg                                                                    <  1:4124398/93‑1 (MQ=255)

                      cGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGc                                                                   >  2:5222657/1‑93 (MQ=255)

                      cGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGc                                                                   >  2:5682413/1‑93 (MQ=255)

                       gACAAACACCGCCGGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCt                                                                  >  1:6438589/1‑93 (MQ=255)

                        aCAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTg                                                                 >  1:4891597/1‑93 (MQ=255)

                         cAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTgg                                                                <  2:5475632/93‑1 (MQ=255)

                          aaaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGt                                                               <  1:5713691/93‑1 (MQ=255)

                           aaCGCCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTg                                                              <  2:48235/93‑1 (MQ=255)

                           aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTg                                                              >  2:2550017/1‑93 (MQ=255)

                           aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTg                                                              >  1:5059751/1‑93 (MQ=255)

                           aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTg                                                              >  2:4374602/1‑93 (MQ=255)

                            acacCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGa                                                             <  1:1022607/93‑1 (MQ=255)

                             cacCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGAc                                                            >  2:625079/1‑93 (MQ=255)

                              cccGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGAcc                                                           <  2:1957050/92‑1 (MQ=255)

                              acCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGAcc                                                           >  2:6337393/1‑93 (MQ=255)

                              acCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGAcc                                                           <  2:1915618/93‑1 (MQ=255)

                                cGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTACGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCAt                                                         <  2:2608389/93‑1 (MQ=255)

                                 gCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCAtt                                                        <  2:2991483/93‑1 (MQ=255)

                                 gCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCAtt                                                        >  1:5291065/1‑93 (MQ=255)

                                 gCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCAtt                                                        >  2:5915049/1‑93 (MQ=255)

                                 gCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCAtt                                                        >  1:148373/1‑93 (MQ=255)

                                  ccAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTg                                                       >  2:2968682/1‑93 (MQ=255)

                                    aGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCa                                                     >  2:2189595/1‑93 (MQ=255)

                                    aGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCa                                                     >  1:3832657/1‑93 (MQ=255)

                                     gACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCaa                                                    >  2:4489413/1‑93 (MQ=255)

                                      aCCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAAc                                                   <  1:2637466/93‑1 (MQ=255)

                                      aCCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAAc                                                   <  1:4040595/93‑1 (MQ=255)

                                         aaCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGAc                                                >  1:4593000/1‑93 (MQ=255)

                                         aaCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGAc                                                >  2:6564148/1‑93 (MQ=255)

                                         aaCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTAGTGACCATTGCAACGAc                                                >  2:6676480/1‑93 (MQ=255)

                                           cGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCa                                              <  1:1112794/93‑1 (MQ=255)

                                             gCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGa                                            <  1:4683105/93‑1 (MQ=255)

                                              cGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGAt                                           >  1:3933374/1‑93 (MQ=255)

                                              cGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGAt                                           <  1:2123732/93‑1 (MQ=255)

                                                tCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCg                                         >  2:5728325/1‑93 (MQ=255)

                                                  cAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGt                                       <  1:4171556/93‑1 (MQ=255)

                                                  cAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGt                                       <  1:1456283/93‑1 (MQ=255)

                                                   aGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTg                                      >  1:1845896/1‑93 (MQ=255)

                                                   aGCACCGGCAACATCCATGTTGGTAAGCGTCTTGGCTGGTGCACTAGCACGATGGTCAGCATGCTGGTGACCATCGCAACGACCAGATCGGTg                                      >  1:6734190/1‑93 (MQ=255)

                                                       ccGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGAt                                  >  2:1335459/1‑93 (MQ=255)

                                                          gCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCtt                               <  2:2519581/93‑1 (MQ=255)

                                                          gCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCtt                               <  1:3433714/93‑1 (MQ=255)

                                                          gCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCtt                               <  2:2757522/93‑1 (MQ=255)

                                                          gCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCtt                               <  2:1808685/93‑1 (MQ=255)

                                                          gCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCtt                               <  2:323706/93‑1 (MQ=255)

                                                           cAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCttt                              <  2:3674417/93‑1 (MQ=255)

                                                             aCATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCa                            >  1:949397/1‑93 (MQ=255)

                                                             aCATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGCGACCATTGCAACGACCAGATCGGTGGGATCTTTCa                            >  1:5508010/1‑93 (MQ=255)

                                                                tacATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCAccc                         <  1:5573915/91‑1 (MQ=255)

                                                                tCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACNAGATCGGTGGGATCTNTCCncc                         >  2:5314181/1‑93 (MQ=255)

                                                                tCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCAccc                         >  1:6277839/1‑93 (MQ=255)

                                                                tCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCAccc                         >  2:5008855/1‑93 (MQ=255)

                                                                tCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATATTTCAacc                         >  2:3322785/1‑93 (MQ=255)

                                                                 ccATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCa                        >  1:675011/1‑93 (MQ=255)

                                                                    ttttGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaaaa                     >  1:6039781/1‑93 (MQ=255)

                                                                    ttttGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaaaa                     <  1:3066505/93‑1 (MQ=255)

                                                                    ttttGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaaaa                     >  2:789630/1‑93 (MQ=255)

                                                                      ttGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATa                   <  2:2260903/93‑1 (MQ=255)

                                                                       tGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATAc                  <  2:2324487/93‑1 (MQ=255)

                                                                        ggTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATAcc                 <  1:788743/93‑1 (MQ=255)

                                                                        ggTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATAcc                 <  1:2743187/93‑1 (MQ=255)

                                                                         gTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCg                <  1:6689924/93‑1 (MQ=255)

                                                                          tAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCgg               <  1:3268651/93‑1 (MQ=255)

                                                                             gCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGcc            >  1:5868260/1‑93 (MQ=255)

                                                                               gTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGt          <  2:4610771/93‑1 (MQ=255)

                                                                                tCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTg         >  1:3509993/1‑93 (MQ=255)

                                                                                tCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTg         >  2:2635417/1‑93 (MQ=255)

                                                                                tCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTg         >  2:4070278/1‑93 (MQ=255)

                                                                                  ttGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGAt       <  1:4505320/93‑1 (MQ=255)

                                                                                  ttGGCGGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGAt       <  2:1580991/93‑1 (MQ=255)

                                                                                    tgctgggACACCAGCACGTTGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATcc     <  1:4543024/92‑1 (MQ=255)

                                                                                    ggctggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATcc     >  1:4790013/1‑93 (MQ=255)

                                                                                    ggctggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATcc     <  1:4974941/93‑1 (MQ=255)

                                                                                    ggctggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATcc     <  1:180319/93‑1 (MQ=255)

                                                                                      ctggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAg   <  2:3601317/93‑1 (MQ=255)

AAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGC  >  NC_000913/177296‑177475

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