BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,406,298	A→G	Q689Q (CAA→CAG)	b4179 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,406,298	0	A	G	99.7%	1136.0 / ‑5.1	352	Q689Q (CAA→CAG)	b4179	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (196/155); total (197/155)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

Only 100 of 366 total aligned reads displayed.

CTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAA  >  NC_000913/4406208‑4406387

cttcGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCg                                                                                         <  1:2522885/93‑1 (MQ=255)

 ttcGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAgcgc                                                                                        >  1:7865353/1‑93 (MQ=255)

 ttcGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAgcgc                                                                                        >  1:7764275/1‑93 (MQ=255)

     tccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGa                                                                                    <  2:1456946/93‑1 (MQ=255)

     tccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGCGCAGCGCCTGa                                                                                    >  1:157744/1‑93 (MQ=255)

     tccgtcTAGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGa                                                                                    <  2:100186/93‑1 (MQ=255)

      ccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGAt                                                                                   <  1:5466605/93‑1 (MQ=255)

      ccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGAt                                                                                   <  2:6970447/93‑1 (MQ=255)

       cgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATg                                                                                  <  2:5813316/93‑1 (MQ=255)

       cgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATg                                                                                  <  1:3523358/93‑1 (MQ=255)

       cgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATg                                                                                  <  2:7398215/93‑1 (MQ=255)

        gtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgg                                                                                 >  2:9293069/1‑93 (MQ=255)

        gtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgg                                                                                 >  2:1034705/1‑93 (MQ=255)

        gtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgg                                                                                 >  2:6500910/1‑93 (MQ=255)

         tcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATggg                                                                                <  1:1646084/93‑1 (MQ=255)

           tGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATggggg                                                                              >  1:6539761/1‑93 (MQ=255)

           tGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATggggg                                                                              <  2:4503511/93‑1 (MQ=255)

            ggACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgggggt                                                                             >  1:3748071/1‑92 (MQ=255)

                gacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATcc                                                                         <  2:6150317/93‑1 (MQ=255)

                gacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATcc                                                                         <  1:8016635/93‑1 (MQ=255)

                gacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATcc                                                                         >  2:4785000/1‑93 (MQ=255)

                  cTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAg                                                                       >  1:7087580/1‑93 (MQ=255)

                   ttGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGc                                                                      <  2:2346747/93‑1 (MQ=255)

                           tGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGAc                                                              >  1:1377995/1‑93 (MQ=255)

                           tGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGAc                                                              >  2:1808981/1‑93 (MQ=255)

                            gATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGAGCCTGATGGGGGAATCCAGCGGCCAGACt                                                             >  1:871576/1‑93 (MQ=255)

                              tggtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTa                                                           <  2:6352220/93‑1 (MQ=255)

                               ggtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTAt                                                          >  1:5832481/1‑93 (MQ=255)

                                gtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATc                                                         <  2:4069697/93‑1 (MQ=255)

                                gtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATc                                                         >  1:7782924/1‑93 (MQ=255)

                                 tctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATAGGGGAATCCAGCGGCCAGACTTATCg                                                        >  1:9687088/1‑93 (MQ=255)

                                   tggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGcc                                                      >  2:7917406/1‑93 (MQ=255)

                                    ggtcCATGTCTCTTCGCTGGACAATGACTACTAGCGCTTGGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCt                                                     <  1:4258662/93‑1 (MQ=255)

                                     gtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTg                                                    <  2:8860294/93‑1 (MQ=255)

                                      tcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTgg                                                   <  1:9379971/93‑1 (MQ=255)

                                         aTGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCg                                                >  1:5694115/1‑93 (MQ=255)

                                          tGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGa                                               >  1:271990/1‑93 (MQ=255)

                                          tGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGa                                               >  2:6099310/1‑93 (MQ=255)

                                           gTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGAt                                              <  2:9111883/93‑1 (MQ=255)

                                           gTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGAt                                              <  1:1002057/93‑1 (MQ=255)

                                               ctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcgcg                                          >  2:7201437/1‑93 (MQ=255)

                                               ctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcgcg                                          >  2:3921024/1‑93 (MQ=255)

                                                ttCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGt                                         >  1:1132074/1‑93 (MQ=255)

                                                ttCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGt                                         >  1:7280451/1‑93 (MQ=255)

                                                ttCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGt                                         >  1:6337735/1‑93 (MQ=255)

                                                   gCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGa                                      <  2:5914429/93‑1 (MQ=255)

                                                   gCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGa                                      >  2:3326392/1‑93 (MQ=255)

                                                    cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGag                                     >  1:3048951/1‑92 (MQ=255)

                                                    cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGaa                                     <  1:6109216/93‑1 (MQ=255)

                                                    cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGaa                                     <  2:2090477/93‑1 (MQ=255)

                                                      ggACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGt                                   >  1:5516836/1‑93 (MQ=255)

                                                       gACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGtt                                  <  1:5412793/93‑1 (MQ=255)

                                                       gACAATGACTACTATCACTTTGACCAGGTAGGGCAGCGCCTGGTGGGGGAATCCAGCGGCCAGACTTATCGCCTGGACGATCGAATGGGAGta                                  <  2:8986719/93‑2 (MQ=255)

                                                        aCAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTc                                 <  2:1645534/93‑1 (MQ=255)

                                                          aaTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcgc                               >  2:8491865/1‑93 (MQ=255)

                                                          aaTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcgc                               <  1:2764227/93‑1 (MQ=255)

                                                           aTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcgcg                              >  1:1277421/1‑93 (MQ=255)

                                                            tGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGt                             <  2:7666288/93‑1 (MQ=255)

                                                             gACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTc                            <  2:9566283/93‑1 (MQ=255)

                                                              actactATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCg                           >  2:6758340/1‑93 (MQ=255)

                                                              actactATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCg                           >  1:380805/1‑93 (MQ=255)

                                                               ctactaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGa                          >  2:5007153/1‑93 (MQ=255)

                                                               ctactaGCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGa                          <  2:511472/93‑1 (MQ=255)

                                                                tactaTCTCTTTGACCAGGGAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGaa                         <  2:8129839/93‑1 (MQ=255)

                                                                tactaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGaa                         >  2:1859518/1‑93 (MQ=255)

                                                                tactaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGaa                         >  2:1290267/1‑93 (MQ=255)

                                                                 actaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAg                        >  1:2137441/1‑93 (MQ=255)

                                                                  ctaTCGCTTTGACCAGGTATGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGc                       >  1:6977815/1‑93 (MQ=255)

                                                                  ctaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGc                       <  1:5996320/93‑1 (MQ=255)

                                                                   taTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGCAGCg                      >  1:4964420/1‑93 (MQ=255)

                                                                   taTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCg                      >  2:4639112/1‑93 (MQ=255)

                                                                   taTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCg                      >  2:2749139/1‑93 (MQ=255)

                                                                   taTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGCTCGCGTCGAAGCg                      >  2:687312/1‑93 (MQ=255)

                                                                    aTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCgg                     >  1:4741041/1‑93 (MQ=255)

                                                                        cTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTat                 >  1:1238655/1‑92 (MQ=255)

                                                                        cTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTaa                 >  1:7681860/1‑93 (MQ=255)

                                                                        cTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTaa                 >  2:4167593/1‑93 (MQ=255)

                                                                        cTTTGACCAGGTAGGGCAGCGCCTGATGGGAGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTaa                 <  2:6660906/93‑1 (MQ=255)

                                                                         tttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAat                <  2:3689556/93‑1 (MQ=255)

                                                                          ttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAata               <  2:2258413/93‑1 (MQ=255)

                                                                          ttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAata               <  1:4182909/93‑1 (MQ=255)

                                                                           tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat              <  1:8181320/93‑1 (MQ=255)

                                                                           tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat              <  1:9472339/93‑1 (MQ=255)

                                                                           tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat              <  2:345311/93‑1 (MQ=255)

                                                                           tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat              <  2:9352227/93‑1 (MQ=255)

                                                                           tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGGTCGCGTCGAAGCGGTTAatat              >  1:1941965/1‑93 (MQ=255)

                                                                              ccAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGa           <  1:75097/93‑1 (MQ=255)

                                                                              acAGGTAGGGCAGCGCCTGATGGCGGAATCCAGCGGCCAGACTTATCGCCTGTGGGACCGCGTGGAAGTTCGCTGCGAAGCGGTTAATATCGa           <  2:5509322/92‑1 (MQ=255)

                                                                               cAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGAc          <  1:6706155/93‑1 (MQ=255)

                                                                               cAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGAc          >  1:9197100/1‑93 (MQ=255)

                                                                                 ggTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGTGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa        >  2:3552993/1‑93 (MQ=255)

                                                                                 ggTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa        >  2:3172805/1‑93 (MQ=255)

                                                                                 ggTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa        <  1:3241928/93‑1 (MQ=255)

                                                                                  gTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAg       >  1:5111133/1‑93 (MQ=255)

                                                                                  gTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAg       >  2:2637162/1‑93 (MQ=255)

                                                                                   tAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgc      >  1:23610/1‑93 (MQ=255)

                                                                                   tAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgc      <  1:4378751/93‑1 (MQ=255)

                                                                                    aGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgcg     >  2:536273/1‑93 (MQ=255)

                                                                                     gggCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgcgc    >  2:827735/1‑93 (MQ=255)

                                                                                      ggCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCa   >  1:2216986/1‑93 (MQ=255)

CTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAA  >  NC_000913/4406208‑4406387

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