BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,100,603	A→G	F90L (TTC→CTC)	yggH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,100,603	0	A	G	99.6%	800.0 / ‑4.8	249	F90L (TTC→CTC)	yggH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (143/104); total (144/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 263 total aligned reads displayed.

GACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAA  >  NC_000913/3100513‑3100692

gACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGt                                                                                         <  1:5926448/93‑1 (MQ=255)

gACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGt                                                                                         <  2:323041/93‑1 (MQ=255)

  cacaTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTcc                                                                                       <  2:1356105/93‑1 (MQ=255)

  cacaTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTcc                                                                                       >  1:1731380/1‑93 (MQ=255)

  cacaTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTcc                                                                                       <  2:2497723/93‑1 (MQ=255)

  cacaTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTcc                                                                                       <  1:4042665/93‑1 (MQ=255)

   acaTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCt                                                                                      <  2:195697/93‑1 (MQ=255)

    caTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTg                                                                                     <  1:4750410/93‑1 (MQ=255)

     aTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGc                                                                                    <  1:5903417/93‑1 (MQ=255)

     aTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGc                                                                                    >  2:5082891/1‑93 (MQ=255)

     aTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGc                                                                                    >  1:3898635/1‑93 (MQ=255)

     aTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGCGAATGCACTTCAATGCCGAGGAGGTCCTGc                                                                                    >  2:2833172/1‑93 (MQ=255)

     aTCACGCGCAGGTTGCTTAAACCGTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGc                                                                                    >  1:5374208/1‑93 (MQ=255)

       cACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTc                                                                                  <  1:1989620/93‑1 (MQ=255)

        aCGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCa                                                                                 >  1:6432722/1‑93 (MQ=255)

        aCGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCa                                                                                 <  2:6439738/93‑1 (MQ=255)

        aCGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCa                                                                                 <  1:2145842/93‑1 (MQ=255)

         cgcgCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAg                                                                                <  2:1126672/93‑1 (MQ=255)

           cgcAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAggg                                                                              <  1:2519268/93‑1 (MQ=255)

             cAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGATGTCCTGCTCAGGGCg                                                                            <  1:3421927/93‑1 (MQ=255)

              aGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGa                                                                           >  2:4333785/1‑93 (MQ=255)

              aGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGa                                                                           >  1:142844/1‑93 (MQ=255)

              aGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGa                                                                           >  2:3351708/1‑93 (MQ=255)

               ggTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGAt                                                                          >  2:1783881/1‑93 (MQ=255)

               ggTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGAt                                                                          <  1:4621028/93‑1 (MQ=255)

               ggTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGAt                                                                          >  2:5344407/1‑93 (MQ=255)

                 ttGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGTGCGATCt                                                                        >  1:4760798/1‑93 (MQ=255)

                 ttGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCt                                                                        >  2:1286803/1‑93 (MQ=255)

                 ttGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCt                                                                        <  2:5767384/93‑1 (MQ=255)

                   gCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCttt                                                                      <  2:4260626/93‑1 (MQ=255)

                    cTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTa                                                                     >  2:5826981/1‑93 (MQ=255)

                    cTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTa                                                                     <  2:5228181/93‑1 (MQ=255)

                    cTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTa                                                                     <  2:3419189/93‑1 (MQ=255)

                    cTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTa                                                                     <  2:6380868/93‑1 (MQ=255)

                     ttAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCTATCTTTAg                                                                    >  1:3525161/1‑93 (MQ=255)

                     ttAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAg                                                                    >  2:1462321/1‑93 (MQ=255)

                     ttAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAg                                                                    >  2:2766789/1‑93 (MQ=255)

                        aaCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCa                                                                 <  2:4816844/93‑1 (MQ=255)

                         aCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCAt                                                                >  2:4950885/1‑93 (MQ=255)

                         aCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGACGAGGAGGTCCTACTCAGGGCGATCTTTAGCCAt                                                                >  1:526347/1‑93 (MQ=255)

                           cttcttCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTg                                                              >  2:840400/1‑93 (MQ=255)

                           cttcttCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTg                                                              <  1:444977/93‑1 (MQ=255)

                            ttcttcATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGc                                                             >  2:5990303/1‑93 (MQ=255)

                             tcttcATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGcc                                                            >  1:990629/1‑93 (MQ=255)

                               ttcATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGccac                                                          >  1:1001940/1‑93 (MQ=255)

                               ttcATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGccac                                                          >  1:309820/1‑93 (MQ=255)

                               ttcATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGccac                                                          >  1:1622590/1‑93 (MQ=255)

                                  aTGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAg                                                       >  2:2262278/1‑93 (MQ=255)

                                   tGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGc                                                      >  1:942834/1‑93 (MQ=255)

                                   tGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGc                                                      >  1:1425412/1‑93 (MQ=255)

                                     cgcAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGa                                                    >  1:3699122/1‑93 (MQ=255)

                                       cAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACg                                                  <  2:5501418/93‑1 (MQ=255)

                                        aGAAGCCAGGCACGCACCAACGCTCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGc                                                 >  2:4712671/1‑93 (MQ=255)

                                        aGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGc                                                 >  1:1033812/1‑93 (MQ=255)

                                        aGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGc                                                 >  1:4636363/1‑93 (MQ=255)

                                         gAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGcc                                                >  2:6179878/1‑93 (MQ=255)

                                         gAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGcc                                                >  1:3840883/1‑93 (MQ=255)

                                         gAAGCCAGGCACGCACCAACGCCCAGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGcc                                                >  1:6640636/1‑93 (MQ=255)

                                         aaaGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGcc                                                <  1:3641795/92‑1 (MQ=255)

                                          aaGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGccc                                               >  2:1657289/1‑93 (MQ=255)

                                           aGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGcccc                                              <  2:3022734/93‑1 (MQ=255)

                                            gCCAGGCGCGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGccccc                                             >  1:832599/1‑93 (MQ=255)

                                            gCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGccccc                                             >  1:6343086/1‑93 (MQ=255)

                                              cAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCTCCCAt                                           >  1:608697/1‑93 (MQ=255)

                                              cAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCAt                                           <  2:2217038/93‑1 (MQ=255)

                                              cAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCAt                                           <  2:1199493/93‑1 (MQ=255)

                                               aGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATg                                          >  1:5093954/1‑93 (MQ=255)

                                               aGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATg                                          <  2:1682845/93‑1 (MQ=255)

                                                     gcacCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCaaaa                                    <  2:6310204/93‑1 (MQ=255)

                                                        ccAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCa                                 >  2:1809608/1‑93 (MQ=255)

                                                           aCGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAAtc                              >  2:681515/1‑93 (MQ=255)

                                                           aCGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAAtc                              <  1:2848120/93‑1 (MQ=255)

                                                           aCGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAAtc                              >  2:473882/1‑93 (MQ=255)

                                                            cGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAAtct                             <  2:360342/93‑1 (MQ=255)

                                                            cGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAAtct                             <  1:4895585/93‑1 (MQ=255)

                                                              cccGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCa                           >  2:3108982/1‑93 (MQ=255)

                                                              cccGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCa                           >  1:5604539/1‑93 (MQ=255)

                                                               ccGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCaa                          <  1:4157611/93‑1 (MQ=255)

                                                                cGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAg                         >  2:2872353/1‑93 (MQ=255)

                                                                  gTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCg                       <  1:2712599/93‑1 (MQ=255)

                                                                    gAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTc                     <  1:4947965/93‑1 (MQ=255)

                                                                    gAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTc                     >  2:4982168/1‑93 (MQ=255)

                                                                       tGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCAcc                  >  2:2582223/1‑93 (MQ=255)

                                                                       tGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCCACGCCCCCATGCCAAAACCAATCTCAAGCGTCAcc                  <  2:6595018/93‑1 (MQ=255)

                                                                         cACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCgg                <  2:3878854/93‑1 (MQ=255)

                                                                          ccTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGgc               <  2:4479899/92‑1 (MQ=255)

                                                                           cTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGgcg              >  2:40579/1‑93 (MQ=255)

                                                                            ttCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGgcgc             >  2:4167330/1‑93 (MQ=255)

                                                                             tCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCt            <  1:2323404/93‑1 (MQ=255)

                                                                               aaTGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTc          <  2:5936832/93‑1 (MQ=255)

                                                                                aTGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCa         >  1:655232/1‑93 (MQ=255)

                                                                                   ccGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACgg      >  1:4699482/1‑93 (MQ=255)

                                                                                    cGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGc     >  2:1318271/1‑93 (MQ=255)

                                                                                     gaggagGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAACCCAATCTCAAGCGTCACCGGCGCTTCACGGcc    <  1:2234665/93‑1 (MQ=255)

                                                                                     gaggagGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGcc    >  1:6283834/1‑93 (MQ=255)

                                                                                     gaggagGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGcc    >  2:3668562/1‑93 (MQ=255)

                                                                                      aggaggTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCa   >  1:5109924/1‑93 (MQ=255)

                                                                                       ggAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCaa  <  1:4245784/93‑1 (MQ=255)

                                                                                       ggAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCaa  >  2:3164412/1‑93 (MQ=255)

                                                                                       ggAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCaa  >  1:5516496/1‑93 (MQ=255)

GACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAA  >  NC_000913/3100513‑3100692

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